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ABSTRACT  OF  THE  DISSERTATION 
Shear,  Strain  and  Thermohaiine  Vertical  Fine  Structure  in  the  Upper  Ocean. 

by 

Steven  Paul  Anderson 
University  of  California,  San  Diego,  1992 
Professor  Robert  Pinkel,  Chairman 


Observations  of  fine  scale  (l-100m)  vertical  shear,  strain  and  thermohaiine 
structure  are  interpreted  in  terms  of  internal  wave  and  frontal  zone  dynamics.  The  data 
were  collected  with  the  Marine  Physical  Laboratory  161  kHz  Doppler  sonar  and  a  pair  of 
CTD's.  The  two  systems  were  operated  from  R/P  FLIP  for  18  days  during  February- 
March,  1990,  as  part  of  the  SWAPP  (35°  N,  127°  W)  experiment . 

These  simultaneous  measurements  of  velocity  and  density  provide  an  unique 
opportunity  to  study  the  vertical  advection  of  shear  layers.  The  observed  shear  variance 
is  dominated  by  near  inertial  waves,  with  35  m  vertical  wavelengths,  These  inertial 
waves  propagate  down  from  the  mixed  layer  and  are  vertically  advected  by  higher 
frequency  waves.  Vertical  wavenumber-frequency  spectra  of  shear  and  strain  are 
examined  in  Eulerian  (fixed  depth)  and  semi-Lagrangian  (isopycnal  following) 
coordinates.  The  semi-Lagrangian  shear  and  strain  variances  are  more  consistent  with 
variances  modeled  by  a  linear  superposition  of  internal  waves. 

The  p.  d.  f.’s  of  buoyancy  frequency  and  Froude  number  (buoyancy  normalized 
shear)  are  determined  using  tne  kinematic  strain  model  of  Pinkel  and  Anderson  (1992). 
The  Froude  number  distribution  shows  good  agreement  with  data  presented  here  and 
historical  data  (PATCHEX;  Kunze  et  al.,  1990). 

The  advection  of  a  frontal  zone  past  the  experiment  site  yielded  detailed 
observations  of  frontal  interleaving.  Two  vertical  interleaving  scales  are  observed.  Five 
to  fifteen  meter  scales  appear  to  be  double  diffusively  driven  and  are  consistent  with  the 
laboratory  model  of  Ruddick  and  Turner  (1979),  40  m  to  60  m  scales  are  likely 
determined  by  internal  wave  or  sub-mesoscale  frontal  dynamics.  The  observed  front 
appears  to  be  in  near  geostrophic  balance. 


Chapter  1 


Introduction 

This  dissertation  consists  of  three  investigations  of  fine  scale  (2-300m)  upper 
ocean  dynamics.  These  Investigations  are  presented  in  the  following  three  chapters. 
Each  chapter  ic  lar  ely  self  contained.  Chapter  2  examines  observations  of  vertical  shear 
and  strain  in  the  context  of  linear  internal  wave  dynamics.  Chapter  3  presents 
observations  of  thermohaline  fine  structure  associated  with  an  oceanic  frontal  zone. 
Chapter  4  examines  the  probability  distribution  of  shear  instabilities.  An  extensive  data 
set  unifies  t.,ese  investigations.  This  chapter  outlines  the  principal  results  from  each 
investigation. 

This  dissertation  presents  data  collected  primarily  by  two  systems  during 
February  and  March,  1990,  as  part  of  the  Surface  Wave  Processes  ProgTarr  (SWAPP; 
Appendix  A).  The  experiment  was  located  approximately  500  miles  west  of  central 
California  (35°  8.2'  N,  126°  59.0’  W;  Figure  A.l)  where  the  instrumentation  was 
deployed  from  R/P  FLIP  (Figure  1.1).  The  Marine  Physical  Laboratory  161  kHz  coded- 
pulse  Doppler  sonar  system  is  used  to  profile  the  ocean  currents  from  50  m  to  325  m 
depth.  The  processing  of  the  sonar  data  involved  several  steps  before  reliable  velocity 
profiles  where  obtained.  The  processing  (see  Appendix  E)  includes  variable  gain 
correction,  band  limited  noise  removal,  estimation  of  mechanical  alignment  error;  and  tilt 
and  heading  correction.  The  resulting  velocity  profiles  are  one  minute  averages  of  east 
and  north  components  of  horizontal  velocity  from  50  m  to  325  m.  The  acoustic  pulse 
length  smoothing  is  5.5  m  in  depth.  The  Marine  Physical  Laboratory  yo-yo  CTD  system 
profiled  the  ocean  density  and  thermohaline  structure  from  the  surface  to  420  m  every 
130  seconds.  A  pair  of  CTD’s  operated  in  tandem  with  the  upper  CTD  profiling  from  the 
surface  to  220  m  and  the  second  from  200  m  to  420  m.  The  CTD  data  processing  (see 
Appendix  C)  involves  phase  matching  the  temperature  and  conductivity  data  to  remove 
salinity  spiking,  filtering  the  pressure  signal  and  merging  the  data  from  both  CTD’s.  This 
results  in  profiles  of  temperature,  salinity  and  density  with  a  vertical  resolution  of  1.5  m. 
The  need  to  properly  align  the  depth  information  is  key  to  the  merging  to  the  sonar  and 
CTD  data  sets.  Depths  are  inferred  from  acoustic  travel  time  for  the  sonar  and  pressure 
for  the  CTD.  It  is  possible  to  cross  check  these  very  different  indicators  of  depth  using  a 
vertical  string  of  vector  measuring  current  meters  (VMCM)  provided  by  R.  Wellei.  The 
two  systems  were  operated  simultaneously  from  R/P  FLIP  for  12  days. 

L  Chapter  2,  observations  of  vertical  advection  of  fine  scale  shear  and  strain  are 
presented.  Tve  observed  shear  variance  is  dominated  by  near  inertial  waves  that 
propagate  down  from  the  mixed  layer  with  35  meter  vertical  wavelengths.  The  near 
inertial  shear  layers  are  vertically  advected  by  the  internal  semi-diurnal  tide  and  higher 


1 


frequency  waves  and  appear  to  ride  up  and  down  with  isopycnal  layers.  Vertical 
wavenumber-frequency  spectra  of  shear  and  strain  are  examined  in  Eulenan  (fixed  depth) 
and  semi-Lagrangian  (isopycnal  following)  coordinates.  This  is  the  first  analyses  of  the 
true  evolution  of  velocity  in  isopycnal  following  coordinates.  Strain  spectra  show-  energy 
propagating  towards  the  surface  in  the  semi-diurnal  tidal  band.  Eulerian  vertical 
wavenumber  spectra  display  a  vertical  wavenumber  dependence  which  changes 
significantly  with  frequency.  This  contrasts  the  spectral  models  of  Garrett  and  Munk 
(1972,  1975)  which  propose  a  vertical  wavenumber  dependence  that  is  independent  of 
frequency.  Interpretation  of  Eulerian  rotary  shear  is  confused  by  "Doppler"  shifting  of 
the  fine  scale  shear  layers  to  higher  frequencies.  The  semi-Lagrangian  spectra  of  shear 
and  strain  appear  more  separable  in  vertical  wavenumber  and  frequency.  The  comparison 
of  shear  and  strain  variances  exhibit  better  agreement  with  linear  theory  in  the  semi- 
Lagrangian  coordinates.  Estimates  of  semi-Lagrangian  counterclockwise  to  clockwise 
rotating  shear  variance  ratios  agree  well  with  linear  theory.  The  main  conclusion  is  that 
the  internal  wave  field  appears  more  consistent  with  a  linear  superposition  of  internal 
waves  when  observed  in  semi-Lagrangian  coordinates. 

Chapter  3  presents  observations  of  the  vertical  structure  of  temperature  and 
salinity  near  a  frontal  zone.  During  the  SWAPP  experiment,  an  oceanic  front  was 
advected  past  the  experiment  site.  The  front  was  characterized  below  200  m  by  warm, 
salty  water  encountered  at  the  Stan  of  the  experiment  and  a  transition  in  the  middle  of  the 
time  series  to  colder,  fresher  water  with  nearly  the  same  density.  The  temperature  and 
salinity  change  on  an  isopycnal  located  at  300  m  is  0.425  °C  and  0.0825  %c.  The  front 
appears  to  be  in  a  geostrophic  balance.  This  is  determined  by  comparing  velocity 
estimates  from  the  sonar  to  geopotential  anomalies  To  study  the  time  and  space  scales  of 
the  intrusive  interleaving  associated  with  the  front,  a  depth/time  series  of  temperature  and 
salinity  anomalies  is  presented.  The  anomalies  are  calculated  by  removing  the  mean 
temperature  and  salinity  at  constant  density  for  each  CTD  profile.  Two  vertical  scales  of 
interleaving  are  observed  near  the  front.  The  small  scale  intrusions  have  vertical  lengths 
of  5-15  meters,  are  double  diffusively  unstable  on  the  top  and  bottom  and  slope  upward 
from  the  warm  side  to  the  cold  side  of  the  front.  The  intrusions  are  consistent  with  scales 
predicted  by  Ruddick  and  Turner  (1979).  Larger  scale  interleaving  is  seen  at  40  to  60  m. 
These  larger  scales  are  may  driven  by  sub-mesoscale  or  internal  wave  forcing. 

In  Chapter  4,  the  probability  density  function  (p.  d.  f.)  of  buoyancy  normalized 
shear,  or  Froude  number,  is  examined.  It  is  conjectured  that  most  of  the  turbulence 
observed  in  the  ocean  interior  results  from  breaking  internal  waves.  Shear  instability  is 
the  hypothesized  cause  of  breaking.  The  Froude  number  is  a  ratio  of  the  inertial  forces  to 
the  buoyancy  forces.  It  is  equal  to  the  inverse  square  root  of  the  gradient  Richardson 
number.  The  potential  for  shear  instabilities  occurs  when  the  Froude  number  is  greater 
than  2.0  (Richardson  number  less  than  0.25).  The  model  p.  d.  f.  of  Froude  numbers  is 
determined  using  modeled  distributions  of  shear  and  strain.  The  buoyancy  frequency 
distribution  is  derived  from  the  kinematic  strain  model  of  Pinkel  and  Anderson  (1992). 
The  model  vertical  shear  distribution  is  given  by  a  Rayleigh  distribution  under  the 
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assumption  that  each  component  of  shear  is  independent  and  normally  distributed.  By 
assuming  that  shear  and  buoyancy  frequency  are  independent,  the  joint  p,  d.  f.  of  vertical 
shear  and  buoyancy  frequency  is  estimated.  Actual  data  from  the  CTD  and  sonar  suggest 
that  even  though  there  is  detectable  dependence,  this  is  a  reasonable  approximation.  The 
p.  d.  f.  of  Froude  numbers  is  calculated  from  the  joint  p.  d.  f.  of  vertical  shear  and 
buoyancy  frequency.  The  resulting  p.  d.  f.  successfully  models  histograms  of  Froude 
numbers  estimated  at  1  m  to  10  m  scales.  The  Froude  number  distribution  scales  with  the 
mean  buoyancy  frequency  normalized  by  the  mean  vertical  shear.  Little  change  in  the 
Froude  number  distribution  is  seen  from  125  m  to  300  m.  in  depth. 
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Figure  1.1  RIP  FLIP  Instrumentation  Diagram.  Upper  ocean  observations  are  made  with  the  Marine 
Physical  Laboratory  161  kHz  Doppler  sonar  system  and  a  pair  of  CTD’s.  This  diagram  shows  the 
deployment  configuration  used  during  SWAPP. 


Chapter  2 


Vertical  Advection  of  Near  Inertial  Shear  Lavers 

2.1.  Introduction  and  Background. 

Mixing  in  the  ocean  interior  is  thought  to  be  caused  by  the  fine  scale  vertical  shear 
and  strain  of  the  internal  wave  field.  Previous  Eulerian  estimates  of  internal  wave 
statistics  are  difficult  to  interpret  due  to  vertical  advection  (commonly  referred  to  as 
finestructure  contamination.)  Phillips  (1971)  describes  how  the  internal  wave  statistics  at 
a  fixed  depth  are  altered  by  advection  of  “frozen”  layers  past  an  instrument  at  fixed  depth. 
This  can  lead  to  “Doppler  shifting”  and  aliasing  of  intrinsic  wave  frequencies.  To 
examine  the  significance  of  internal  wave  self  advection,  shear  and  strain  are  studied  here 
in  semi-Lagrangian  and  Eulerian  coordinates.  Vertical  wavenumber-frequency  spectra 
are  estimated  and  compared.  Linear  theory  consistency  tests  of  counterclockwise  to 
clockwise  shear  variance  and  shear  to  strain  variance  ratios  are  studied  in  both 
coordinates.  The  results  show  that  the  semi-Lagrangian  statistics  give  better  agreement 
with  linear  internal  wave  theory  than  Eulerian  statistics. 

The  new  coordinate  system,  defined  by  isopycnal  surfaces,  is  referred  to  as  semi- 
Lagrangian.  Finestructure  contamination  can  be  minimized  by  tracking  isopycnal 
surfaces  in  depth  and  time  (Pinkel  et.  al.;  1991).  Fluid  dynamics  is  classically  formulated 
in  either  Eulerian  or  Lagrangian  variables.  Eulerian  variables  track  the  fluid  flow  from  a 
particular  point  in  three  coordinate  space.  Lagrangian  variables  divide  the  fluid  into 
macroscopically  small  parcels  that  are  typically  identified  by  their  initial  or  undisturbed 
positions  and  track  the  position  of  each  parcel.  Most  ocean  observations  and  model 
studies  (including  the  internal  wave  field)  are  made  in  terms  of  Eulerian  variables.  The 
semi-Lagrangian  statistics  are  not  true  Lagrangian  statistics  since  horizontal  advection  is 
neglected.  Much  of  the  Doppler  shifting  of  the  internal  wave  field  can  be  removed  in 
semi-Lagrangian  coordinates,  thus  yielding  a  better  estimate  of  the  intrinsic  wave 
frequencies. 

Previous  measurements  of  the  oceanic  internal  gravity  wave  field  have  shown 
common  spectral  form.  Garrett  and  Munk  (1972,  1975,  hereafter  GM)  developed  a 
model  of  the  internal  wave  variance  spectrum  by  consolidating  a  variety  of  oceanic 
temperature  and  velocity  measurements  into  a  universal  spectral  structure.  GM  internal 
wave  spectrum  is  an  empirical  model  that  assumes  that  the  velocity  and  temperature 
fluctuations  in  the  ocean  are  due  to  a  random  superposition  of  internal  waves.  The  GM 
spectrum  has  been  a  valuable  reference  for  comparing  various  internal  wave 
measurements.  It  has  had  modifications  for  high  wave  number  roll  off  (Munk,  1981  and 
Gregg  and  Kunze,  1991)  in  light  of  more  recent  observations.  The  GM  hypothesize  a 
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nearly  white  spectra  of  vertical  shear  and  strain  at  low  wave  numbers  and  a  roll  off  at 
vertical  wavenumbers  above  0.1  cpm.  GM  also  assume  that  at  all  internal  wave 
frequencies,  shear  and  strain  vertical  wave  number  spectra  have  the  same  shape. 
However,  not  all  observations  fit  the  hypothesized  spectral  form  and  the  model  lacks  a 
theoretical  explanation  for  its  spectral  shape. 

Recent  studies  suggest  that  self  advection  of  the  internal  wave  field  may  be 
significant,  thus  invalidating  the  assumption  that  the  internal  wave  field  is  a  linear 
superposition  of  waves.  Allen  and  Joseph  (1989)  model  the  oceanic  internal  wave  field 
as  a  superposition  of  linear  internal  waves  in  Lagrangian  coordinates.  By  examining  the 
spectra  in  Eulerian  coordinates  they  find  that  many  different  spectral  shapes  in 
Lagrangian  coordinates  will  result  in  the  same  Eulerian  spectra.  By  using  measured 
Eulerian  spectra,  the  can  not  determine  the  underlying  Lagrangian  spectrum.  Hines 
(1991a,  1991b)  has  used  Lagrangian  coordinates  to  model  the  atmospheric  internal  wave 
spectrum  with  similar  results.  Observations  of  vertical  wavenumber-frequency  shear 
spectra  by  Pinkel  (1985)  and  Sherman  and  Pinkel  (1991)  dispute  the  GM  assumption  of 
spectral  separability  of  wavenumber  and  frequency  in  Eulerian  coordinates.  They 
observe  red  vertical  wavenumber  dependence  of  shear  spectra  at  low  frequency  and  white 
wavenumber  dependence  at  high  frequency.  The  vertical  wavenumber-frequency 
spectrum  of  strain  in  semi-Lagrangian  coordinates,  as  estimated  by  Sherman  and  Pinkel 
(1991),  does  appear  to  be  more  separable.  The  wavenumber  dependence  is  nearly  white 
out  to  the  cutoff  at  all  frequencies.  Using  the  GM  model  spectrum  in  semi-Lagrangian 
coordinates,  Sherman  and  Pinkel  (1991)  show  that  vertical  advection  can  lead  to  the 
observed  Eulerian  shear  spectrum. 

Tests  for  consistency  of  shear  and  strain  variances  with  linear  theory  often  show 
more  shear  than  expected  at  high  frequencies  and  more  strain  than  expected  near  the 
inertial  band  (Eriksen,  1978;  Muller  et  al.,  1978;Kunze  et  al.,  1990).  Holloway  (1983) 
suggests  that  this  discrepancy  may  be  due  to  fine-scale  potential  vorticity  anomalies. 
Kunze  et  al.  (1990)  using  data  taken  with  a  neutrally  buoyant  float,  claim  that  the  excess 
shear  and  strain  is  due  to  Doppler  shifting  and  aliasing  of  internal  waves  with  intrinsic 
near  inertial  frequencies.  Sherman  (1989)  compares  observed  counterclockwise  to 
clockwise  rotating  shear  variances  to  linear  theory  and  finds  excess  counterclockwise 
shear  in  the  inertial  band.  Sherman  (1989)  also  reports  that  his  estimates  are  distorted  by 
vertical  advection  of  the  internal  waves. 

The  data  presented  in  this  chapter  were  taking  with  the  MPL  rapid  profiling  CTD 
system  and  the  161  kHz  Doppler  sonar  system.  The  161  kHz  Doppler  sonar  system  is 
used  to  study  the  ocean  currents  from  50m  to  325m  depth.  A  pair  of  rapid  profiling 
CTD's  provided  density  profiles  from  the  surface  to  420  meters.  The  two  systems  were 
operated  for  12  days  during  March,  1990  as  part  of  the  SWAPP  experiment.  The 
experiment  was  located  approximately  500  miles  west  of  central  California  (35°  8.2'  N, 
126°  59.0'  W  ;Figure  A.l).  The  simultaneous  measurement  of  velocity  and  density 
provides  a  unique  opportunity  to  study  the  vertical  advection  of  shear  layers.  Vertical 
wavenumber-frequency  spectra  of  shear  and  strain  are  examined  in  Eulerian  (fixed  depth) 


and  semi-Lagrangian  (isopycnal  following)  coordinates.  Shear  and  strain  statistics  are 
then  compared  to  linear  internal  wave  theory. 

2.2.  Stratification  and  Current  measurements. 

This  section  describes  the  merging  of  the  CTD  and  Doppler  sonar  data  into  one 
depth/time  series.  Two  profiling  CTD’s  measured  density  from  ~5  to  420  meters  in  depth 
with  a  vertical  resolution  of  1.5  meters  and  a  cycle  period  of  130  seconds  (Appedix  C).  A 
Doppler  sonar  system  measured  velocity  profiles  from  -40  to  320  meters  in  depth  with  a 
vertical  resolution  of  5.5  meters  every  60  seconds  (Appendix  E).  The  need  to  properly 
align  the  depth  information  is  key  to  the  merging  to  the  data  sets.  Depths  are  inferred 
from  acoustic  travel  time  for  the  sonar  and  pressure  for  the  CTD.  It  is  possible  to  cross 
check  these  very  different  indicators  of  depth  using  a  vertical  string  of  vector  measuring 
current  meters  (VMCM)  provided  by  R.  Weller. 

The  VMCM’s  were  deployed  from  the  aft  boom  of  FLIP:  10  meters  horizontally 
from  the  sonar,  which  was  mounted  on  the  thruster,  and  15  meters  from  the  upper  CTD 
on  the  starboard  boom.  The  VMCM  data  set  consists  of  there  15  minute  averaged 
temperature  and  vector  averaged  velocities  (depths  72.5  m  (vml2),  100.25  (vm!3)  and 
132.5  m  (vm20))  from  February  25,  1990  to  March  17,  1990.  Hour  averaged 
temperatures  from  vml2  and  vm20  are  compared  with  corresponding  hour  averaged 
temperature  profiles  from  the  upper  CTD  sampled  at  1.0  meter  increments.  The  CTD 
calibrated  current  meter  depths  are  determined  by  using  the  depths  of  maximum 
correlation  and  are  72±0.5  and  13410.5  meters  for  vml2  and  vm20  respectively  (Figure 
2.3c).  The  VMCM  calibrated  sonar  depths  are  found  by  comparing  15  minute  averaged 
sonar  velocity  profiles  with  the  corresponding  current  meter  records.  The  vertical  length 
of  each  acoustic  range  bin  is  2.91  meters  (Appendix  E).  Maximum  East/West  velocity 
correlations  occur  at  range  bin  numbers,  bn,  20  (73.2  m),  31  (105.2  m)  and  42  (137.2  m) 
for  vml2,  vml3  and  vm20  respectively  (Figure  2.3a).  Maximum  North/South  velocity 
correlations  occur  at  range  bin  numbers  21  (76.1  m) ,  31  (105.2  m)  and  42  (137.2  m)  for 
vml2,  vml3  and  vm20  respectively  (Figure  2.3b).  The  vml2  is  located  directly  between 
range  bins  20  and  21  and  the  first  depth  bin  is  located  at  15  meters  depth.  The  CTD 
calibrated  center  depth  of  each  range  bin  is  then  ( 12.1 +2.9 lbn)±2  meters  where  bn  is  the 
range  bin  number. 

Before  combining  the  velocity  and  density  profiles,  the  spatial  resolution  of  each 
density  profile  is  degraded  to  correspond  to  that  of  the  velocity  data.  This  is  done  using  a 
running  boxcar  filter  of  length  5.5  meters  in  depth.  Then  both  data  sets  are  sampled 
using  linear  interpolation  at  integer  meters  from  40  m  to  350  m.  These  are  low  pass 
filtered  using  a  5  pole  Butterworth  filter  with  a  cutoff  at  0.2  cpm  and  zero  phase  shift. 
The  velocity  data  are  then  smoothed  over  3  minutes  in  time  and  sampled  30  seconds  after 
the  start  of  the  CTD  profile  using  linear  interpolation.  This  results  in  the  combined 
velocity  and  density  data  which  is  sampled  every  1  meter  in  depth  from  40  to  350  meters 
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and  every  130  seconds  in  time.  The  time  when  data  was  obtained  on  both  systems  is 
from  yearday  63  through  year  day  75. 

The  Eulerian  depth/time  series  used  in  the  spectral  analysis  comes  from  the 
combined  velocity  and  density  data  from  yearday  63  to  75.  Both  velocity  and  density  are 
sampled  in  depth  every  2  meters  from  50  meters  to  306  meters  using  a  cubic  spline 
interpolation.  Three  profiles  are  averaged  and  missing  data  points  at  each  depth  are 
linearly  interpolated  in  time  to  create  a  profile  every  6.5  minutes.  To  create  the  Semi- 
Lagrangian  (isopycnal  following)  depth/time  series,  128  isopycnals  are  chosen  that  have 
a  mean  separation  of  2  meters  in  depth  from  50  meters  to  306  meters.  The  instantaneous 
isopycnal  depths  and  the  velocity  at  each  isopycnal  depth  are  found  by  linear 
interpolation  in  depth  for  each  profile.  The  Semi-Lagrangian  velocities  are  then  averaged 
over  3  profiles  and  missing  data  on  an  isopycnal  is  linearly  interpolated  in  time. 

2.3.  Strain  Spectra. 

Eulerian  and  semi-Lagrangian  strain  are  calculated  from  the  density  profiles 
obtained  from  the  profiling  CTD  system.  Eulerian  strain  is  obtained  by  first  differencing 
density  in  depth,  normalizing  each  profile  by  the  first  differenced  mean  density  profile 
and  subtracting  1.0  (Appendix  B).  Semi-Lagrangian  strain  is  calculated  by  first 
differencing  the  isopycnal  displacements  in  depth,  normalizing  by  the  mean  isopycnal 
separation  and  subtracting  1.0.  Spectral  variance  is  estimated  from  9  blocks  of  data,  each 
128  points  in  depth  and  512  points  in  time  corresponding  to  50-306  meters  depth  and 
55.4  hours  in  time.  The  blocks  are  overlapped  by  50%  in  time  and  tapered  in  space  and 
time  with  a  triangle  window  before  fourier  transforming. 

To  correct  high  wavenumber  variances  that  are  degraded  by  vertical  smoothing, 
the  spectra  are  multiplied  in  wavenumber  by  sinc2(5.5*7tm).  This  is  not  the  proper 
correction  for  the  semi-Lagrangian  strain  since  the  depth  smoothing  is  done  in  Eulerian 
coordinates.  It  is  applied  to  the  semi-Lagrangian  spectra  for  comparison  with  the 
Eulerian  spectra.  The  correction  is  approximately  the  same  if  straining  is  small  at  the 
scales  considered.  A  strain  vertical  wavenumber  spectrum  with  a  vertical  resolution  of 
1.5  meters  (Figure  2.5)  is  calculated  for  comparison  with  the.  vertically  smoothed 
spectrum.  The  high  resolution  spectrum  is  flat  from  low  wavenumbers  to  ~0. 1  cpm  The 
vertical  wavenumber,  cumulative  spectra,  corrected  at  high  wavenumber  shows  a  similar 
spectral  shape  out  to  0.1  cpm.  Thus,  the  correction  restores  the  true  spectra  shape. 

The  spectra  are  averaged  into  6  frequency  bands  and  plotted  versus  wavenumber 
out  to  0.1  cpm  (Figures  2.6  and  2.7).  The  Eulerian  and  semi-Lagrangian  strain  spectra 
look  very  similar  to  those  of  Sherman  and  Pinkel  (1991).  The  Eulerian  bands  do  not 
show  a  common  spectral  shape  but  do  converge  to  a  nr 2  slope  at  high  wavenumber.  The 
two  lowest  frequency  bands  are  red  at  all  observed  wavenumbers.  The  high  frequency 
bands  rise  from  low  wavenumber  then  roll  off  as  they  meet  the  common  high 
wavenumber  slope.  In  contrast,  the  semi-Lagrangian  spectrum  is  nearly  white  at  all  but 
the  highest  frequency  bands.  The  high  frequency  bands  rise  up  as  they  approach  0.1  cpm. 
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The  frequency  dependence  of  the  spectrum  can  be  seen  by  averaging  into  four 
vertical  wavenumber  bands.  Each  of  the  four  wavenumber  bands  is  smoothed  over  3 
frequency  bands  above  0.27  cph  and  12  frequency  bands  between  0.27  and  0.52  cph 
before  plotting  versus  frequency  (Figures  2.8  and  2.9) .  The  Eulerian  wavenumber  bands 
again  do  not  show  common  spectral  shape.  The  longest  wave  lengths  (128-32m)  are  red 
and  have  a  nr2  slope.  There  is  an  upward  energy  peak  in  the  128-64m  band  at  the  semi¬ 
diurnal  tide.  Pinkel  et  al.  (1987)  also  observed  net  upward  tidal  energy  off  the  coast  of 
California.  The  shortest  wavelengths  (32- 10m)  are  nearly  flat  at  low  frequency  with  less 
variance  than  the  longer  waves  and  roll  off  at  higher  frequency  with  more  variance.  No 
semi-diurnal  peak  is  present  in  the  short  wavelengths.  The  semi-Lagrangian  spectra  show 
congruent  form  at  all  wavenumbers.  The  longest  waves  resolved  (128-32m)  have  an 
upward  energy  semi-diurnal  peak.  The  cumulative  spectrum  has  a  nr1  slope  at  low 
frequencies  and  a  ox 2  slope  above  0.2  cph. 

2.4.  Shear  spectral  analysis. 

Eulerian  and  semi-Lagrangian  shear  are  calculated  from  the  velocity  profiles 
obtained  by  the  Doppler  sonar  system.  Eulerian  shear  is  obtained  by  first  differencing  the 
complex  horizontal  velocity  (u  =  u  +  iv)  at  2  meter  depth  intervals.  Semi-Lagrangian 
shear  is  calculated  by  first  differencing  the  complex  velocity  found  on  a  previously 
selected  set  of  isopycnals  that  are  spaced  2  meters  apart  on  average.  Spectral  variance  is 
estimated  from  9  blocks  of  data,  each  128  points  in  depth  and  512  points  in  time 
corresponding  to  50-306  meters  depth  and  55.4  hours  in  time.  The  blocks  are  overlapped 
by  50%  in  time  and  tapered  in  space  and  time  with  a  triangle  window  before  fourier 
transforming.  Rotary  spectra  are  found  by  averaging  the  spectra  in  each  quadrant  of 
wavenumber-frequency  space. 

The  high  wavenumber  variances  have  been  corrected  in  the  same  manner  as  the 
strain  spectra.  To  correct  high  wavenumber  variances  that  are  degraded  by  vertical 
smoothing,  the  spectra  are  multiplied  in  wavenumber  by  sinc2(5.5*7tm).  This  is  not  the 
proper  correction  for  the  semi-Lagrangian  strain  since  the  depth  smoothing  is  done  in 
Eulerian  coordinates.  It  is  applied  to  the  semi-Lagrangian  spectra,  however,  for 
comparison  to  the  Eulerian  shear  spectra  and  the  corrected  strain  spectra. 

The  shear  spectra  are  averaged  into  the  same  6  frequency  bands  as  the 
corresponding  strain  spectra  and  plotted  versus  wavenumber  out  to  0.1  cpm  (Figures 
2.10b  and  2.11b).  The  Eulerian  shear  spectrum  is  similar  to  those  of  Pinkel  (1985)  and 
Sherman  and  Pinkel  (1991).  The  highest  frequency  band  (1-0.5  hours)  is  red.  All  the 
bands  appear  to  have  a  slope  change  at  ~0.3  cpm.  The  lowest  four  frequency  bands  (55-4 
hours)  converge  at  high  wave  number  to  a  nr dependence.  All  frequency  bands  appear 
to  converge  near  0. 1  cpm.  The  semi-Lagrangian  shear  spectrum  has  some  of  the  same 
characteristics  as  that  of  the  modeled  spectrum  of  Sherman  and  Pinkel  (1991).  The 
highest  frequency  bands  (4-0.5  hours)  rise  slowly  towards  0.1  cph.  There  is  no 
convergence  of  the  frequency  bands  at  high  wavenumber.  The  lowest  frequency  bands 
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(55-8),  however,  still  show  a  peak  at  ~0.3  cpm  and  roll  off  at  higher  wavenumber.  As 
will  be  shown  later,  this  variance  peak  is  due  to  local  forcing  of  near  inertial  waves  at  the 
surface. 

To  examine  the  frequency  dependence,  spectra  are  averaged  into  4  vertical 
wavenumber  bands  and  smoothed  in  frequency  above  0.3  cph  the  same  as  the  strain 
spectra  (Figures  2.10a  and  2.11a).  The  lowest  wavenumbers  (128-32m)  in  the  Eulerian 
spectrum  rise  to  a  peak  at  the  inertial  frequency  {f=  .0478  cph).  Note  that  the  64-32m 
wave  length  band  has  a  second  peak  at  co  -  0.13  cph  which  is  the  sum  of  the  inertial  and 
semi-diurnal  tide  frequencies.  Kunze  et  al.  (1991)  also  report  a  vertical  shear  "super 
tidal"  peak  in  data  taken  with  the  RiNo  float  during  PATCHEX  which  took  place  1° 
south  of  SWAPP.  The  shorter  wave  lengths  have  no  inertial  peak  but  instead  are  nearly 
flat  until  rolling  off  at  higher  frequency  to  a  co  *2  slope.  This  is  the  same  general  shape 
observed  by  Sherman  and  Pinkel  (1991)  where  the  low  wavenumbers  have  a  <y*^  slope 
starting  at  the  inertial  frequency,  while  the  higher  wavenumbers  appear  more  flat  and 
have  higher  energy  density  at  high  frequency.  The  semi-Lagrangian  spectra  (Figure 
2.1  la)  have  quite  a  different  shape.  All  but  the  shortest  wavelengths  have  a  similar  shape 
and  show  an  inertial  peak.  The  spectrum  falls  with  a  or4  slope  until  a  slope  change  at 
-0.3  cph  and  the  64-32  meter  wavelength  band  approaches  a  or1  slope.  This  is  different 
from  the  modeled  spectra  of  Sherman  and  Pinkel  (1991)  who  observed  a  co  slope. 
There  is  no  longer  a  super  tidal  peak.  The  highest  wavenumbers  (16- 10m)  exhibit  no 
inertial  peak. 


2.5.  Comparison  with  GM76  model  variance. 


In  this  section,  measured  shear  and  strain  variances  are  compared  with  GM 
modeled  variances.  There  is  some  question  as  to  which  version  of  the  GM  model  to  use. 
Gregg  and  Kunze  (1991)  show  that  shear  and  strain  variances  obtained  by  using 
Munk(1981)  are  2/tc  times  smaller  than  those  obtained  from  GM76  (Garrett  and  Munk, 
1975  as  modified  by  Carins  and  Williams,  1976).  GM76  is  used  here  since  it  appears  to 
be  the  more  commonly  used  of  the  models.  From  Gregg  and  Kunze  (1991),  the  GM76 
spectra  of  vertical  strain,  fa,  and  vertical  shear,  fa,  as  functions  of  vertical  wavenumber, 
P  (given  in  rad/sec),  are 
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where  N0  =  .00524  rad  s'1  is  the  reference  buoyancy  frequency,  £  =  6.3  x  10'5  is  the 
dimensionless  energy  level  and  b  -  1300  m  is  the  scale  depth  of  the  thermocline.  The 
vertical  wavenumber,  p,  is  related  to  the  mode  number,),  by 
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and  for  j*  =  3  then  /?*  =  0.0073  ( N/N0 )  rad  m*1.  To  compare  the  GM76  variances  to  the 
measured  variance,  (2.1)  and  (2.2)  are  integrated  over  the  sampled  wavenumber 
bandwidth  yielding 


(2.4) 


(2.5) 


where  /3 /  and  pu  are  the  lower  and  upper  vertical  wavenumber  bandwidth  limits 
respectively.  The  GM76  shear  variance  to  strain  variance  ratio 
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is  independent  of  vertical  wavenumber  bandwidth. 


The  Semi-Lagrangian  and  Eulerian  variances  are  found  by  integrating  the 
wavenumber- frequency  spectra  over  vertical  wavelength  band  0.00391  to  0.1  cpm  (256  to 
10  meter  wavelengths)  and  frequency  band  0.0182  to  2.0  cph  (55  to  0.5  hours)  and  are 
displayed  in  Table  2.1.  Using  the  depth  average  N  =  4.4  cph  (Figure  2.12),  (2.4)  and 

(2.5)  over  estimate  the  both  Eulerian  and  semi-Lagrangian  total  strain  variance  and  under 
estimate  the  total  shear  variance.  Eulerian  and  semi-Lagrangian  strain  variances  are  85% 
and  75%  of  GM76  respectively.  Eulerian  and  semi-Lagrangian  shear  variances  are  130% 
and  1 16%  of  GM76.  The  excess  shear  and  low  strain  in  consistent  with  the  observations 
of  strong  near  inertial  waves.  The  expected  shear  to  strain  variance  ratio  from  linear 
theory  for  a  near  inertial  wave  is 


(*2) 

(A2)/V2 


&)2  +  /2 
co2~f2 


(2.7) 


Assuming  the  total  variance  is  dominated  by  a  single  near  inertial  wave,  the  wave 
frequency  can  be  estimated  from  the  variance  ratio  using  (2.6).  Using  N ,  the  Eulerian 
and  semi-Lagrangian  variance  ratio  can  be  accounted  for  by  a  near  inertial  wave 
frequency  of  1.26/ .  ( It  is  not  clear  whether  N  or  the  maximum  buoyancy  frequency. 
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Nm lt=6.2  cph.,  provides  the  appropriate  scaling  because  of  the  large  variation  in  the 
observed  buoyancy  frequency.  Using  N^,  the  frequency  becomes  1.61/which  is  slightly 
higher  than  the  observed  spectral  frequency  peaks.) 

2.6.  Rotary  shear  spectra!  analysis. 

Rotary  spectra  are  found  by  averaging  the  spectra  found  in  each  quadrant  of 
wavenumber-frequency  space.  The  spectra  are  divided  into  4  quadrants:  counter 
clockwise  downward  energy  (CCWDN)  ( +m,-(0 ),  clockwise  downward  energy  (CWDN) 
C -m,+co ),  counter  clockwise  upward  energy  (CCWUP)  (- m,-co )  and  clockwise  upward 
energy  (CWUP)  (+m,+co).  The  processing  and  averaging  are  the  same  as  discussed  with 
the  strain  and  shear  spectra.  The  Eulerian  rotary  shear  spectra  versus  vertical 
wavenumber  reveals  the  observed  slope  break  at  0.03  cpm  at  all  frequencies  in  the 
CWDN  and  CCWUP  spectra(Figure  2.13a-d).  This  slope  change  is  not  apparent  in  the 
CWUP  and  CCWDN  Eulerian  spectra.  The  semi-Lagrangian  spectra  only  have  a  peak  at 
0.03  cpm  in  the  lowest  wavenumber  bands  (55-8  hours). 

Vertical  wavenumber  is  averaged  into  4  bands  and  plotted  against  frequency  for 
the  Eulerian  and  semi-Lagrangian  spectra  (Figures  2.15  and  2.16).  The  super  tidal  peak 
is  only  located  in  the  64-32  m  band  of  the  Eulerian  CWDN  spectra.  At  high  frequency  in 
the  Eulerian  spectra,  the  highest  wavenumber  bands  have  the  largest  spectral  density.  An 
inertial  peak  is  observed  in  the  CWDN  and  CWUP  spectra  in  the  low  wavenumber  bands 
(128-32m).  The  highest  wavenumber  band  (16-10m)  is  white  at  low  frequency  and  rolls 
off  above  0.3-0.5  cph  in  all  four  quadrants.  The  wavenumber  bands  look  much  more 
similar  in  the  semi-Lagrangian  spectra.  They  have  aw1  power  law  at  all  frequencies  in 
the  CCWDN  and  CCWUP  spectra  with  the  lowest  wavenumbers  (128-64m)  having 
slightly  less  spectral  density  at  high  frequency.  An  inertial  peak  occurs  in  all  but  the 
highest  wavenumber  bands  in  the  CWDN  and  CWUP  spectra. 

The  total  shear  variance  is  dominated  by  near-inertial  waves  propagating  away 
from  the  surface.  The  CWDN  contains  43%  total  Eulerian  shear  variance  and  54%  of  the 
semi-Lagrangian  shear  variance  (Table  2.1).  The  prominent  wavelength  associated  with 
these  waves  is  35  meters.  There  is  noticeably  more  CW  semi-Lagrangian  than  Eulerian 
shear  variance.  Also  evident  is  the  decrease  in  CCWUP  semi-Lagrangian  variance  where 
vertical  advection  of  the  near  inertial  waves  has  the  largest  effect  (see  Section  2.8). 

2.7.  WKB  approximation. 

Linear  internal  wave  theory  assumes  that  the  background  density  gradient  is 
constant.  For  cases  where  the  buoyancy  frequency  varies  slowly  in  depth,  the  vertical 
wavenumber,  m,  also  becomes  a  slow  function  of  depth.  This  approximation  is  called  the 
Liouville-Green  or  WKB  approximation  (Gill,  1982).  As  a  wave  propagates  from  a 
region  of  large  to  small  N,  its  characteristic  vertical  wavelength  will  grow  and  horizontal 
velocity  amplitudes  will  shrink.  Following  the  approach  of  Gill  (1982),  assuming  that  co 
is  small  compared  with  N  and  that  the  vertical  scale  nr1  is  small  compared  with  the  scale 
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over  which  m  varies,  a  new  vertical  coordinate  can  be  found  so  that  m  is  constant  with 
depth.  The  vertically  stretched  vertics'  coordinate,  zm,  is  defined  by 
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where  Nm  is  the  maximum  value  of  the  buoyancy  frequency  and  N(z)  is  the  mean 
buoyancy  frequency  profile.  The  horizontal  velocity,  u,  scales  with  depth  according  to 
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where  um  is  the  WKB  vertical  velocity.  This  scaling  allows  the  comparison  of  wave 
energy  and  vertical  scales  measured  at  different  depths. 

The  WKB  scaling  is  applied  to  the  data  using  the  mean  buoyancy  profile  from  the 
combined  velocity  and  density  data  set  (Figure  2.12).  The  amplitude  of  the  velocity  at 
each  depth  is  adjusted  using  (2.9)  then  resampled  at  stretched  depths,  zm.  The  scaling  and 
stretching  is  applied  to  both  the  Eulerian  and  Semi-Lagrangian  velocity  field  data  sets. 
Using  the  same  routines  to  calculate  rotary  shear  spectra  as  described  in  Section  2.7 
yields  the  WKB  shear  spectra  (Figures  2.17  and  2.18).  The  high  wavenumber  correction 
has  been  applied  to  both  the  spectra. 

One  change  in  the  vertical  wavenumber  spectra  is  in  the  inertial  band.  The 
clockwise  downward  near  inertial  energy  from  the  non-WKB  scaled  Eulerian  and  Semi- 
Lagrangian  spectra  (Figures  2.13b  and  2.14b;  55-16  hour  period)  has  one  peak  at  0.009 
cpm  and  a  second  peak  at  0.03  cpm  In  contrast,  the  WKB  scaled  spectra  (Figures  2.17b 
and  2.18b;  55-16  hour  period)  have  only  one  peak  at  0.02  cpm.  This  difference  suggests 
that  the  near  inertial  waves  are  to  the  stretching  by  the  changing  buoyancy  frequency. 
(The  vertical  group  velocity  for  these  waves  is  30  meters  per  day 
(o)  =  1.26/)  and  it  would  take  nearly  10  days  for  this  near  inertial  wave  group  to  reach 
300  meters  from  the  surface.)  There  is  also  a  shape  difference  between  the  non-WKB 
and  WKB  spectra  at  high  wavenumbers.  This  is  expected  since  the  vertical  smoothing  is 
nonlinear  in  stretched  coordinates  and  has  not  been  correctly  removed.  The  non-WKB 
spectra  will  be  used  instead  since  there  are  no  significant  spectral  shape  changes  at  low 
wavenumber  and  the  effects  of  vertical  smoothing  are  better  understood. 

2.8.  Inertial  Wave  Vertical  Advection  Model. 

A  simple  model  is  presented  here  to  illustrate  inertial  wave  spectra  distortion  by 
vertical  advection.  Consider  the  complex  velocity  field  of  an  inertial  wave  in  vertically 
advective  coordinates  and  given  by 
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where  £  *s  defined  to  follow  vertically  advected  layers  (isopycnals),  Uq  is  the  velocity 
amplitude,  m  is  the  vertical  wavelength  and /is  the  inertial  frequency.  The  layer  vertical 
acceleration  is  found  from  the  total  derivative  of  the  vertical  velocity  field  on  that  layer. 
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The  term  of  the  left  hand  side  of  the  equation  is  the  vertical  acceleration  of  the  layer.  The 
first  term  on  the  right  side  is  the  local  Eulerian  acceleration  of  the  water  and  the  last  3 
terms  are  the  convective  accelerations.  Using  the  small  amplitude  wave  assumption,  the 
second  order  terms  are  dropped  and  the  equation  is  then  linearized.  Thus,  the 
acceleration  of  a  layer  is  approximated  by  the  acceleration  found  at  its  mean  depth  in 
Eulerian  coordinates.  Then  the  displacement  of  the  layer  is  found  by  integrating  the 
vertical  velocity  at  a  given  depth.  Using  the  mean  depth  of  a  layer  as  a  Lagrangian  label 
for  that  layer,  the  semi-Lagrangian  coordinate,  £,  for  a  single  internal  gravity  wave  is 
related  to  the  Eulerian  coordinates  by 

£(z,r)  =  z-/tsin(mz-<yr)  =  z-  T]  (2.12) 

where  T]  is  the  vertical  displacement,  m  and  d)  are  the  vertical  wavenumber  and 
frequency  of  the  displacement  wave  and  A  is  the  amplitude  of  displacement.  The 
Eulerian  shear  field  is  related  to  the  semi-Lagrangian  shear  by 

u,  =  —-uf  =  0-  7],)uf  (2. 13' 


If  the  strain  is  small  ( r\1  «  1),  solving  for  the  Eulerian  shear  yields 

uI(z,r)  =  m(//'8-/'VH) 


(2.14) 


If  the  vertical  advection  is  small  compared  to  the  vertical  wavelength  of  the  inertial  wave 
Mm'1  «  1)  then 

ui(z,r)  =  rnU>‘(ml'“,(l-fm7])  (2.15) 

Using  (2.12)  and  solving  for  the  Eulerian  vertical  shear  rotary  spectrum  yields 
0£(m,G>)  = 

m*  (2.16) 

Osx  +  — A2[<5(m  +  m  —  m)6(f  +  cb  +  6))  +  8(m  —  m  —  m)5(f  -  6)  +  <£>)] 
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where  the  semi-Lagrangian  or  intrinsic  wave  spectrum  is 

2 

0c,  (rn,cb)  —  — Ul5(m  ~m)8(co  +  d)) 

2  (2.17) 


The  Eulerian  spectrum  contains  “Doppler”  shifted  frequencies  and  wavenumbers  due  to 
the  vertical  advection.  The  variance  is  shifted  to  the  summed  and  differenced 
wavenumber  and  frequency  and  scales  with  the  displacement  variance  normalized  by  the 
inertial  vertical  wavelength. 

Assuming  that  the  vertical  wavelength  of  vertical  advection  is  long  compared  to 
the  inertial  wave  ( m»m ),  the  Eulerian  rotary  spectrum  will  be  altered  at  the  vertical 
wavenumber,  m,  in  both  the  CWDN  and  CCWUP  spectra  for  an  inertial  wave  traveling 
downward.  No  energy  from  this  wave  is  transferred  to  the  CWUP  and  CCWDN  spectra. 
The  observed  Eulerian  rotary  spectra  show  this  to  be  the  case  (Figure  2.13).  The  inertial 
shear  variance  is  dominated  by  m  =  0.03  cpm  seen  at  all  frequencies  in  Figures  2.13b  and 
2.13d  corresponding  to  the  CWDN  and  CCWUP  spectra.  The  distortion  is  absent  from 
the  CCWDN  and  CWUP  (Figures  2.13a  and  2.13c). 

The  frequency  dependence  of  the  distortion  will  rely  on  the  frequency  spectrum  of 
the  displacement,  0^.  Modeling  the  displacement,  r\,  as  a  spectrum  of  frequencies,  the 
Eulerian  shear  rotary  spectrum  is  written  as 

0£(m,d>)  =  0sx  +^-Uo2[0n(-/-d))  +  0n(/  +  ft))]5(m-m).  (2.18) 


Typical  0„  have  a  peak  at  the  semi-diurnal  tide,  M2,  and  fall  as  co  * 2  out  to  the  buoyancy 
cutoff  at  N  (Figure  2.19).  The  strongest  distortion  is  expected  at  co  =  -(/  +  M2)  and 
(0  =-(/ -  M2)  which  would  lead  to  a  super  tidal  and  sub-inertial  spectral  peaks.  The 
super  tidal  peak  is  observed  in  the  Eulerian  CWDN  shear  spectium  in  the  64-32  m  band 
(Figure  2.15b)  as  expected  and  not  anywhere  else.  The  sub-inertial  peak  can  not  be  seen 
because  the  records  are  too  short.  At  high  frequency  (co  » /),  the  distortion  will  have  a 
slope  of  <y  *2.  The  shear  frequency  spectrum  also  has  a  slope  of  nearly  co  * 2  suggesting 
that  observable  distortion  occurs  at  all  frequencies. 


2.9.  Counter  Clockwise  to  Clockwise  Shear  Variance  Ratio. 

The  CCW  to  CW  shear  variance  ratio  is  calculated  by  dividing  the  sum  of 
CCWDN  and  CCWUP  spectra  by  the  sum  of  CWDN  and  CWUP  spectra.  From  linear 
theory,  the  expected  CCW/CW  ratio  is 


0 

cc 

0. 


w 


/ 


\2 


for  co  >  f 


(2.19) 


15 


which  is  only  a  function  of  frequency.  As  co  approaches  the  inertial  frequency,/,  the 
waves  should  be  dominated  by  clockwise  rotation  (in  the  nonhem  hemisphere)  as  the 
ratio  approaches  zero.  At  higher  frequency,  rotation  becomes  less  imponant  and  CCW 
variance  should  equal  CW  variance. 

The  CCW/CW  ratio  is  computed  for  the  4  wavenumber  bands  and  smoothed  in 
frequency  as  discussed  for  both  Eulerian  and  semi-Lagrangian  shear  (Figure  2.20a-b). 
All  the  wavenumber  bands  in  the  semi-Lagrangian  coordinates  have  greater  CW  rotation 
near  co  =  f.  Only  the  lowest  wavenumbers  in  Eulerian  coordinates  appear  to  be  affected 
by  rotation  and  the  CCW/CW  ratio  for  the  high  wavenumbers  does  not  deviate 
significantly  from  1.0  at  any  frequency.  To  better  illustrate  the  deviations  from  linear 
theory,  the  CCW/CW  ratios  are  normalized  by  expected  values  (2.19)  and  plotted  with 
95%  confidence  limits  (Figure  2.21).  The  semi-Lagrangian  shear  fits  linear  theory  at  all 
wavenumbers  and  frequencies  better  than  the  Eulerian  shear.  At  high  frequency  (co  >  1 
cph),  CCW  shear  appears  to  be  overestimated  in  both  coordinates.  One  possibility  for 
this  discrepancy  is  measurement  noise  which  would  equally  corrupt  CW  and  CCW  shear 
variance  thus  pushing  the  ratio  closer  to  1.0  as  the  signal  to  noise  ratio  gets  small.  At  low 
frequencies,  the  CCW/CW  shear  ratio  measured  at  low  wavenumbers  is  better  predicted 
by  linear  theory  than  higher  wavenumbers.  This  also  could  be  caused  by  measurement 
noise  or  failure  to  prefectly  resolve  the  ~~mi-Lagrangian  coordinates  . 


2.10.  Strain  Variance  To  Shear  Variance  Ratio 


The  expected  strain  variance  to  shear  variance  ratio  from  linear  theory  is 
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where  is  the  strain  spectrum,  d>j  is  the  shear  spectrum,  N  is  the  buoyancy  frequency 
and /is  the  inertial  frequency.  Independent  of  wavenumber,  near  inertial  waves  have 
excess  shear  variance  and  high  frequency  waves  have  excess  strain  variance.  Observed 
ratios  are  found  in  Eulerian  and  semi-Lagrangian  coordinates  using  the  estimated  shear 
and  strain  spectra  from  four  wavenumber  bands  (Figure  2.22a-b).  The  Eulerian  ratio 
exhibits  excess  shear  at  the  two  middle  wavenumber  bands  (64- 16m)  except  at  the  inertial 
frequency.  This  excess  shear  has  been  observed  by  Sherman  (19S9)  and  Kunze  et  al. 
(1990).  They  both  claim  that  the  vertical  advection  of  low  mode  internal  waves  Doppler 
shifts  the  intrinsic  wave  frequency  to  higher  frequencies. 

The  semi-Lagrangian  strain  variance  to  shear  variance  ratio  fits  ratios  modeled  by 
linear  theory  very  well  from  (0  =/  to  (0  ~  0.5  cph  at  all  but  the  highest  wavenumber  band. 
Above  0.5  cph,  excess  shear  is  found  at  all  wavenumbers.  This  is  associated  with  the 
break  in  slope  of  the  semi-Lagrangian  shear  frequency  spectrum  (Figure  2.9a)  since  the 
semi-Lagrangian  strain  frequency  spectrum  has  no  observed  slope  change  (Figure  2.7 c). 
This  excess  shear  is  most  likely  due  to  possible  errors  in  estimating  the  semi-Lagrangian 
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shear.  First  of  all,  the  sonar  overestimates  the  horizontal  velocity  associated  with  high 
frequency  waves  due  to  beam  separation  biasing  (Appendix  D).  Secondly,  the  depth 
smoothing  of  the  density  profiles  has  removed  high  wavenumber  straining.  For  a  high 
strain  event,  the  estimated  isopycnals  depths  will  under  resolve  the  true  semi-Lagrangian 
coordinates.  Thirdly,  the  velocity  profiles  have  been  averaged  in  Eulerian  coordinates  for 
three  minutes  and  interpolated  to  a  single  profile  of  density.  This  will  cause  aliasing  of 
the  true  semi-Lagrangian  velocities.  By  not  fully  resolving  the  semi-Lagrangian 
coordinates  in  depth  and  time,  the  high  frequency  semi-Lagrangian  shear  is  still  corrupted 
by  vertical  advection  and  Doppler  shifting  of  the  intrinsic  wave  frequencies.  All  of  these 
may  lead  to  excess  shear  at  high  frequency. 

2.11.  Discussion  and  Conclusions. 

Vertical  advection  of  high  mode,  near  inertial  w-aves  by  higher  frequency  (co  > 
M2),  low  mode  waves  significantly  distorts  Eulerian  measurements  of  shear  and  strain. 
The  Eulerian  shear  to  strain  ratios  reveal  excess  shear  at  all  frequencies  above  the  near 
inertial  band.  A  simple  model  of  a  vertically  advected,  downward  propagating,  inertial 
wave  demonstrates  how  the  wave  energy  is  moved  to  higher  frequency  when  sampled  in 
fixed  depth  coordinates.  The  predicted  spectral  distortion  scales  with  the  variance  of  the 
displacement  normalized  by  the  vertical  wavelength  of  the  inertial  wave  squared  and  is 
found  at  the  sum  and  difference  inertial  frequency  and  advecting  frequencies.  A  super 
tidal  peak  in  the  Eulerian  shear  spectrum  is  seen  in  the  wavenumber  band  [0.0161- 
0.0313]  cpm.  The  displacement  variance  is  dominated  by  the  semi-diurnal  internal  tide 
and  this  super  tidal  peak  corresponds  to  the  sum  of  the  inertial  frequency  and  the  semi¬ 
diurnal  tide.  Both  observed  displacement  spectrum  and  the  Eulerian  shear  spectrum 
decay  in  frequency  with  a  slope  of  or2  above  the  super  tidal  peak  (Figure  2.25).  This 
indicates  distortion  of  the  shear  is  important  at  all  observed  frequencies  when  sampled  in 
Eulerian  coordinates. 

Semi-Lagrangian  shear  and  strain  ratios  closely  fit  those  predicted  by  linear 
theory  for  the  frequency  band  [f  -  0.5]  cph  and  vertical  wavelengths  [128-16]  cpm. 
Excess  shear  is  observed  above  0.5  cph  .  The  semi-Lagrangian  shear  frequency  spectrum 
exhibits  a  slope  change  from  or^-ar4  above/ to  or* -or2  above  0.5  cph  which  is  not  seen 
in  the  strain  spectrum.  Thus  the  excess  shear  may  due  to  Doppler  shifting  since  the  semi- 
I  .agrangian  coordinate  is  not  fully  resolved  in  time  and  space. 

The  total  strain  variance  is  75-85%  of  GM76  strain  variance  and  the  total  shear 
variance  is  1 16-130%  of  GM76  shear  variance.  The  deficiency  in  strain  and  excess  shear 
can  be  accounted  for  assuming  that  the  shear  and  strain  is  dominated  by  a  near  inertial 
wave  of  frequency  1.26/  Semi-Lagrangian  rotary  shear  spectrum  show  54%  of  the  shear 
variance  in  the  CW  downward  energy  corresponding  to  near  inertial  waves  propagating 
away  from  the  surface.  The  dominate  vertical  wavenumber  of  these  waves  is  0.03  cpm. 

Since  most  internal  w-ave  spectra  are  from  either  single  vertical  profiles  or  from  a 
moored  instrument  at  fixed  depth  yielding  a  time  series  it  is  useful  to  look  at  cumulative 
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spectra  for  comparison  (Figure  2.24).  The  cumulative  Eulerian  and  semi-Lagrangian 
shear  frequency  spectra([128-10jm)  have  some  distinct  differences.  The  semi- 
Lagrangian  spectrum  reveals  nearly  twice  energy  at  the  inertial  peak  and  much  less 
energy  at  all  frequencies  above  the  Eulerian  super  tidal  peak.  The  cumulative  Eulerian 
and  semi-Lagrangian  shear  vertical  wavenumber  spectra  ([55-0.5]  hours)  are  very  similar 
(Figure  2.24).  This  is  expected  since  the  distortion  by  vertical  advection  does  not  shift 
intrinsic  vertical  wavenumber.  The  vertical  wavenumber  shear  spectra  do  not  fit  the  form 
of  the  composite  shear  spectrum  presented  by  A.  E.  Gargett  et  al.  (1981).  The  spectra 
appear  more  like  those  of  Duda  and  Cox  (1989)  who  measured  vertical  shear  off  the  coast 
of  California  with  the  Cartesian  Diver  and  found  no  reproducible  power  law  in  the  band 
[0.02-0.1]  cpm.  The  spectra  rise  from  low  wavenumber  to  a  peak  at  0.03  cpm.  and  then 
decay  with  a  m  'l-$  power  law  out  to  0.1  cpm.  The  observed  spectral  shape  is  likely  due 
to  the  local  forcing  of  the  near  inertial  waves. 

Eulerian  shear  and  strain  vertical  wavenumber  and  frequency  spectra  show  similar 
characteristics  and  do  not  have  the  same  wavenumber  dependence  at  all  frequencies  as 
predicted  by  the  GM  model.  The  spectra  are  have  red  vertical  wavenumber  spectra  at  low 
frequency  and  blue  wavenumber  dependence  at  high  frequency  with  all  frequencies 
converging  near  0.1  cpm.  All  frequency  bands  in  the  Eulerian  shear  spectrum  show  a 
break  in  slope  at  0.03  cpm  corresponding  to  Doppler  shifting  of  a  wave  with  near  inertial 
intrinsic  frequency  to  higher  frequency.  The  semi-Lagrangian  shear  and  strain  vertical 
wavenumber  and  frequency  exhibit  more  separability  into  wavenumber  and  frequency. 
Both  spectra  have  slightly  blue  wavenumber  dependence  at  high  frequency,  [0.25-2]  cph. 
At  low  frequency,  [.0182-0.0625]  cph,  the  strain  spectrum  has  white  wavenumber 
dependence  and  the  shear  spectrum  are  slightly  red  with  a  peak  at  0.03  cpm. 

The  results  here  show  that  self  advection  of  the  internal  wave  field  is  significant 
on  vertical  scales  of  10  to  128  meters.  The  SWAPP  experiment  emphasis  this  fact  since 
the  dominant  vertical  scale  of  the  shear  variance  is  on  the  same  order  as  vertical 
displacement.  At  longer  wavelengths,  advection  should  not  be  as  important  and  Eulerian 
and  Lagrangian  statistics  should  be  the  same  (Allen  and  Joseph,  1989).  The  observation 
that  the  waves  obey  linear  internal  wave  dynamics  best  in  semi-Lagrangian  coordinate 
may  be  useful  in  developing  new  internal  wave  models. 
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Tune  on  year  day  65.  1990 


0:00  2.00  4  00  6  00  8:00  10  0?  12-00 
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East  Velocity  Shear  (sec  i 
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Figure  2.1  East  Velocity  Shear  on  yearday  6$.  East  velocity  shear  from  the  161  kHz  sonar  averaged 
over  3  minutes  is  contoured  to  a  depth  of  550  meters  from  0:00  to  12:00  :PST'  on  yearday  65.  1992. 
Overlaid  (black  curves!  are  tsopycnals  separated  by  10  meters  in  depth  on  average.  The  lack  of  acoustic 
scatterers  during  the  night  and  early  morning  limits  the  sonar  range  to  500  meters.  Alter  7:00.  reliable 
velocity  estimates  are  made  to  550  meters.  Strong  *hear  layers  are  visible  in  the  75  to  200  meiers  depth 
ranee.  Note  that  the  shear  lasers  appear  to  '  ride'  up  and  down  with  the  isopwnaK 
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Figure  2.2  Schematic  of  W.H.0.1.  aft  Current  Meter  Sting.  During  the  SWAPP  experiment,  W.  H, 
O.  I.  deployed  19  current  meters  on  2  strings  from  two  booms  on  FLIP.  The  port  deployment  extended 
to  41  meters  and  the  aft  deployment  (shown  in  this  schematic)  extended  to  132.25  meters  depth.  15 
minute  averaged  velocity  time  series  for  the  three  deepest  current  meters  (cm  28,  cm  38,  and  cm  18)  was 
obtained  from  W.  H.  O.  I.  for  comparison  with  the  M.  P.  L.  161  kHz  ADCP.  Also  shown  in  this 
diagram  is  the  depth  of  the  ADCP  which  was  mounted  on  FLIP’S  directional  thruster  and  has  a  6.5  meter 
range  resolution.  The  horizontal  spacing  between  the  current  meter  string  and  the  ADCP  was 
approximately  10  meters  and  is  not  drawn  to  scale  on  the  schematic. 
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Table  2.1  Shear  and  Strain  Variances.  This  table  shows  the  shear  and  strain  variance  found  by 
integrating  the  measured  shear  and  strain  spectra  over  the  vertical  wavelength  band  0.00391  to  0.1  cpm 
(256  to  10  meter  wavelengths)  and  frequency  band  0.0182  to  2.0  cph  (55  to  0.5  hours).  Variances  at  high 
wavenumber  have  been  amplified  to  compensate  for  roll  off  due  to  vertical  smoothing.  (Numbers  in 
parentheses  are  percent  of  total.)  <\2gm>  and  <S2gm>  are  found  by  integrating  the  GM76  spectrum 
over  the  same  wavenumber  band.  Shear  to  strain  variance  has  been  scaled  by  both  the  depth  averaged 
buoyancy  frequency  (N=4.4  cph)  and  the  maximum  buoyancy  frequency  (Nmax=6.2  cph)  for  comparison. 


vertical  wavenumber  (cpm)  vertical  wavenumber  (cpm) 


Figure  2.7  Semi-Lagrangian  Strain  Vertical  Wavenumber  Spectra.  Data  used  for  these  spectral 
estimates  are  from  the  combined  velocity  and  density  data.  The  spectra  have  been  corrected  at  high 
wavenumbers  for  vertical  smoothing,  (a),  (b),  and  (c)  show  variances  in  six  frequency  bands  plotted 
against  vertical  wavenumber,  (a)  contains  just  the  downward  propagating  energy,  (b)  contains  the  upward 
energy,  (c)  contains  the  total  energy,  (d)  is  the  cumulative  spectrum. 
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frequency  (cph)  frequency  (cph) 


Figure  2.8  Eulerian  Strain  Vertical  Frequency  Spectra.  Data  used  for  these  spectral  estimates  arc  from 
the  combined  velocity  and  density  data.  The  spectra  have  been  corrected  at  high  wavenumbers  for  vertical 
smoothing,  (a),  (b),  and  (c)  show  variances  in  four  wave  vertical  wavenumber  bands  plotted  against 
frequency,  (a)  contains  just  the  downward  propagating  energy,  (b)  contains  the  upward  energy,  (c) 
contains  the  total  energy,  (d)  is  the  cumulative  spectrum. 
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frequency  (cph)  frequency  (cph) 


Figure  2.9  Semi-Lagrangian  Strain  Vertical  Frequency  Spectra.  Data  used  for  these  spectral  estimates 
are  from  the  combined  velocity  and  density  data.  The  spectra  have  been  corrected  at  high  wavenumbers  for 
vertical  smoothing,  (a),  (b),  and  (c)  show  variances  in  four  wave  vertical  wavenumber  bands  plotted 
against  frequency,  (a)  contains  just  the  downward  propagating  energy,  (b)  contains  the  upward  energy, 
(c)  contains  the  total  energy,  (d)  is  the  cumulative  spectrum. 
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vertical  wavenumber  (cpm) 


IQ'3  10'2  10'1  10° 
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Figure  2.10  Eulerian  Shear  Vertical  Wavenumber-Frequency  Spectrum.  Euierain  shear  is  found  by 
first  differencing  the  velocity  profile  over  2  meter  depth  intervals,  (a)  the  average  spectra  density  is 
calculated  for  4  wavenumber  bands  and  plotted  against  frequency.  Spectra  have  been  smoothed  in 
frequency  above  0.27  c.p.h.  (b)  the  averaged  spectral  density  is  found  for  6  frequency  bands  and  plotted 
against  wavenumber.  >.c)  the  spectum  shown  in  3-D  format. 
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vertical  wavenumber  (cpm) 


Figure  2.11  Semi-Lagrangian  Shear  Vertical  Wavenumber— Frequency  Spectrum.  Semi-Lagrangian 
shear  is  found  by  first  differencing  in  depth  the  velocity  profile  found  on  isopycnals  separated  by  2  meters 
on  average,  (a)  the  average  specira  density  is  calculated  for  4  wavenumber  bands  and  plotted  against 
frequency.  Spectra  have  been  smoothed  in  frequency  above  0.27  c.p.h.  (b)  the  averaged  spectral  density  is 
found  for  6  frequency  bands  and  plotted  against  wavenumber,  (c)  the  spectum  shown  in  3-D  format. 
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Buoyancy  Frequency 


Figure  2.12  Mean  Buoyancy  Frequency  Profile.  Density  profiles  from  yearday  63 
averaged  together.  The  mean  buoyancy  profile  is  found  by  first  differencing  the  tim 
profile. 
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Figure  2.13  Eulerian  Rotary  Shear  Spectra  Vs.  Vertical  Wavenumber.  The  spectra  have  been  corrected 
at  high  wavenumbers  for  vertical  smoothing.  Six  frequency  bands  are  plotted,  (a)  contains  downward 
CCW  propaging  energy,  (b)  contains  downward  CW  propaging  energy,  (c)  contains  upward  CW 
propaging  energy,  (d)  contains  upward  CW  propaging  energy. 


32 


Lagrangian  Shear  (1/sty{cpm*cph)  Semi-Lagrangian  Shear  (l/sfyfcpm'cph) 


Figure  2.14  Semi-Lagrangian  Rotary  Shear  Spectra  Fs  Vertical  Wavenumber.  Same  as  Figure  2.13 
except  semi-Lagrangian  shear  spectra  are  plotted. 


frequency  (cph)  frequency  (cph) 


Figure  2.15  Eulerian  Rotary  Shear  Spectra  Vs.  Frequnecy.  The  spectra  have  been  corrected  at  high 
wavenumbers  for  vertical  smoothing.  Four  wavenumber  bands  are  plotted.  Approximate  95% 
confindence  limits  scale  with  frequency  smoothing,  (a)  contains  downward  CCW  propaging  energy,  (b) 
contains  downward  CW  propaging  energy,  (c)  contains  upward  CW  propaging  energy,  (d)  contains 
upward  CW  propaging  energy. 
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Semi-Lagrangian  Shear  (1/sJ)/(cpm’cph)  Semi-Lagrangian  Shear  (1/ss)/(cpm*cph) 


frequency  (cph) 


frequency  (cph) 


frequency  (cph)  frequency  (cph) 


Figure  2.16  Semi-Lagrangian  Rotary  Shear  Spectra  Ks.  Frequency.  Same  as  Figure  2.15  except  semi 
Lagrangian  shear  spectra  are  plotted. 


vertical  wavenumber  (cpm)  vertical  wavenumber  (cpm) 


Figure  2.17  WKB  Eulerian  Rotary  Shear  Spectra.  Velocity  profiles  are  stretched  and  amplitude 
normalized  according  to  WKB  theory  using  the  time  averaged  buoyancy  frequency  profile  prior  to  shear 
spectra  calculations  The  spectra  have  been  corrected  at  high  wavenumbers  for  vertical  smoothing.  Four 
wavenumber  bands  are  plotted.  Approximate  95%  confindence  limits  scale  with  frequency  smoothing, 
(a)  contains  downward  CCW  propaging  energy,  (b)  contains  downward  CW  propaging  energy,  (c) 
contains  upward  CW  propaging  energy,  (d)  contains  upward  CW  propaging  energy. 
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vertical  wavenumber  (cpm) 


vertical  wavenumber  (cpm) 


Figure  2.18  WKB  Semi-Lagrangian  Rotary  Shear  Spectra.  Same  as  Figure  2.17  except  semi 
Lagrangian  velocity  profiles  are  used. 


Frequency  (cph) 


Figure  2.19  Displacement  Frequency  Spectrum.  The  displacement  of  an  isopycnal  whose  mean  depth 
is  200m  is  tracked  for  1 5  days.  Shown  is  the  displacement  frequency  spectrum.  The  spectrum  has  been 


smoothed  in  frequency  above  1.0  c.p.h.  The  spectrum  has  a  peak  at  the  semi-diurnal  tidal  frequency,  has  a 
slope  of  co  down  to  the  buoyancy  frequency  cutoff  at  -  4.0  c.p.h. 
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Figure  2.20  Measured  CCW  to  CW  Shear  Variance  Ratio  .  The  ratio  of  counterclockwise  (CCW)  to  clockwise  (CW)  shear  variance  found 
in  the  rotary  shear  spectra  is  plotted.  The  heavy  line  is  predicted  CCW/CW  variance  ratio  using  linear  theory,  (a)  measured  Eulcrian 
CCW/CW  shear  ratio  averaged  over  4  wavenumber  bands,  (b)  measured  Scmi-Lagrangian  CCW/CW  shear  ratio  averaged  over  4  wavenumber 
bands. 
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Figure  2.21  Measured  CCW  to  CW  Shear  Ratio  Fit  to  Linear  Theory.  The  ratio  of  counterclockwise  (CCW)  to  clockwise  (CW)  shear 
variance  found  in  the  rotary  shear  spectra  normalized  by  that  expected  from  linear  theory  is  plotted,  (a)  normalized  Eulcrian  CCW/CW  shear 
ratio  averaged  over  4  wavenumber  bands,  (b)  normalized  Scmi-Lagrangian  CCW/CW  shear  ratio  averaged  over  4  wavenumber  bands. 
Approximate  95%  confidence  interval  for  the  longest  wavelength  band  is  drawn  around  1.0  and  tapers  towards  higher  frequencies  due  to 
frequency  smoothing. 
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Figure  2.22  Measured  Strain  to  Shear  Variance  Ratio.  The  ratio  of  measured  strain  variance  to  measured  shear  variance  normalized  by  the 
average  bouyancy  frequency  (N=6.15  cph).  The  heavy  line  is  predicted  strain  to  shear  ratio  using  linear  theory,  (a)  measured  Eulcrian  strain 
shear  ratio  averaged  over  4  wavenumber  bands,  (b)  measured  Scmi-Lagrangian  strain  to  shear  ratio  averaged  over  4  wavenumber  bands. 
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Figure  2.23  Noise  Corrected  Strain  to  Shear  Ratio  .  Same  as  Figure  2.22  except  that  an  assumed  white  noise  frequency  spcctium  has 
been  subtracted  from  each  shear  spectrum  wavenumber  band  before  calculating  the  strain  to  shear  ratio.  Removal  of  the  noise  decreases  the 
the  shear  variance  at  high  frequencies  (>0.4  eph)  and  measured  Semi-Lagrangian  strain  to  shear  ratio  is  more  closely  approximated  by  linear 
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Figure  2.24  Eulerian  and  Semi-Lagrangian  Cumulative  shear  spectra.  Vertical  wavenumber-frequency 
shear  spectra  are  integrated  from  low  frequency  to  successively  higher  frequency  and  plotted  for  (a) 
Eulerian  shear  and  (c)  semi-Lagrangian  shear.  Vertical  wavenumber-frequency  shear  spectra  are  integrated 
from  low  wavenumber  to  successively  higher  wavenumber  and  plotted  for  (b)  Eulerian  shear  and  (d)  semi- 
Lagrangian  shear. 
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Figure  2.25  Shear  Frequency  Spectrum.  Eulerian  and  semi-Lagrangian  shear  frequency  spectra  are 
found  for  the  vertical  wavenumber  band  [0.0078-0.1]  c.p.m.  The  almost  twice  the  variance  is  found  in 
the  semi-Lagrangian  near  inertial  band  than  in  the  Eulerian.  A  super  tidal  peak  is  observed  in  the  Eulerian 
spectrum  and  at  higher  frequencies  more  variance  is  found  in  the  Eulerian  spectrum.  Vertical  advection 
has  Doppler  shifted  waves  with  near  inertial  intrinsic  frequency  leading  to  contamination  of  the  Eulerian 
frequency  spectrum. 
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Chapter  3 


Thermohaline  Finestructure 


3.1.  Introduction 

In  this  chapter  the  fine  scale  vertical  structure  of  temperature  and  salinity 
associated  with  an  oceanic  frontal  zone  is  examined.  An  oceanic  front  is  a  boundary 
between  two  water  masses  with  different  thermohaline  characteristics.  Fronts  play  an 
important  roll  in  the  transformation  of  mesoscale  inhomogeneities  into  thermohaline 
finestructure.  Intrusive  layers  (a  layer  of  one  water  type  surrounded  by  water  with 
different  thermohaline  characteristics)  as  well  as  temperature  and  salinity  inversions  are 
commonly  observed  in  thermohaline  fine  structure  near  oceanic  fronts  (Home  1978; 
Georgi,  1978;  Joyce  et  al.  1978;  Posmentier  and  Houghton,  1978,  Fedorov,  l98t>,.  The 
formation  of  finestructure  is  an  important  process  since  the  interleaving  of  the  two  water 
masses  acts  as  stirring  which  enlarges  the  surface  area  between  the  two  water  masses  and 
thus  enhances  cross  frontal  fluxes. 

There  are  two  driving  mechanisms  for  interleaving  near  a  frontal  zone. 
Horizontal  advection  by  internal  waves  may  move  a  water  parcel  from  one  side  of  the 
front  to  the  other.  If  horizontal  temperature  gradients  are  strong  enough,  this  may  lead 
the  salinity  or  temperature  inversion  and  double  diffusive  instabilities.  Turner  (1978) 
observed  in  laboratory  models  that  double  diffusion  alone  can  drive  interleaving  at  the 
interface  of  two  water  masses  with  different  thermohaline  characteristics. 

The  advection  of  a  frontal  zone  past  the  SWAPP  site  during  the  operation  of  the 
CTD  system  provides  a  unique  opportunity  to  study  the  finestructure  associated  with  the 
front.  The  MPL  rapid  profiling  CTD  system  profiles  from  the  surface  to  420  m  every’  130 
seconds  during  yeardays  60  to  75, 1990  (Appendix  C).  The  data  are  examined  here  in  the 
context  of  traditional  fine  scale  frontal  analysis.  The  data  are  lacking  ir.  that  they  are  only 
a  time  series,  not  a  survey,  and  spatial  gradients  are  not  directly  known.  They  are 
estimated  by  time  integrating  velocity  data  from  the  MPL  161  kHz  Doppler  sonar  system. 
The  time  resolution  of  the  CTD  data  allows  the  tracking  of  intrusive  layers  in  isopycnal 
coordinates.  The  front  is  examined  for  a  geostrophic  balance  by  comparing  geopotential 
anomalies  to  velocity  shears  from  the  Doppler  sonar.  Intrusions  are  tracked  by  removing 
the  mean  temperature  and  salinity  on  isopycnais  yielding  a  depth/time  series  of 
thermohaline  anomalies. 


45  1 


3.2.  Description  of  Data 


An  XBT  survey  of  the  waters  surrounding  the  SWAPP  experiment  site,  taken  on 
yearday  55,  1992  (see  Appendix  A),  reveals  a  front  in  the  200  to  400  meter  depth  range. 
The  front  is  initially  located  to  the  northwest  of  FLIP’S  location  and  is  characterized  by 
warm  water  to  the  southeast  and  cold  water  to  the  northwest.  This  is  the  only  spatial 
survey  available  of  the  deep  frontal  structure.  Sea  surface  temperature  estimates  from 
satellite  data  were  limited  by  cloud  cover  and  the  surface  water  temperatures  do  not 
appear  correlated  with  the  deep  water.  The  SWAPP  experiment  was  not  originally 
designed  to  be  a  study  of  local  hydrographic  features.  However,  the  data  from  the  lower 
of  the  two  CTD’s  in  operation  provides  a  detailed  observation  of  frontal  structures  as  the 
water  masses  are  advected  past  the  experiment  site.  The  rapid  profiling  CTD  system 
provides  a  time  series  consisting  of  9200  profiles  from  200  to  420  meters  with  a  cycle 
period  of  130  seconds  and  was  in  operation  from  yearday  60  to  75. 

The  CTD  profiles  clearly  show  two  different  water  masses.  A  scatter  plot  of 
temperature  and  salinity  (Figure  3.1)  from  200  to  420  m  reveals  a  clustering  of  data 
points  along  two  different  characteristic  T-S  lines.  Few  data  points  are  found  to  lie  near 
the  cruise  averaged  T  -  S  curve.  A  picture  of  the  frontal  structure  can  be  formed  by 
tracking  the  temperature  along  an  isopycnal  in  depth  and  time  (Figures  3.2  and  3.3).  It  is 
useful  at  this  point  to  divide  the  time  series  into  three  segments.  The  first  segment 
contains  yeardays  60  to  66  (drop  numbers  0  to  3750).  It  is  characterized  by  warmer  and 
saltier  water  than  average.  The  water  mass  is  relatively  free  from  intrusions  except  for  a 
small  intrusion  of  cold  water  that  is  found  below  350  meters  from  day  60  to  63.  The 
temperature  (and  salinity)  along  all  the  isopycnals  slowly  decreases  in  time.  The  second 
segment  contains  yeardays  66  to  70  (drop  numbers  3750  to  6250).  It  is  characterized  by  a 
rapid  drop  in  temperature  (and  salinity)  along  isopycnals.  The  deepest  isopycnals 
(420>  z  >300  m)  continually  decrease  in  temperature  during  this  time  while  the  shallower 
isopycnals  (300>  z  >200  m)  go  from  cold  to  warm  to  cold.  Many  temperature  inversions 
occur  at  all  depth  ranges.  The  third  segment  contains  yeardays  70  to  75  (drop  numbers 
6250  to  9500)  and  is  characterized  by  colder  and  fresher  water  than  average  at  all  depth 
ranges.  Many  temperature  inversions  can  still  be  seen.  The  temperature  change  along  an 
isopycnal  increases  with  depth  and  thus  the  average  vertical  gradients  of  temperature  and 
salinity  are  larger  at  end  of  the  experiment  than  at  the  beginning. 

Estimating  the  horizontal  gradients  of  temperature  and  salinity  across  the  front 
can  be  done  it  two  ways.  The  XBT  survey  revealed  a  horizontal  temperature  gradient  in 
the  200-400  meter  depth  range  (Appendix  A).  The  estimated  cross  frontal  temperature 
gradient  is  1.85xl0*5  °C/m.  Assuming  that  the  cross  frontal  density  gradient  is  small,  the 
salinity  gradient  can  be  estimated  by 


dS  __  a  dT 
dy  P  dy 


(3.1) 
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where  /3  =  p 

y°  dS 

and  y  is  the  across  front  coordinate,  a  and  (5  are  nearly  constant  from  200  to  420  meter 
and  are  equal  to  a  =  1.508xl0'4  0C_i  and  p  =  7.712xl0"4  %o’1.  Using  (3.1),  estimated 
salinity  gradient  is  3.62xl0*6  %o/m.  This  method  may  underestimate  the  horizontal 
gradients  due  to  the  low  spatial  resolution  of  the  survey. 

The  second  method  is  to  compare  a  time  series  of  temperature  on  an  isopycnal 
with  the  time  integrated  velocity  field.  The  temperature  and  salinity  on  three  isopycnals 
with  average  depths  of  250,  300  and  350  m  are  tracked  (Figure  3.3).  Between  yeardays 
66  and  68,  the  temperature  and  salinity  decrease  by  approximately  0.425  °C  and  0.0825 
%o  on  the  300m  isopycnal.  The  largest  temperature  change  occurs  on  the  deepest 
isopycnal  which  continues  to  decrease  until  yearday  70.  The  vertical  temperature 
gradient  becomes  larger  on  the  cold  side  of  the  front.  The  horizontal  range  of 
temperatures  measured  is  still  not  as  large  as  that  seen  in  depth  averaged  temperatures 
from  the  XBT  survey  due  to  the  temperature  gradient  change.  The  horizontal 
displacement  is  found  by  time  integrating  the  depth  averaged  velocity  from  the  161  kHz 
sonar  (Figure  3.4).  The  resulting  trajectory  assumes  the  velocities  result  from  a  frozen 
ocean  moving  past  an  anchored  FLIP.  The  depth  averaged  velocities  are  dominated  by  a 
current  of  nearly  9  cm/s  flowing  towards  the  east.  Using  the  displacement  from  yearday 
66.7  to  67.2  (-7.1  km)  yields  temperature  and  salinity  gradients  of  5.27xlO-5  °C/m  and 
1.03  xlO*5  %c/ m.  It  is  unlikely  that  the  velocity  field  is  perpendicular  to  the  from,  so  the 
width  of  the  front  is  overestimated.  Even  though  these  cross  frontal  gradients  are  larger 
than  those  estimated  from  the  XBT  survey,  they  may  still  be  smaller  than  the  true  cross 
frontal  gradients. 

3.3.  Geopotential  Anomalies 

Often  oceanic  fronts  are  characterized  by  a  horizontal  density  gradient  across  the 
front  and  an  associated  geostrophic  velocity  balance.  Since  no  spatial  survey  was 
conducted  of  the  potential  density  field,  the  only  clue  of  a  geostrophic  balance  comes 
from  a  time  series  of  geopotential  anomaly,  <E>’ .  The  geopotential  anomaly  at  some 
pressure,  p,  is  defined  by 


-i  dp 
a-p  — 
BT 


(3.2) 


where  the  specific  volume  anomaly,  5,  is  defined  in  terms  of  the  specific  volume,  vs, 

5  =  vt(S,T,p)-vJ(35,0,p).  (3.3) 

Geopotential  anomalies  over  four  100  meter  depth  ranges  ([0-100],  [100-200],  [200-300], 
[300-400])  are  found  using  hour  averaged  density  profiles  (Figure  3.5).  The  signal  at  all 
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depth  ranges  is  dominated  by  the  semi-diurnal  internal  tide  vertically  advecting  the  mean 
density  profile.  This  tidal  signal  can  be  removed  by  time  averaging  O'  Geopotential 
anomalies  in  the  0  to  100  meter  depth  ranges  are  much  toO  large  to  be  accounted  for  by  a 
geostrophic  balance  and  are  likely  due  to  local  mixed  layer  dynamics.  The  [100-200]m 
water  gets  steadily  more  dense  with  time,  the  [200-300]m  water  does  not  have  any  net 
change  in  density  across  the  front  and  the  [300-400]m  water  gets  steadily  more  dense 
with  time.  Thus  the  cold,  fresh  water  encountered  at  the  end  of  the  experiment  is  more 
dense  than  the  warm,  salty  water  found  at  the  same  depths  at  the  beginning  of  the  cruise. 


It  is  useful  to  estimate  the  geostrophic  flow  needed  to  support  such  a  horizontal 
pressure  gradient.  The  geostrophic  velocity  at  pressure  p\  relative  to  the  velocity  at  p2  is 
calculated  by 


v,  (A )  ~  v,  (Pi )  =  /"'  “  {4>'(A )  “  $'(p2 )} , 

(3.4) 

ut  <A )  ~  ( A )  -  ’  37  {$'(  A )  -  O \p2 )} . 

dy 

(3.5) 

were  x  and  y 
components, 
estimated  by 

are  the  east  and  north  coordinates  and  u  and  v  are  the  east  and  north  velocity 
Using  the  time  series  of  O’ ,  the  magnitude  of  the  relative  velocity  is 

!V4G0  V30o|  ~  [^^400  (^2 )  ^300(^2))  (O400(t,)  O300(r,))] 

(3.6) 

where  the  AX  is  the  time  integrated  current  trajectory  determined  by 

'2 

-J  udt 

AX= 

\t2  ~  h ) 

(3.7) 

The  depth  average  velocity  vector,  u,  is  determined  using  velocity  profiles  obtained  from 
the  MPL  161  kHz  sonar  (Appendix  E).  Velocity  estimates  are  averaged  over  6.5  minutes 
in  time  and  100  to  350  meters  in  depth  and  integrated  in  time  to  get  AX  (Figure  3.4).  The 
velocities  are  dominated  by  an  eastward  flowing  current  of  nearly  9  cm/sec  from  yearday 
63  to  yearday  73.  The  northward  component  of  velocity  changes  direction  on  yearday 
68.  The  total  estimated  displacement  from  yearday  63  to  73  is  85  km  to  the  west  with 
negligible  southward  displacement.  The  net  change  in  geopotential  anomaly  from 
yearday  63  to  73  is  approximately  0.120  m2/s2  for  [O200-O100L  ~0  0  m2/s2  for  [O300- 
O200]  and  -0.024  m2/s2  for  [0400-0300]-  Since  the  total  displacement  is  towards  the 
west,  the  northward  component  of  the  current  is  calculated  using  (3.6).  The  geostrophic 
velocity  at  100  m  relative  to  200  m  is  1.7  cm/s  northward,  at  200  m  relative  to  300  m  is 
0.0  cm/s  northward  and  at  300  m  relative  to  400  m  is  0.4  cm/s  southward.  For 
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comparison,  the  northward  component  of  sonar  estimated  currents  at  100  m,  200  m,  and 
300  m  are  averaged  from  yearday  63  to  73.  The  mean  velocity  difference  is  2.03  cm/s 
between  100  and  200  m  and  0.03  cm/s  between  200  and  300  m.  These  are  consistent 
with  the  estimated  geostrophic  velocity  calculations.  Fedorov  (1986)  says  that  small- 
scale  fronts  are  likely  to  be  exceptions  to  geostrophic  balance.  Fedorov  refers  to  this  as 
“semi-geostrophic”  phenomena  and  is  the  result  of  mesoscale  forcing  leading  to  the  front 
f  'rmation  (frontogenesis).  Nothing  is  known  about  the  history  of  this  front  but  the 
rjrthward  component  of  velocity  appears  to  be  in  geostrophic  balance.  Without  a 
north/south  section,  it  can  not  be  determined  if  the  strong  eastward  current  observed  is  in 
a  geostrophic  balance. 


3.4.  Vertical  Interleaving 


Typical  finestructure  thermohaline  profiles  from  the  lower  CTD  reveal  many 
temperature  and  salinity  inversions  on  scales  of  3  to  20  meters  (Figure  3.6).  Temperature 
and  salinity  inversions  are  matched  such  that  the  resulting  density  profile  is 
hydrostatically  stable.  The  energy  stored  in  the  inversions  is  released  through  double 
diffusive  mixing  processes. 


The  difference  in  the  molecular  diffusivity  between  salt  and  heat  in  the  ocean 
which  has  a  significant  effect  on  vertical  fluxes  when  both  salinity  and  temperature 
increase  or  decrease  with  depth.  The  energy  driving  these  fluxes  comes  from  the 
thermohaline  property  that  has  a  destabilizing  effect  on  the  density  gradient.  The 
diffusive  case  occurs  when  temperature  and  salinity  increase  with  depth.  It  is 
characterized  by  convective  layers  and  by  more  heat  flux  than  salt  flux  through  the  layer. 
The  salt  fingering  case  occurs  when  temperature  and  salinity  decrease  with  depth.  It  is 
characterized  by  intrusive  ‘fingers’  propagating  vertically  across  the  interface  with 
vertical  salt  fluxes  dominating  heat  fluxes.  The  potential  for  one  of  these  processes  to 
take  place  is  quantified  by  the  density  ratio 


(3.8) 


or  the  Turner  angle 


Tu  =  tan'! 


far,-/*?,) 

W.+PS,) 


=  tan-1(/?p)-45°. 


(3.9) 


where  T  is  the  temperature,  S  is  the  salinity  and  z  is  the  vertical  coordinate.  The  two 
cases  are  then  described  by 
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0<flp<l 

-90°  <  Tu<  -45° 
i  *Rp 

45°  <  Tu  <  90°  I 


Diffusive 


Fingers. 


A  typical  profile  of  Tu  reveals  that  much  of  the  water  column  (200-420m)  is  subject  to 
double  diffusive  processes  (Figure  3.7). 

The  percent  of  the  water  column  that  is  subject  to  salt  fingering  (Tu  >  45°)  is 
typically  twice  as  large  as  the  percentage  of  the  water  column  that  is  subject  to  diffusive 
convection  (Tu  >  45°)  (Figure  3.8).  During  the  first  part  of  the  time  series  (yeardays  60 
to  65),  diffusively  unstable  water  occurs  less  than  3%  of  the  time,  whereas  fingering 
occurs  in  10%  to  20%  of  the  time.  On  yearday  65  occurrences  of  temperature  and 
salinity  inversion  increases  dramatically.  This  increase  takes  place  before  the  large 
isopycnal  temperature  drop  that  takes  place  at  yeardays  66  to  67  (Figure  3.3).  The 
incidents  of  temperature  and  salinity  inversions  reach  a  maximum  at  the  time  of  the  large 
temperature  drop  (yeardays  66-67)  and  again  at  yearday  69  when  the  isopycnals  in  the 
200  to  300  meter  depth  range  warm  up  for  a  half  a  day.  As  much  as  50%  of  the  water 
column  from  200  to  420  meters  is  double  diffusively  unstable.  After  yearday  70,  the 
percentage  of  temperature  inversions  decreases  to  5%  then  slowly  rises  to  10%  at  yearday 
75  and  the  percentage  of  salinity  inversions  stays  large  at  30%  to  35%. 

The  likely  cause  for  the  large  number  of  temperature  and  salinity  inversions  cross 
frontal  advection.  With  the  small  vertical  mixing  rates  typically  found  in  the  ocean 
interior,  these  intrusions  can  carry  some  of  the  thermohaline  characteristics  of  their 
source  waters  for  a  long  time.  This  may  be  why  the  number  of  inversions  rises  before  the 
actual  front  is  encountered.  There  are  two  kinematic  effects  which  would  cause  the 
horizontal  displacement  of  water  parcels  across  a  front.  The  first  is  horizontal  velocities 
from  the  internal  wave  field.  The  second  is  horizontal  intrusions  or  noses  that  are  driven 
by  the  vertical  flux  of  heat  and  salt  by  double  diffusion.  To  study  the  intrusions  and 
possible  driving  mechanisms,  the  salinity  and  temperature  signals  associated  with 
mesoscale  as  well  as  internal  wave  advection  and  straining  must  be  removed  leaving  just 
the  finestructure  thermohaline  anomalies  associated  with  the  intrusions. 

Following  Joyce  (1977)  and  Fedorov  (1986),  three  scales  are  considered  to 
resolve  the  finestructure  thermohaline  anomalies.  Thus  any  scalar  quantity,  0,  is 
described  by 

0  =  0  +  0  +  0'  (3.10) 


where 

0  =  0'=O 
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The  overbar  denotes  the  mean  field  and  is  associated  with  the  large  (mesoscale)  features. 
The  tilde  indicates  those  fine  scale  fluctuations  associated  with  vertical  interleaving 
which  are  anisotropic.  The  accent  denotes  the  micro  scale  fluctuations  associated  with 
dissipation  which  are  turbulent  and  isotropic  in  nature.  Temperature  and  salinity 
anomalies,  T  and  S ,  are  determined  by  removing  the  cruise  averaged  mean  temperatures 
and  salinity,  T  and  S  (Figures  3.9  and_3.10).  To  eliminate  the  effects  of  vertical 
advection  due  to  internal  waves,  T  and  5  are  found  as  functions  of  potential  density 
using  a  cubic  fit  to  the  cruise  average.  The  cruise  averaged  functions  used  are 

f(<T1)  =  a3(26-cr()3  +  fl2(26-a,)2  +  fl1(26-a,)1+a0  n  in 


S (cr, )  =  b, (26  -  a,  )3  +  6,(26  -  a,  )2  +  b,  (26  -  a, )’  +  b 0 


(3.12) 


where 


a0  =  8.5861;  a\  =  4.285;  a2  =  -12.72;  =  5.438; 

b0  =  33.491;  b\  =  2.143;  b2  *  -2.651;  and  b2  =  1.264. 

The  scale  of  the  turbulent  fluctuations  is  not  resolved  by  the  CTD.  Thus,  the  turbulent 
fluctuations  are  assumed  to  be  zero.  The  temperature  and  salinity  anomalies,  f  and  5 , 
are  then  determined  for  each  hydrographic  profile  by 

fn(z)  =  Tn(z)-T((Jln(z)) 


and 


SH(z)  =  Sm(z)-S(ollt(z)) 


where  Tn,  5n,  and  <r,n  and  the  temperature,  salinity  and  density  profiles  for  profile 
number  n.  This  technique  uses  density  as  a  semi-Lagrangian  coordinate  to  remove  the 
effects  cf  internal  wave  vertical  advecdow  md  straining  ieavmg  only  the  temperamre  and 
salinity  signals  associated  with  intrusions  of  different  water  masses. 


Horizontal  displacement,  due  to  internal  waves,  of  water  parcels  from  two 
different  water  masses  is  a  possible  cause  of  vertical  interleaving  near  a  front.  Joyce  et  al. 
(1978)  and  Georgi  (1978)  relate  the  intrusion  temperature  signal,  f,  to  the  horizontal 
displacement  across  the  front  by  internal  waves  and  the  horizontal  temperature  gradient. 
This  relationship  can  be  written  in  vertical  wavenumber  frequency  space  as 


3T  V(co,m) 
dy  ico 


(3.13) 


where  y  is  the  cross  frontal  coordinate  and  V  is  the  cross  frontal  velocity.  An  upper 
bound  on  the  vertical  wavenumber  spectrum  of  T  is  found  by  assuming  are  the  velocity 
variance  is  at  the  inertial  frequency.  The  vertical  wavenumber  spectrum  of  T  is  then 
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where  is  the  spectrum  of  f  and  Ck  is  the  spectrum  of  V  .  Since  the  internal  wave 
field  is  nearly  horizontally  isentropic,  the  velocity  spectrum  can  be  calculated  from  one 
component  of  velocity  estimated  from  the  161  kHz  Doppler  sonar  system  (Appendix  E). 
The  velocity  spectrum  is  determined  using  6.5  minute  averaged  Eulerian  east  velocity, 
from  depths  50  to  306  meters  and  ensemble  averaged  from  yeardays  63  to  78.  For 
comparison,  the  temperature  anomaly  spectrum  (instantaneous  temperature  minus  the 
strain  corrected,  mean  temperature  profile)  from  depths  200-420  meters  is  ensemble 
averaged  from  yeardays  60  to  75.  Since  internal  wave  velocities  should  scale  with 
buoyancy  frequency,  a  velocity  spectrum  calculated  from  shallower  depths  will  over 
estimate  the  velocity  variance  at  200-420  m.  Comparing  the  temperature  anomaly 
spectrum  and  the  velocity  spectrum  reveals  that  internal  wave  horizontal  advection  can 
only  account  for  up  to  5%  of  the  total  observed  temperature  anomaly  signal  (Figure  3.1 1). 
This  calculation  is  sensitive  to  the  horizontal  ’emperature  gradient  estimate.  If  the  true 
gradient  is  twice  as  large  as  the  estimated  gradient,  displacements  from  internal  waves 
could  still  account  for  only  20%  of  the  temperature  variance.  Thus,  the  internal  wave 
field  is  not  strong  enough  to  account  for  the  observed  interleaving  unless  the  horizontal 
gradients  are  much  larger  than  estimated.  This  is  the  same  conclusion  reached  by  Joyce 
et  al.  (1978)  and  Georgi  (1978). 


Since  the  vertical  interleaving  is  not  a  result  of  internal  waves  advection,  it  is 
suggested  that  the  interleaving  is  double  diffusively  driven.  The  potential  for  the  double 
diffusive  processes  to  occur  has  been  shown,  but  no  direct  measurements  of  salt  fingering 
or  diffusive  convection  were  made.  Two  typical  features  of  double  diffusively  active 
intrusions  can  be  compared  to  observations  for  consistency.  These  are  the  expected 
vertical  scale  of  the  intrusion  and  the  slope  of  the  intrusion  across  isopycnals. 

Ruddick  and  Turner  (1979)  used  a  laboratory  model  to  study  mixing  at  an  oceanic 
front  across  which  there  are  large  T-S  anomalies  but  a  small  net  horizontal  density 
difference.  They  found  that  the  vertical  scales  of  the  intrusions  formed  is  directly 
proportional  to  the  horizontal  concentration  gradient  and  inversely  proportional  to  the 
vertical  density  gradient.  They  conclude  that  the  vertical  length  scale  of  cross-frontal 
intrusions,  H,  can  be  determined  to  within  a  factor  of  two  by 


pAS 
1  dp 


p  dz 


(3.15) 


where  n  is  the  density  flux  ratio  (0.56  for  heat/salt  fingers), 
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A5  is  the  salinity  contrast  measured  across  the  front  along  an  isopycnal, 
1  dp 

-  1  do 

P  &  is  the  vertical  density  gradient  and  /?  =  — — 


Ruddick  dnd  Turner  (1979)  successfully  predict  historical  intrusion  scales  near  frontal 
zones.  Using  the  mean  density  gradient  from  200  to  400  meters  (2.36xl0*6  nr1)  and  the 
salinity  change  across  the  front,  the  expected  vertical  length  scale  of  the  intrusions  from 
(3.15)  is  17.8  m. 

Toole  and  Georgi  (1981)  contend  that  it  is  incorrect  to  extrapolate  the  laboratory 
results  of  Ruddick  and  Turner  (1979).  Their  argument  is  that  in  the  laboratory,  the 
interleaving  takes  on  the  same  vertical  scales  as  the  salt  fingers  suggesting  the  intrusion 
scale  is  determined  by  the  energy  in  the  salt  fingering.  In  the  ocean,  the  expected  scales 
of  the  interleaving  is  much  larger  than  the  vertical  scale  of  the  salt  fingers  (which  may  be 
as  large  as  1  m  (Williams,  1975))  and  thus  are  unimportant  to  the  behavior  of  large 
intrusions. 


Toole  and  Georgi  (1981)  use  linear  stability  analysis  of  intrusions  in  a  frontal 
zone  and  look  for  the  fastest  growing  mode.  Their  model  predicts  the  fastest  growing 
mode  as  a  function  of  lateral  salinity  gradient,  the  vertical  density  gradient  of  the  mean 
field,  the  vertical  viscosity,  haline  diffusivity  and  flux  and  the  local  value  of  the  Coriolis 
acceleration.  The  model  uses  the  dimensionless  salinity  gradients,  £x  and  £y,  and  vertical 
wavenumber,  m,  defined  by 


_  (1  ■ -  Y)gPS, 
N2 

_d  -Y)gPSy 

N2 


(3.16) 

(3-17) 

(3.18) 


where  yis  the  ratio  of  the  buoyancy  contributions  of  the  heat  and  salt  fluxes,  As  is  the 
haline  vertical  eddy  diffusivity  and  m  is  the  vertical  wavenumber  with  dimensions  m*1. 
The  fastest  growing  wavenumber  is  given  by 


/(£,) 


m2oi 


(3.19) 


where  /(£,)  is  obtained  from  Fig.  7.  in  Toole  and  Georgi  (1981)  and  cr is  the  ratio  of  the 
haline  to  density  vertical  eddy  diffusivities.  Unlike  the  model  of  Ruddick  and  Turner 
(1979)  which  assumes  a  step  function  in  salinity  across  the  front,  Toole  and  Georgi 
(1981)  use  the  horizontal  salinity  gradient  which  is  not  well  know  for  the  SWAPP 
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experiment.  Only  the  lower  bound  on  the  horizontal  salinity  gradient  is  known.  The 
gradients 


PS,  =7.94X  lO'9  nr1, 

PS,  =  *9.25  X  10-7  m-1  and 
N  =  5.24  X  103  s'1 

and  the  constants  suggested  by  Toole  and  Georgi  (1981) 

/T  =  1.0  X  10-4  m2/s 
a  =0.1 

A1  =  a  A*  =  1.0  X  10*3  m2/s 
y  =0.6 

are  used  to  obtain  to  dimensionless  variables 

£x=  1.13  X  lO*3- 
£z  =  -1.32x  10* 1  and 
/(£,)=0.6. 


Solving  for  the  dominate  vertical  wavelength  of  the 

,5  V* 


H  =  ~  =  2jtoyt 
mr 


f 


exf(£z)N 


intrusions  yields. 
58  meters. 


(3.20) 


This  vertical  scale  is  more  than  twice  as  large  as  those  estimated  using  Ruddick  and 
Turner.  There  are  to  parameters  in  question  when  using  the  Toole  and  Georgi  (1981) 
interleaving  model  which  may  lead  to  this  overestimation.  The  first  is  the  cross  frontal 
salinity  gradient  which  may  be  underestimated.  The  second  is  the  vertical  eddy 
diffusivity,  A*  .  Measurements  taking  by  Gregg  (1989)  suggest  that  A'  may  be  as  small 
as  5  X  10~6  m2/s.  Using 

Av  =  1  x  10-5  m2/s, 

the  Toole  and  Georgi  (1981)  model  yields  a  vertical  intrusion  scale  of 
H  -  18.3  meters. 

This  scale  is  nearly  the  same  as  the  17.8  meter  scale  predicted  by  the  model  of  Ruddick 
and  Turner  (1979)  in  spite  of  the  complaint  by  Toole  and  Georgi  (1981). 

The  observed  intrusion  scales  are  estimated  using  vertical  wavenumber  spectra. 
Half  day  ensemble  averaged  spectra  are  calculated  for  f  and  S  over  the  depth  range  220 
to  420  meters  using  triangle  windowing  on  two  overlapping  1 28  meter  depth  segments 
(Figure  3.13).  These  are  convened  to  gradient  spectra  by  multiplying  by  (2?tm)2,  where 
m  is  the  venical  wavenumber,  and  corrected  for  high  wavenumber  smoothing  by  the 
finite  CTD  resolution.  The  temperature  and  salinity  anomaly  variance  fluctuates 
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throughout  the  cruise  (Figure  3.12)  reaching  a  maximum  on  yearday  68.  The  spectra 
from  the  half  day  averages  corresponding  to  the  times  of  maximum  (yearday  62-62.5)  and 
minimum  (yearday  68.5-69)  variance  are  plotted  along  with  the  15  day  averaged 
spectrum  for  comparison  (Figure  3.14).  The  cruise  average  spectra  are  nearly  white  out 
to  0.1  cpm  and  then  fall  off  slightly  out  to  0.7  c.p.m.  The  minimum  variance  spectrum  is 
red  at  all  wavenumbers  except  for  a  sight  rise  at  0.2  c.p.m.  The  maximum  spectrum  rises 
out  to  0.1  c.p.m.  then  falls  off  at  higher  wavenumber.  The  maximum  spectrum  has 
excess  variance  in  the  wavenumber  band  [0.05-0.2]  c.p.m.  while  the  minimum  spectrum 
lack  variance  in  the  same  band  relative  to  the  mean.  The  expected  wavelength  of  the 
intrusions,  as  predicted  by  Ruddick  and  Turner  (1979)  and  Toole  and  Georgi  (1981),  falls 
within  this  band.  Thus  the  observed  excess  variance  in  the  expected  wavenumber  band 
associated  with  the  intrusions  is  consistent  with  double  diffusive  driving  of  the  intrusions. 

The  slope  of  the  intrusions  across  isopycnals  can  also  be  an  indication  of  active 
double  diffusive  processes.  Turner  (1978)  used  laboratory  experiments  to  show  that  cross 
frontal  intrusions  will  slope  across  isopycnals  surfaces.  His  laboratory  experiments 
demons  jate  that  for  diffusively  unstable  intrusions,  the  increase  in  density  below  a  salt 
fingering  region  dominates  the  decrease  above  the  diffusive  regime.  Thus  warm  salty 
intrusions  will  tend  to  rise  when  surrounded  by  cold,  less  salty  waters.  The  opposite  is 
true  as  well:  cold,  fresh  intrusions  will  tend  to  sink  when  surrounded  by  warm,  salty 
waters.  Similar  conclusions  were  reached  theoretically  by  Stern  (1967).  Ocean 
observations  by  Joyce,  Zenk  and  Toole  (1978)  and  Gregg  (1979)  have  show  intrusions 
crossing  isopycnals  in  the  predicted  manner. 

Tracking  intrusions  along  isopycnals  is  done  here  by  sampling  temperature 
anomalies  on  isopycnal  surfaces.  This  is  the  semi-Lagrangian,  or  isopycnal  following, 
coordinate  system  (Appendix  B).  Here,  isopycnals  that  are  spaced  0.5  m  apart  on  average 
are  tracked.  The  temperature  anomaly  on  an  isopycnal  is  then  just  the  difference  between 
the  average  temperature  on  that  isopycnal  and  the  instantaneous  temperature.  A  cascade 
plot  of  temperature  anomalies  reveals  many  intrusions  which  propagate  across  isopycnals 
(Figure  3.15).  Without  a  cross  section  of  profiles  directly  across  the  front,  it  is  difficult  to 
say  whether  the  intrusions  are  sloping  across  the  from  in  the  correct  sense.  The 
horizontal  advection  of  intrusions  below  FLIP  confuses  the  orientation  of  intrusion 
slopes.  Assuming  the  intrusions  slope  up  from  the  warm  water  encountered  at  the 
beginning  of  the  cruise  to  the  cold  water  encountered  at  the  end  of  the  cruise,  any  mean 
water  velocity  reversal  will  change  the  direction  of  the  slope  in  time.  For  this  reason, 
only  a  subset  of  the  data  where  the  temperature  anomalies  change  most  rapidly  during 
yeardays  [66.7-67.5]  is  investigated  for  intrusion  slope  (Figure  3.16).  By  eye,  there 
appears  to  be  equal  signal  from  upward  and  downward  propagating  intrusions  so  fourier 
analysis  is  used  to  investigate  the  slopes.  Only  the  sense  of  the  slope  can  be  determined 
from  the  time  series.  The  depth  range  [300-464]  m.  is  used  since  temperatures  decrease 
steadily  during  this  time.  A  depth/time  series  consisting  of  128  points  in  depth  and  512 
points  in  time  is  windowed  with  a  triangle  taper  in  both  directions  before  two 
dimensional  Fourier  transformation.  Upward  and  downward  sloping  phase  contributions 
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are  separated  and  summed  over  all  frequency  bands  to  make  vertical  wavenumber 
temperature  anomaly  spectra.  Gradient  spectra  are  found  by  multiplying  the  anomaly 
spectra  by  {Inm)2  (Figure  3.18).  The  ‘upward  sloping’  spectrum  corresponds  to 
intrusions  sloping  upward  from  the  warm,  salty  water  to  the  cold,  fresh  water  (the 
expected  result).  The  upward  sloping  intrusions  contain  two  thirds  of  the  total  variance, 
(f1}-  The  upward  sloping  vertical  wavenumber  spectrum  shows  more  variance  than  that 
downward  sloping  spectrum  at  all  wavenumber  below  0.225  c.p.m.  This  is  the 
wavenumber  band  predicted  to  contain  the  scales  of  the  interleaving  intrusions.  These 
results  are  consistent  with  the  hypothesis  that  the  intrusions  are  actively  double  diffusive. 

3.5.  Finestructure  at  the  Frontal  Interface 

It  is  worth  while  to  take  a  closer  look  at  the  structure  of  the  intrusions  during  the 
yeardays  [66.7-67.5]  when  the  temperature  and  salinity  characteristics  change  most 
rapidly.  There  appear  to  be  two  vertical  scales  associated  with  the  temperature  anomalies 
on  isopycnals  (Figure  3.16).  The  larger  scale  is  40-60  meters  and  takes  the  form  of  noses 
of  warm  water  protruding  into  cold  water.  Two  of  the  warm  noses  ''an  be  seen  at  depths 
220-270m,  day  66  ,  hour  18:00-22:00  and  at  depths  310-360  m,  day  67;  hour  2:00-4:00. 
The  cold  noses  lie  above  and  below  the  warm  noses.  The  small  scale  comprises  the 
structure  on  the  noses  which  has  vertical  scales  of  5  to  15  meters. 

The  potential  for  double  diffusion  to  occur  is  mapped  using  Turner  angle,  Tu, 
calculated  during  yeardays  [66.7-67.5],  Following  the  same  isopycnals  (separated  by  0.5 
m.  in  depth  on  average),  7m  is  calculated  for  each  isopycnal  from  the  temperature  and 
salinity  differences  across  isopycnals  that  are  separated  by  3.0  meters  in  depth  on 
average.  This  yields  a  time/depth  series  of  Tu  (Figure  3.17).  The  characteristic  scale  for 
the  diffusive  and  fingering  interfaces  is  5  to  15  meters.  Alternating  diffusive  and 
fingering  interfaces  are  observed  at  the  edges  of  the  larger  scales  noses  which  appear  at 
nearly  constant  density  or  rise  slightly  with  time.  Diffusive  and  fingering  layers  are 
absent  at  depths  200  to  270  meters  from  5:00-12:00;  day  67  where  the  temperature 
anomalies  suggest  only  one  water  mass  is  present  ( f =-0.3°C).  Away  from  the  noses,  the 
alternating  layers  are  not  as  apparent  and  diffusive  and  fingering  interfaces  tend  to  rise 
and  fall  across  isopycnals  (depths  300-420,  before  to  2:00  hours  ;day67).  The  smaller 
layers  seen  in  the  temperature  anomalies  (5- 15m)  are  double  diffusively  active  above  and 
below  and  their  scales  are  likely  set  by  the  processes  as  described  by  Ruddick  and  Turner 
(1979)  and  Toole  and  Georgi  (1981). 

The  physics  governing  the  larger  noses  (40-60m)  is  not  as  clear.  One  possibility 
is  that  the  noses  are  driven  by  submesoscale  dynamics.  The  dynamics  that  created  the 
front  in  the  first  place  may  have  led  to  this  larger  scale  distortion  on  its  interface.  A 
second  possibility  is  that  strong  shears  for  near  inertial  internal  waves  are  distorting  the 
front.  Schmitt  (198?)  shows  that  salt  fingers  still  form  in  the  presence  of  strong  vertical 
shear  so  the  smaller  layers  are  still  double  diffusively  active  in  the  presence  of  larger 
scale  inertial  shears.  The  shear  analysis  performed  in  Chapter  2  shows  that  the  vertical 


56 


shear  ’  anance  is  dominated  by  near  inertial  waves  with  vertical  wavelengths  of  35-45  m 
which  is  nearly  consistent  with  the  scale  of  the  noses. 

A  third  possibility  is  that  larger  scales  are  linked  in  some  way  to  the  double 
diffusive  processes.  Fedorov  (1986)  performed  an  interesting  experiment.  He  placed  a 
single  (warm  and  salty)  intrusive  wedge  into  a  surrounding  body  of  cold  water.  Salt 
fingers  began  to  form  at  the  bottom  interface  of  the  wedge  and,  in  time  a  counter  flow 
(opposite  to  the  direction  of  the  moving  wedge)  was  observed.  This  led  to  the  formation 
of  a  second  wedge  below  the  original.  Eventually  the  intrusion  split  into  several  laminae. 
The  resulting  structure  was  characterized  by  the  upper  most  layer  protruding  the  farthest 
from  the  source  and  deeper  layers  protruding  successively  less.  This  has  many  of  the 
characteristics  of  the  two  scales  observed  here.  Two  larger  scale  noses  (located  at  [220- 
280]m  and  [3l0-370]m)  are  warm  and  salty,  the  top  of  the  nose  protrudes  the  deepest  into 
the  cold,  fresh  side  of  the  front.  Beneath,  several  smaller  layers  are  found  with  the 
general  shape  of  the  nose  sloping  away  from  the  cold  water.  This  process  still  needs  an 
outside  force  to  define  the  larger  scale  noses  but  may  explain  the  generation  of  the 
smaller  scale  layers  and  their  link  to  the  larger  scales. 

Newly  formed  intrusive  lenses,  referred  to  as  ‘calving’  by  Fedorov  (1986),  are 
another  feature  seen  during  yeardays  [66.7-67.5].  Calving  is  the  separation  of  an 
intrusive  tongue  from  the  frontal  interface  and  its  transformation  into  an  isolated  lens. 
Several  lenses  can  be  seen  in  Figure  3.16.  One  such  feature  occurs  at  300  m.  at  20:00; 
day  66;  where  a  lens  of  warm  water  is  surrounded  by  colder  water.  Without  a  three 
dimensional  survey  i.  is  impossible  to  tell  whether  this  lens  is  completely  separated  form 
it  origin.  Separated  lenses  near  frontal  zones  have  been  observed  before  (Fedorov,  1986). 
These  lenses  are  likely  formed  from  motions  in  the  horizontal  plane.  The  calving  of 
interleaving  layers  will  lead  to  increased  variability  in  isopycnal  temperatures  near  a  front 
and  lead  to  enhanced  horizontal  and  vertical  mixing. 

3.6.  Estimates  of  Horizontal  Heat  and  Salt  Fluxes 

In  order  to  maintain  a  statistically  steady  state  of  interleaving  temperature  and 
salinity  fields  in  the  presence  of  small-scale  dissipation,  fluxes  of  heat  and  salt  must  be 
supplied.  The  interleaving  acts  as  stirring  which  enlarges  the  surface  area  between  the 
two  water  masses.  Inversions  in  temperature  and  salinity  due  to  the  interleaving  permit 
small  scale  enhanced  vertical  transfer  of  heat  and  salt  by  double  diffusive  processes. 
Horizontal  temperature  and  salinity  fluxes,  Fj  and  F$,  across  the  from  are  estimated  by 

F7=-{vf)=V(?,’)f,"'  (3-21) 

Fs=-(V5}=4/(S,,)JJ-'.  (3.22) 

where  V  is  the  cross  frontal  velocity  and  Aj  and  Af  are  vertical  eddy  diffusivities  of 
temperature  and  salt  (Joyce,  1977).  ($,2\  and  (^V)  are  the  variances  of  vertical 
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gradients  of  salinity  and  temperature  associated  with  vertical  interleaving.  (There  is  a 
slight  confusion  as  to  which  scales  to  use.  Gargett  (1978)  asserts  that  Joyce’s  model  uses 
the  microscale  fluctuations  of  T  and  S.  However,  Joyce  (1977)  and  Fedorov  (1986) 
clearly  state  that  the  balance  is  with  the  interleaving  scales.  Joyce  et.  al  (1978)  use  the 
interleaving  scales  to  estimate  horizontal  fluxes  with  Joyce’s  (1977)  model.)  At  some 
small  vertical  scale  it  is  expected  that  dissipation  scales  will  become  important  and  using 
a  vertical  eddy  diffusivity  is  no  longer  valid.  Measurements  by  Gregg  (1975)  suggest  that 
dissipation  dominates  at  scales  smaller  that  0.5  c.p.m.  The  CTD  deployed  during 
SWAPP  has  a  vertical  resolution  of  0.7  c.p.m  (Appendix  C).  (s,2^  and  y2]  over  the 
wavenumber  band  [0.0078-  0.5]  c.p.m.  are  estimated  from  half  day  ensemble  averaged 
spectra  of  f  and  S  (Figure  3.12).  Using  the  maximum  values 

(■^L'2-72*10-’  (VJ 
{f-2L=8-81*10"  C %S 

and  Aj  -  A /  =  .1  cm2/s  (the  value  suggested  by  Joyce  et  al.  (1978)),  the  corresponding 
horizontal  temperature  and  salt  fluxes  are 

Fr  =  1.67x10^  (°Cm/s)  (3.23) 

and  Fs  =  2.62 x  10'5  (%ct/{).  (3-24) 

These  estimates  are  considered  as  upper  bounds  on  the  fluxes  since  the  cross  frontal 
gradients  are  likely  underestimated  and  (s2\  and  (j2^  are  taken  at  their  maximums.  The 
temperature  and  salt  fluxes  are  20%  ana  4%  respectively  of  those  found  at  the  Antarctic 
polar  front.  The  uncertainty  of  the  size  of  the  vertical  eddy  diffusivities  raises  some 
question  as  to  the  proper  bound  of  the  horizontal  fluxes. 

Assuming  all  the  dissipation  is  from  shear  instabilities,  the  vertical  eddy 
diffusivity  of  density,  Kp,  can  be  estimated  using  Gregg’s  (1989)  ten  meter  shear  variance 
scaling  of  dissipation.  K p  is  estimated  by 

^<^^  =  1.4xlO-10N0-2^r  (m2  /  s)  (3.25) 

where  eiw  is  the  estimated  dissipation,  N  is  the  buoyancy  frequency,  Sjo  is  the  measured 
ten  meter  shear,  Sgm  is  the  GM76  shear  and  No  is  the  normalized  buoyancy  frequency 
(5.2X10'3  1/sec).  Using  the  shear  variance  estimates  from  the  161  kHz  Doppler  sonar 
system  (Table  2.1),  yields  Kp  <  8.7x1  O'6  mtys.  This  value  is  nearly  equal  to  the  eddy 
diffusivities  suggested  by  Joyce  et  al.  (1978).  The  ten  meter  shear  scaling  does  not 
include  double  diffusively  driven  mixing.  For  individual  intrusions,  Aj  and  A/  may  be 
as  large  as  l.OxlO-3  rr^/s  (Fedorov;  1986;  Horne  1978)  and  assuming  that  50%  of  the 
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water  column  is  subject  to  double  diffusive  fluxes  (Figure  3.8)  at  the  peak  of  interleaving, 
the  vertically  averaged  vertical  diffusivities  might  be  as  large  as  5.0X10"4  m-/s.  (It  is 
important  to  note  that  if  the  double  diffusive  fluxes  are  dominating  then  it  is  likely  that 
\s  >  Aj  as  reflected  by  the  intrusion  slopes  across  isopycnals.)  Thus  it  is  possible  that 
the  upper  bounds  on  the  horizontal  fluxes  may  be  even  larger  but  since  the  cross  frontal 
temperature  and  salinity  gradients  a/e  underestimated,  (3.23)  and  (3.24)  are  taken  as 
conservative  upper  bounds. 

The  calculation  of  f  and  S  presented  here  effectively  removes  the  mean 
temperature  and  salinity  on  isopycnals.  If  the  water  was  horizontally  homogeneous,  T 
and  5  would  always  be  zero.  Joyce  et  al.  (1978)  removes  the  mean  temperature  and 
salinity  at  depth  and  doing  so  still  leaves  a  signal  from  internal  wave  straining  and 
advection  of  the  mean  profile.  Even  in  a  horizontally  homogeneous  ocean,  Joyce  et  al. 
(1978)  would  have  non  zero  f  and  S  and  would  predict  a  ‘cross  front’  flux. 

The  question  of  buoyancy  fluxes  across  a  frontal  interface  has  been  a  subject  of 
some  confusion.  Joyce  et  al.  (1978)  used  the  above  model  (Equations  (3.21)  and  (3.22)) 
to  estimate  salt  and  heat  fluxes  across  the  polar  front  at  Drakes  passage.  In  turn,  they 
estimated  the  buoyancy  flux  by 


Ft  =  -(vf)  =  P'(pF,-aFT). 


(3.26) 


and  found  that  the  buoyancy  flux  was  dominated  by  the  salt  flux  working  against  the 
mean  horizontal  density  gradient.  Toole  and  Georgi  (1981)  argue  that  these  results  are 
just  an  artifact  of  overestimating  the  salinity  variance  due  to  uncertainties  in  the  CTD 
sensor  responses.  The  dynamical  model  of  Toole  and  Georgi  (1981)  can  not  yield  steady 
state  fluxes  across  a  front  but  their  results  suggest  that  the  density  flux  must  be  down 
gradient  in  the  Antarctic  Polar  Front.  Fedorov  (1986)  disputes  the  results  of  Joyce  et  al. 
(1978)  suggesting  that  the  whole  concept  of  eddy  diffusivities  is  inappropriate.  Using  the 
upper  bounds  on  the  horizontal  fluxes,  (3.26)  yields 


F P  =Po(2-0  x  10-s  -  2.5  x  10'8)  kg/ 


m  sec 


=  -5  x  10-6  kg/  2 

/  nr  sec 


(3.27) 


Here,  the  buoyancy  flux  is  dominated  by  the  heat  flux  contribution.  The  water  in  the  300 
to  400  meter  depth  range  got  steadily  more  dense  as  well  as  colder  with  time.  Thus  the 
net  buoyancy  flux  is  up  gradient  rather  than  down  gradient  yielding  the  same  results 
reached  by  Joyce  et  al.  (1978)  at  the  Antarctic  Polar  Front.  This  result  does  not  seem 
physically  realistic  since  the  horizontal  density  gradient  appears  to  be  supported  by  a 
geostrophic  balance.  The  buoyancy  flux  calculation  is  presented  here  just  as  an 
illustration.  It  is  not  taken  as  a  physically  meaningful  result. 


This  raises  some  question  about  the  appropriateness  of  using  (3.26).  The 
assumption  that  Aj  =  Af  is  likely  incorrect  especially  if  there  is  significant  double 
diffusively  driven  vertical  fluxes.  For  diffusive  interfaces  on  intrusions,  heat  fluxes  are 
larger  than  salt  fluxes  ( Aj  >  Avs),  and  for  salt  fingering  interfaces,  salt  fluxes  are  larger 
than  heat  fluxes  ( Aj  <  Aj).  Since  salt  fingering  interfaces  occur  at  least  twice  as  often 
as  diffusive  ones  (Figure  3.8),  it  is  likely  that  the  average  vertical  haline  eddy  diffusivity 
is  larger  (Aj  <  A/).  It  would  not  take  much  more  haline  vertical  diffusivity  to  change 
the  up  gradient  buoyancy  flux  to  a  down  gradient  flux.  Consider  a  case  where  semi- 
Lagrangian  (isopycnal)  coordinates  are  used  to  track  temperature  and  salinity  anomalies. 
The  ratio  of  /3Fs  to  aFj  is  equal  to  one  since  aAT  =  /3AS  on  constant  density  surfaces 
and  the  horizontal  buoyancy  flux  must  be  zero.  Assuming  that  Aj  *  Af ,  would  lead  to  a 
net  buoyancy  transport  along  isopycnals  which  is  physically  impossible.  Only  in  the 
case  of  using  depth  coordinates  with  sloping  isopycnals  can  a  net  buoyancy  flux  exist. 
This  is  because  the  depth  coordinates  cross  isopycnals  themselves  and  the  model 
correctly  reflects  this  ‘flux’  of  isopycnals. 

3.7,  Summary 

Observations  of  increased  variance  in  thermohaline  finestructure  during  the 
SWAPP  experiment  is  linked  physical  processes  associated  with  a  frontal  zone.  The 
front  appears  to  be  in  geostrophic  balance.  The  front  is  characterized  in  the  depth  range 
[200-420]  m.  by  warm  water  to  the  southeast  and  warm  water  to  the  northwest. 
Horizontal  salinity  and  temperature  gradients  across  the  front  are  at  least  1.0  xlO'5  %d m 
and  5.3x10-5  °c/m. 

The  fine  scale  thermohaline  structure  is  sampled  as  a  CTD  time/depth  series  taken 
from  a  moored  FLIP  while  the  frontal  zone  is  advected  past  the  sampling  location.  The 
time  series  is  characterized  by  warm,  salty  water  starting  at  yearday  60,  1990.  A  rapid 
drop  in  isopycnal  temperatures  occurs  near  yearday  67  as  the  water  characteristics 
become  colder  and  fresher.  Individual  temperature  and  salinity  profiles  are  characterized 
by  many  salinity  and  temperature  inversions.  Estimates  of  Turner  angles  suggest  the 
presence  of  double  diffusive  processes.  The  potential  for  salt  fingering  (Tu  >  45°)  is 
larger  than  the  potential  for  diffusive  instabilities  (Tu  <  -45°)  and  the  number  of 
temperature  and  salinity  inversions  rises  yearday  67  and  thus  precedes  the  front. 

Temperature  and  salinity  anomalies  are  calculated  by  removing  the  mean 
temperature  and  salinity  on  isopycnals.  This  technique  allows  the  tacking  of  intrusions. 
The  intrusions  have  two  dominant  vertical  scales  near  the  front.  The  small  scale 
intrusions  have  vertical  lengths  of  5-15  meters,  are  driven  by  double  diffusive  processes, 
slope  upward  from  the  warm  side  to  the  cold  side  of  the  front  and  is  consistent  with  scales 
predicted  by  Ruddick  and  Turner  (1979)  and  Toole  and  Georgi  (1981).  The  larger  scale 
is  40-60  m  and  is  likely  determined  by  submesoscale  or  internal  wave  forcing. 
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Figure  3.1  T-  S  Diagram  from  200-420m.  Scatter  plot  of  salinity  plotted  versus  potential  temperature. 
A  point  is  plotted  for  every  3  meters  from  200-420m.  One  profile  from  the  lower  CTD  is  plotted  every  6 
hours.  Two  distinct  water  masses  can  be  seen.  The  cruise  averaged  profile  (solid  line)  and  a  cubic  curve 
fit  (dashed  line)  are  also  plotted. 
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Constant  Potential  Density  Hour  averaged  tcmpcraturc(a)  and  saiiniiy(b)  arc  found  along  three  isopycnals  (z  =  250,  300 
versus  time. 


Northward  Displacement  (Km) 


Figure  3.4  Integrated  Current  Trajectory  Hour  averaged  velocity  profiles  are  depth  averaged  from  200 
to  350  m  depths  and  integrated  in  time. 
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jre  3.5  Geopotential  Anomalies.  Hour  averaged  (light  line)  and  daily  averaged  (dark  line)  gcopolcntial  anomalies  arc  plotted  verses  lime, 
anomalies  arc  calculated  over  4  separate  depth  ranges:  (0-100|m,  I100-2001m,  (200-300Jm  and  [300-400)m. 
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Figure  3.6  Typical  Temperature  and  Salinity  Profile  from  Yearday  67.  A  typical  temperature  and 
corresponding  salinity  profile  are  plotted  versus  pressure.  Temperature  inversions  are  matched  by 
salinity  inversions  to  yieid  a  stable  density  profile.  Typical  scales  of  temperature  inversions  are  5  to  15 
meters. 
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Figure  3.7  Typical  Turner  Angle  Profile  from  Yearday  67.  Same  Figure  3.6  except  Turner  angle,  Tu. 
plotted.  Tu  is  computed  from  temperature  and  salinity  differences  taken  over  5  meters  in  depth. 
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Figure  3.8  Percent  of  Profile  that  is  Subject  to  Double  Diffusive  Convection.  Turner  angle,  Tu,  is 
computed  from  temperature  and  salinity  differences  taken  over  5  meters  in  depth,  The  percent  of  each 
profile  that  is  subject  to  diffusive  instabilities  (Tu  <  -45°;  dashed  linej  and  salt  fingering  (Tu  >  45°;  solid 
line)  is  averaged  over  half  day  intervals  and  plotted  versus  time.  During  strong  interleaving  on  yeardays 
66  to  70,  as  much  as  50%  of  the  water  column  is  subject  to  double  diffusive  fluxes. 
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Figure  3.10  <3%  -  T  Diagram  from  200-420m.  Scatter  plot  of  temperature  plotted  versus  potential 
density.  A  point  is  plotted  for  every  3  meters  from  200-4 20m.  One  profile  from  the  lower  CTD  is  plotted 
every  6  hours.  Two  distinct  water  masses  can  be  seen.  The  cruise  averaged  profile  (solid  line)  and  a 
cubic  curve  fit  (dashed  line)  are  also  plotted. 


Figure  3.11  Internal  Wave  Temperature  Anomalies.  Temperature  anomaly  spectrum,  ®f,  is  calculated 
using  temperature  fluctuations  from  depth  ranges  220-420  meters  and  ensemble  averaged  over  yeardCj  * 
[60-75].  Maximum  internal  wave  temperature  anomaly  specr  m  is  calculated  using  horizontal 
displacement  spectrum  (from  6.5  minute  averaged  velocity  profiles,  over  depth  ranges  50-30f  n  and 
ensemble  averaged  over  yeardays  [63-77])  and  scaling  by  the  cross  frontal  temperature  gradient  squared. 
Internal  waves  can  only  account  for  up  to  5%  of  observed  vertical  interleaving. 
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Figure  3.13  Vertical  Wavenumber  Spectra  of  Temperature  and  Salinity  Anomalies.  Spectra  are 
Iculated  over  the  depth  ranges  220-420  meters.  Three  ensemble  averaged  spectra  arc  shown.  Two  half 
->ay  averages  from  yeardays  [62-62.5]  and  [68.5-69]  and  a  15  day  average.  Confidence  intervals  are 
approximate  for  half  day  ensembles. 
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Figure  3.14  Vertical  Wavenumber  Spectra  of  Gradient  Temperature  and  Salinity  Anomalies.  Spectra  are 
calculated  over  the  depth  ranges  22CM20  meters.  Three  ensemble  averaged  spectra  arc  shown:  yeardays 
{62-62.5],  [68.5-69]  and  a  15  day  average.  Half  day  ensembles  are  smoothed  in  frequency  with  5  pL 
triangle  window. 


Isopycnal  Mean  Depth  (meters) 


Temperature  Difference  (°C) 


Figure  3.15  Temperature  Anomaly  Cascade  Plot  .  Temperature  anomaly  on  an  isopycnal  is  plotted 
versus  mean  isopycnal  depth.  Cascade  represents  1 .8  hours  in  lime. 
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Figure  3.18  Sloping  Temperature  Gradient  Anomaly  Wavenumber  Spectrum.  Temperature  anomalies  on 
isopycnals  are  tracked.  Upward  and  downward  phase  wavenumber  spectra  are  extracted  from  vertical 
wavenumber/  frequency  spectra.  Wavenumber  temperature  gradient  anomaly  spectrum  corresponding  to 
upward  (heavy  line)  and  downward  (light  line)  propagating  phases  arc  plotted.  Data  is  from  yeardays 
[66.7-67.5]  and  depths  [300-364]  meters.  Upward  sloping  intrusions  correspond  to  warm,  salty  (cold, 
fresh)  water  intrusions  rising  (sinking)  across  isopycnals  and  they  move  into  cold,  fresh  (warm,  salty) 
water. 
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Chapter  4 


Froude  Number  Distribution 


4.1.  Introduction 


Turbulence  observed  in  the  ocean  interior  is  thought  to  be  caused  by  breaking 
internal  waves  .  Shear  instability  is  the  hypothesized  cause  of  breaking  (Garrett  and 
Munk,  1972;  Eriksen,  1978;  Munk,  1981,  Kunze,  1990).  A  linear  superposition  of 
internal  waves  will  lead  to  local  vertical  shearing  and  vertical  straining  of  the  water 
column.  Shear  and  strain  levels  are  occasionally  sufficient  to  cause  a  shear  instability. 
The  distribution  of  these  shear  instabilities  is  not  well  known.  Slow  profiling 
microstructure  instruments  may  undersample  the  underlying  distribution  of  turbulent 
events  due  to  the  intermittency  of  wave  breaking  events  (Baker,  1985).  This 
undersampling  will  likely  lead  to  an  underestimate  of  vertical  mixing  due  to  internal 
wave  breaking.  This  chapter  examines  the  underlying  probability  density  functions  of 
shear  and  strain  and  presents  a  simple  statistical  model  for  the  distribution  of  shear 
instabilities,  or  more  specifically,  Froude  and  Richardson  numbers. 


The  gradient  Richardson  number,  Ri,  is  often  used  to  describe  ratio  of  buoyancy 
forces  to  inertial  forces.  It  is  typically  defined  by 


gP,  _ 
Pc'7,2  ~  S2  ‘ 


N2 


(4.1) 


where  g  is  the  local  acceleration  of  gravity,  po  is  the  mean  density,  pz  is  the  vertical 
potential  density  gradient,  Vz  or  5  is  the  absolute  value  of  the  vertical  shear  and  N  is  the 
buoyancy  frequency.  For  Ri  <  0.25,  a  fluid  has  the  potential  for  a  shear  instability  leading 
to  turbulence.  Since  shear  variance  is  typically  larger  than  N  variance,  it  is  useful  to  use 
to  buoyancy  normalized  shear  or  Froude  Number,  Fr,  defined  by 

-i  y 

Fr  =  Ri  2  =  —  (4.2) 

N 


to  reduce  the  variance  in  the  dimensionless  ratio  of  inertial  to  buoyancy  forces. 

The  Froude  number  distribution  model  is  formed  here  using  the  known 
distributions  of  shear  and  strain.  This  is  very  similar  to  the  derivation  of  a  pdf  of  Ri  by 
Desaubies  and  Smith  (1982).  The  main  difference  is  the  underlying  distribution  of 
buoyancy  frequency  used.  Desaubies  and  Smith  (1982)  assume  that  N  is  normally 
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distributed.  However,  this  assumption  is  not  consistent  with  resent  observations.  Pinkel 
and  Anderson  (1992)  show  that  fine  scale  strain  statistics  can  be  modeled  using  Poisson 
statistics.  From  this  model,  the  probability  density  of  N  is  estimated  and  compared  to 
histograms  derived  from  the  SWAPP  CTD  observations.  Thus,  the  pdf  of  N  is  not  a 
normal  distribution.  The  shear  distribution  used  here  is  the  same  as  that  assumed  by 
Desaubies  and  Smith  (1982).  Individual  horizontal  components  of  vertical  shear  as  seen 
by  the  sonar  appear  normally  distributed  and  independent.  This  suggests  a  Rayleigh 
distribution  of  5.  is  the  appropriate  model  Joint  histograms  of  5  and  N  show  that  they  are 
nearly  independent  which  allows  the  direct  estimation  of  the  Froude  number  distribution. 
The  modeled  Froude  number  distribution  is  compared  to  histograms  of  Fr  derived  from 
the  SWAPP  data. 


The  density  and  velocity  data  presented  in  this  chapter  are  collected  with  the  MPL 
rapidly  profiling  CTD  and  161  kHz  Doppler  sonar  systems  deployed  from  R/P  FLIP 
during  the  SWAPP  (1990)  experiment  (Appendix  A).  The  CTD  system  provides  density 
profiles  from  the  surface  to  420  m.  every  130  seconds  with  a  vertical  resolution  of  1.5  m 
(Appendix  C).  The  sonar  system  yields  1  minute  averaged  velocity  profiles  from  50  to 
325  m.  with  a  vertical  pulse  length  of  5.5  m  (Appendix  E).  The  data  sets  are  depth 
calibrated  as  described  in  Section  2.2.  After  depth  calibration,  the  velocity  profiles  are 
smoothed  in  time  with  a  3  minute  running  mean  filter  and  resampled  at  the  30  seconds 
after  the  start  of  each  CTD  profile.  Then  both  data  sets  are  combined  yielding  more  that 
9200  profiles  of  simultaneous  velocity  and  density  estimates  from  50  to  325  m.  every  130 
seconds  from  yearday  63  to  75.  In  this  study,  nine  25  m.  depth  ranges  are  investigated 
from  75  to  300  m.  Vertical  shear  and  buoyancy  frequency  are  estimated  from  this 
combined  data  set.  Shear,  5,  is  defined  as 
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where  z  is  the  vertical  coordinate  (zero  at  the  surface  and  increasing  with  depth),  Az  is 
the  vertical  differencing  interval,  u  is  the  east  component  of  velocity  and  v  is  the  north 
component  of  velocity.  The  buoyancy  frequency,  N,  is  defined  as 


I  r  .  ‘ 


2 

g 

( 

f  All 

_  rt 

fz 

y 

.Po  dz. 

_P0Az 

Z  f  — 

l  2  J 

P 

l*  2  , 

)  _ 

(4.4) 


where  p  is  the  density  profile,  g  is  the  acceleration  of  gravity  and  po  is  the  average 
density.  Fr  is  calculated  from  the  combined  velocity  and  density  profiles  at  depth  z  by 


(4.5) 
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where  S  and  N  are  taken  as  defined  in  (4.3)  and  (4.4). 

4.2.  Distribution  of  Buoyancy  Frequency  and  Vertical  Shear 

Joint  histograms  of  shear,  S,  and  buoyancy  frequency,  N ,  are  constructed  for  nine 
25  m  depth  ranges  from  75  to  300  m.  5  and  N  are  calculated  using  (4.3)  and  (4.4)  and  a 
difference  interval,  A z,  of  10  m.  Each  histograms  contains  5  estimates  of  S  and  N  in 
depth,  and  9200  estimates  in  time.  The  joint  histograms,  are  normalized  such 

that 

£'X/(W,.S,) ANAS  =  1 

'=0  >= 0  (4.6) 

where  N.  =  (i  +  ^)AN,  Sj  =  (;'  +  -^)A5 

and  i  and  j  are  the  bin  numbers.  The  joint  histogram,  are  contour  plotted  for 

each  of  the  nine  depth  ranges  (Figure  4.1).  It  appears  that  the  total  variance  in  the 
histograms  decreases  with  depth.  The  mean  shear  and  buoyancy  frequency  decreases 
with  depth  below  125  m.  A  tilting  of  the  contours  from  low  S  and  N  to  high  S  and  N  can 
bee  seen  in  the  histograms  above  250  m.  This  suggests  some  dependence  of  shear  on 
buoyancy  frequency. 

Histograms  of  buoyancy  frequency,  fN(N),  are  found  by  integrating  the  joint 
histogram  over  all  the  shear  bins, 

/«(^)=IM,S;)AS,  (4.7) 

and  plotted  for  each  of  the  nine  depth  ranges  (Figure  4.4).  The  mean  and  variance  of  the 
distributions  are  given  in  Table  4.1.  The  mean  has  a  maximum  at  [100-125]  m.  depth 
range  and  the  variance  decreases  steadily  with  depth  below  75  m. 

Histograms  of  vertical  shear,  fs{S),  are  found  by  integrating  the  joint  histogram 
over  all  the  buoyancy  frequency  bins.  fs(S )  is  calculated  by 

fsiS^^fW^AN  (4.8) 

i=0 

and  plotted  for  each  of  the  depth  ranges  (Figure  4.5).  The  mean  and  variance  (Table  4.1), 
like  the  buoyancy  frequency,  reaches  a  maximum  at  [100-125]  m.  and  decreases  steadily 
with  depth.  An  increase  in  shear  is  seen  at  [275-300]  m.  ranges  where  velocity  estimates 
have  reduced  signal  to  noise  at  the  far  range  of  the  sonar. 
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The  dependence  of  shear  and  buoyancy  can  be  tested  using  the  joint  histograms. 
Some  dependence  is  anticipated  since  both  the  mean  shear  and  buoyancy  frequency 
decrease  with  depth.  Vertical  advection  through  the  fixed  depth  range  bins  will  cause 
some  dependence  due  to  the  vertical  gradients  of  S  and  N .  Since  both  decrease  with 
depth,  a  positive  correlation  is  expected.  S  and  N  are  said  to  be  independent  if 

f(N,S)  =  f\N)fs(S)  (4.9) 

Contour  plots  of  fs  (N)f 5 (S)  for  the  nine  depth  ranges  looks  qualitatively  nearly  the 
same  as  those  of  f(N,S)  except  for  the  slight  tilt  in  the  contours  (Figure  4.2).  The 
differences  are  better  illustrated  by  subtracting  the  independent  distribution  from  the  true 
distribution.  By  contouring  only  the  regions  where  f(N,S)  >  fN(N)fs(S)  ,  the  positive 
correlation  between  5  and  N  is  clearly  illustrated  (Figure  4.3).  Typically  a  Chi-squared 
test  can  be  applied  to  determine  if  two  variables  are  independent  but  these  tests  are 
extremely  sensitive  to  oversampling  of  a  distribution  which  is  the  case  here  (Pinkel  and 
Anderson,  1992).  The  correlation  coefficients  are  low  (<0.15)  but  positive  for  all  the 
depth  ranges.  Clearly  there  is  a  slight  dependence  of  S  and  N,  however  this  dependence 
appears  to  be  very  small.  Kunze  (1990)  find  that  shear  and  strain  are  uncorrelated  and 
suggest  that  most  of  the  fine  scale  shear  comes  from  near  inertial  waves. 

The  distribution  of  N  can  be  modeled  from  statistics  of  fine  scale  strain.  From 
Pinkel  and  Anderson  (1992),  strain,  y,  sampled  in  fixed  depth  coordinates,  is  modeled  as 
a  Poisson  process.  The  resulting  distribution  of  yis 

pE(y\Xz)=a?^J*'-  (4.10) 


where 


A  =  1.1m'1,  A*  =  AzA,  y  =  ^ZA-  -  t\  +1 

/A  z  1 


and  Tjt  is  the  Eulerian  strain  (Appendix  B).  To  estimate  the  corresponding  buoyancy 
frequency  distribution,  the  strain  is  defined  by 


y  =  t?2  +  1  =  ~r 
P, 


(4.11) 


where  N  is  the  mean  buoyancy  frequency  and  N  is  calculated  from  (4.4).  This 
transformation  assumes  the  N  does  not  vary  with  depth  This  yields  the  distribution 
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Note  that  Az  is  the  diffe_rencing  interval  used  in  (4.4).  Taking  N  from  each  of  the  depth 
ranges  (Table  4.1)  and  Az=10.0  m.,  is  fit  to  the  histograms  of  N  (Figure  4.4).  Good 
qualitative  agreement  with  sampled  distributions  at  almost  all  depths  can  be  seen  by 
plotting  the  modeled  pdf’s  with  the  histograms  of  N.  The  [75-100]  m.  depth  range  does 
not  yield  good  agreement  to  the  Poisson  model.  Apparently,  mixed  layer  dynamics  are 
significantly  different  than  those  in  the  thermocline.  The  [125-150]  m.  depth  range  does 
not  appear  to  fit  well  with  the  model.  This  is  likely  due  to  rapid  changes  in  N  at  these 
depths. 

To  model  the  distribution  of  S  it  is  useful  to  examine  the  individual  components 
of  shear,  Sx  and  Sy,  where 

S=luJ=(S,2+S;p  (4.13) 


The  two  components  of  shear  sampled  by  the  sonar  have  normal  distributions,  have  small 
correlation  coefficients  and  are  thus  independent.  This  independence  and  normal 
distribution  of  the  shear  components  has  been  seen  by  Gregg  (1991)  and  was  assumed  in 
the  Garret  and  Munk  internal  wave  models.  The  square  root  of  the  sum  of  two 
independent  and  normally  distributed  variables  is  described  by  the  Chi  distribution  with  2 
degrees  of  freedom  (Papoulis,  p.105  &  p.250).  The  Chi  distribution  for  the  case  2 
degrees  of  freedom  (also  called  the  Rayleigh  distribution)  is  defined  for  S  as 


P5(S)  = 


(4.14) 


where 


S  = 


and 


a 

2 


The  shear  distribution,  ps,  is  calculated  for  each  of  the  nine  depth  ranges  using  the  mean 
from  the  corresponding  histograms  (Figure  4.5).  Good  qualitative  agreement  to  the 
sampled  shear  histograms  can  be  seen  at  all  but  the  [75-100]  m.  depth  range  where  mixed 
layer  dynamics  may  be  playing  an  important  role  in  shear  statistics. 

4.3.  Froude  Number  Distribution 

Histograms  of  Fr  are  calculated  over  the  same  nine  depth  ranges  as  used  for  the 
shear  and  buoyancy  frequency  distributions  (Figure  4.6).  Between  100  and  250  meters, 
the  distributions  of  Fr  appear  nearly  identical  with  little  change  in  mean  and  variance 
(Table  4.1).  Increased  mean  and  variance  is  observed  below  250  m.  and  above  100  m. 

The  distribution  of  Fr  is  determined  using  the  joint  distribution  of  S  and  N.  For 
two  variables,  x  andy,  for  which  the  joint  distribution,  f(x,y),  is  known,  the  distribution  of 
the  ratio  of  the  two  variables,  v,  is  given  by  (Papoulis,  p.  197) 
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(4.15) 


A(v)=  Jbl/C'O'.y)^ 


where  v  =  — 

y 


Since  the  Froude  number  is  the  ratio  of  S  and  N,  its  distribution  is  given  by 

pFr(Fr)  =  J  Np(FrN,N)dN 
o 


(4.16) 


where  p(S,N)  is  the  joint  distribution  of  S  and  N.  Assuming  that  S  and  N  are 
independent,  the  joint  distribution  is 


p(S,N)  =  ps(S)pN(N ) 
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(4.17) 


The  joint  histograms  demonstrated  that  these  variables  are  not  completely  independent, 
but  the  dependence,  especially  below  125  m.,  appears  to  be  small  and  correlation  factors 
estimated  for  S  and  N  at  constant  depth  are  less  than  0.15.  Substituting  (4.17)  into  (4.16) 
yields 
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(4.18) 


By  solving  the  integral  in  (4.18)  and  reducing,  we  obtain 


pFr(Fr) 


2qFr{X  + 1) 
(qFr2  +  \f*2 


(4.19) 


where 

N2  _  n  N2  _  \  N2 
q  2a2/.*  4 A*  S2  4 A*  cr] 


(4.20) 


The  probability  density  distribution  function  pFr ,  is  controlled  by  the  parameters  A*  and 
q.  The  A*  parameter  (Pinkel  and  Anderson;  1991)  determines  the  general  shape  of  the 
distribution  and  contains  an  strain  scale  factor  A  =  1.1  m'1,  and  the  vertical  differencing 
interval,  Az.  The  q  parameter  scales  the  distribution  and  is  a  function  of  A*,  the  mean 
buoyancy  frequency,  N ,  and  the  shear  '-ariance,  cr2.  The  probability  of  the  exceeding 
the  shear  instability  threshold  is  found  by  integrating  (4.19) 


84 


P(Fr  >  2)  =  J  pF'(Fr)dFr 

=  (4q  +  \){l+l)  (4.21) 

From  (4.21)  it  can  be  seen  that  increasing  the  shear  variance  or  decreasing  the 
differencing  interval  will  increase  the  likelihood  of  a  shear  instability. 

The  modeled  Froude  number  distributions  are  fit  to  the  histograms  in  the  nine 
depth  ranges  using  the  mean  Froude  numbers  estimated  from  the  histograms  (Table  4. 1 
and  Figure  4.6).  Good  model  agreement  can  be  seen  from  125  to  275  meters,  and  the 
Froude  number  and  shape  of  the  distributions  are  nearly  constant  with  depth.  At 
shallower  depths,  [75-125]  m.,  the  model  fails  since  shear  and  strain  depart  from  their 
modeled  distributions  and  have  increased  dependence. 

The  model  should  b*  'alid  for  different  scales  assuming  that  shear  and  strain 
estimates  are  independent,  the  model  distribution  from  Pinkel  and  Anderson  (1992) 
includes  scales  down  to  2  m.  The  individual  components  of  shear  are  still  likely  to  be 
independent  and  normally  distributed  over  a  broad  range  of  scales.  To  demonstrate  this 
scaling,  the  model  is  compared  to  data  collected  by  Kunze  (1990)  using  the  neutrally 
buoyant  float,  RiNo,  during  PATCHEX  (1985)  of  the  coast  of  central  California.  The 
float,  deployed  at  depths  180-200  m  from  October  12  to  22,  yielded  simultaneous  shear 
and  strain  estimates  calculated  over  1  m,  2.5  m,  and  5  m.  separations.  Using  the  shear 
and  strain,  Kunze  (1990)  calculate  Froude  Numbers  and  present  histograms 
corresponding  to  1  m,  2.5  m,  and  5  m  separations.  Kunze  (1990)  shows  that  the 
distribution  of  Froude  numbers  changes  form  during  the  presence  of  an  inertial  wave 
group.  The  model  distribution  is  fit  to  those  of  Kunze  (1990)  by  using  Figure  13  (from 
Kunze)  to  estimate  the  mean  shear  to  mean  buoyancy  frequency  ratio  for  each  separation 
and  then  calculating  the  corresponding  q  parameter.  The  separations,  Az,  are  1  m.  2.5  m. 
and  5  m  and  A  =  1.1  (Pinkel  and  Anderson,  1992).  It  is  not  clear  whether  averages 
reported  by  Kunze  (1990)  are  from  the  full  duration  of  the  deployment,  however,  good 
model  agreement  can  be  seen  from  the  distributions  on  18-22  October  (Figure  4.7). 
Kunze  (1990)  report  11%  of  the_estimates  of  Fr  at  1  m  spacing  are  greater  than  2  and 
unstable.  From  (4.21)  with  Az  =  1  and  q  =  0.4425,  the  model  estimates  that 
P(Fr  >  2)  =  0.125.  The  model  is  consistent  with  the  observations. 

4.4.  Conclusions 

The  probability  distribution  function  of  Froude  number  is  modeled  using  known 
distributions  of  vertical  strain  and  vertical  shear.  The  derivation  assumes  that  shear  and 
buoyancy  are  independent  variables.  Joint  histograms  of  buoyancy  frequency  and  shear 
demonstrate  that  there  is  small  dependence  on  buoyancy  frequency  of  shear  with  a 
positive  correlation  but  that  the  assumption  of  independence  is  an  appropriate 
approximation.  The  distribution  of  strain  is  Poisson  as  modeled  by  Pinkel  and  Anderson 
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(1991)  and  related  directly  to  the  distribution  of  buoyancy  frequency.  The  distribution  of 
the  absolute  value  of  shear  is  Rayleigh  distribution  with  each  component  of  shear  having 
a  normal  distribution  and  being  independent.  The  resulting  Froude  number  distribution 
has  two  parameters.  The  first  is  X'  -  AAr  from  the  Pinkel  and  Anderson  (1992)  model 
where  X  -  1.1  and  A z  is  the  vertical  scales  over  which  the  Froude  number  is  calculated. 
The  second  parameter,  q,  scales  with  the  square  of  the  ratio  of  the  mean  buoyancy 
frequency  to  the  mean  shear. 
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Figure  4.1  Joint  Histograms  of  Vertical  Shear  and  Buoyancy  Frequency.  Shear  and  buoyancy 
frequency  calculations  are  made  over  10  meter  depth  intervals  and  histograms  represent  nine  25  m.  depth 
bin  ranges.  Normalized  histograms  are  contour  plotted  for  each  depth  range.  Each  histogram  represents 
approximately  50000  estimates  of  shear  and  strain. 
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Figure  4.2  Probaotiity  Distribution  of  Independent  Shear  and  Buoyancy  Frequency. 
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Figure  4.3  Difference  Between  Independent  and  Sampled  Joint  Distributions. 


Table  4.1  Probability  Distribution  Parameters. 
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ixlO--3  3.31x10'  2.43x10-3  2.22xiO'6  0.568  0.167  0.221  0.541 


depth  range  75  to  100m  J  l  depth  range  100  lo  125m  J  L  <kplh  ,,n*e  125  to  iSOm 
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Figure  4.4  Histograms  of  Buoyancy  Frequency.  Buoyancy  frequency  histograms  arc  found  by  inicgrating  joini  histograms.  Buoyancy 
frequency  calculations  arc  made  over  10  meter  depth  intervals  and  represent  nine  25  m.  depth  bin  ranges.  The  solid  curves  arc  corresponding 
model  distributions  from  Pinkcl  and  Anderson  (1991). 


Figure  4.5  Histograms  of  Vertical  Shear.  Shear  histograms  arc  found  by  integrating  joint  histograms  over  buoyanev  frequency  .  Shear 
calculations  arc  made  over  10  meter  depth  intervals  and  represent  nine  25  m.  depth  bin  ranges.  The  solid  curves  arc  die  corresponding  Rayleigh 
probability  distributions. 


75  to  100  ml!  fill  1 1  N  100  to  125  mil  l\  |  (V,  125tol50m 
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Figure  4.6  Histograms  of  Froude  Number  Froude  number  calculations  arc  made  over  10  10  meter  depth  intervals  and  represent  nine  25  m. 
depth  bin  ranges..  The  solid  curves  arc  corresponding  modeled  distribution.  The  model  distribution  assumes  that  shear  and  strain  arc 
independent. 
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Figure  4.7  Distributions  of  Froude  Number  from  a  Neutrally  Buoyant  Float  (Kunze  et  al.;  1990). 
Froude  histograms  are  from  Kunze  et  al.  (1990)  and  were  collected  the  RiNo  float  during  PATCHEX 
(October,  1986).  The  instrument  was  deployed  at  a  depth  of  180-200m.  The  histograms  are  from  18  -22 
Oct.  and  represent  3  vertical  differencing  intervals:  1  m,  2.5  m,  and  5  m,  The  solid  line  overlaid  on  the 
histograms  is  the  modeled  distribution  using  mean  shear  to  buoyancy  frequency  ratio  estimated  from 
Figure  13.b,  Kunze  et  al.  (1990). 
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Appendix  A 


Surface  .Wave.  Processes  PrQsram.(SWAPP)  tom 

The  Surface  Wave  Processes  Program  (SWAPP)  was  conceived  as  an  experiment 
to  study  mixed  layer  dynamics,  Langmuir  circulation  and  surface  wave  growth.  The 
SWAPP  experiment  was  located  at  35°  8.2'  N,  126°  59.0’  W  (Figure  A.l)  where  FLIP 
was  moored  from  February  22  to  March  18,  1990.  The  CSS  PARIZEAU  was  in  the 
region  from  February  26  to  March  18,  1990. 

Principle  investigators  participating  in  the  experiment  included  J.  Smith  (SIO),  R. 
Pinkel  (SIO),  A.  Plueddemann  (WHOI),  R.  Weller  (WHOI),  W.  Crawford  (IOS)  and  D. 
Farmer  (IOS).  Mixed  layer  current  measurements  were  made  from  FLIP  using  two 
strings  of  Vector  Measuring  Current  Meters  (WHOI).  Surface  currents  and  waves  were 
sampled  using  a  surface  scattering  acoustic  Doppler  array  and  a  resistance  wave  wire 
array  (SIO).  Two  down  looking  Doppler  sonar  systems  (161  kHz  and  67/80  kHz) 
mounted  on  FLIP  monitored  the  internal  wave  activity  (SIO).  Density  measurements 
were  made  with  two  continuously  profiling  CTD's  (SIO).  Profiles  of  turbulence  in  the 
mixed  layer  were  made  with  the  Fast  Light  Yo-yo  II  (FLY  II;  IOS)  deployed  from  the 
CSS  PARIZEAU.  Meteorological  measurements  were  taken  from  FLIP  using  a  Vector 
Averaging  Wind  Recorder  (VAWR;  WHOI),  an  Improved  Meteorological  (IMET; 
WHOI)  system  and  sonic  anemometers  (WHOI). 

Meteorological  Data  was  compiled  by  the  group  from  W.  H.  O.  I.  (Figure  A.l). 
Wind  speed  data  was  taken  from  the  mast  of  FLIP.  Wind  velocities  greater  that  10  m/s 
occurred  on  yeardays  60,  64,  67  and  70.  During  the  experiment,  the  sea  surface 
temperature  stayed  between  13°  and  14°C  w’hile  the  air  temperature  ranged  from  9°  to 
14°C.  The  total  heat  flux  was  dominated  by  a  diurnal  cycle  with  daytime  heating 
matched  by  nighttime  cooling. 

Strong  frontal  activity  in  the  upper  ocean  occured  during  the  experiment  . 
Satellite  IR-imaging  of  the  area,  although  hindered  by  cloud  cover,  showed  strong  fronts 
in  surface  temperature  to  the  north  and  south  of  the  SWAPP  location.  An  extensive 
expendable  bathc-rhsrmogram  (XBT)  survey  taken  from  the  USS  NAVAJO  from 
February  20  to  February  28,  1990  also  shows  strong  frontal  activity  in  the  area.  Thirty 
out  of  fifty  XBT  temperature  profiles  yielded  reliable  data  to  400m  depth  and  w'ere  used 
in  the  contour  plots  in  Figure  A. 3.  The  survey  covered  a  100  km  square  area  around 
FLIP.  The  average  temperature  from  the  surface  to  20m  (Figure  A. 3a)  reveals  a  "S” 
shaped  front  with  temperatures  decreasing  1.6°  from  the  northwest  to  the  southeast.  The 
average  temperature  from  the  surface  to  20-l00m  (Figure  A. 3b)  shows  a  similar  shape 
with  temperatures  decreasing  2.0°  from  the  northwest  to  the  southeast.  This  near  surface 
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front  is  not  apparent  at  100-200m  (Figure  A. 3c)  where  temperatures  range  from  9.5°  to 
10.1°.  At  200-400m  (Figure  A. 3d),  a  deep  front  occurs  with  temperatures  decreasing  1.0C 
from  the  southeast  to  the  northwest,  opposing  the  surface  temperature  gradient. 
Temperature-sai’  -*ty  scatter  plots  also  show  evidence  of  two  distinct  water  masses  in  the 
200-400m  r.„^r  the  experiment  site  depth  range  (Figure  3.1). 

The  vertical  density  structure  measured  during  SWAPP  is  typical  of  that 
historically  found  in  the  area  (Figure  A.4).  The  depth  of  the  mixed  layer  is  -100  meters 
which  is  deeper  than  found  in  PATCHEX(Oct.-Nov.,  1986)  and  MILDEX(Oct.,  1983): 
two  upper  ocean  studies  conducted  in  the  area  (Figure  A. 5).  This  is  expected  since 
SWAPP  occurred  in  the  spring  when  the  mixed  layer  has  nearly  reached  its  maximum 
depth.  However,  cruise  averaged  density  profiles  from  MILDEX,  PATCHEX  and 
SWAPP  all  converge  at  about  200  meters  depth. 
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Figure  A.l  The  S^'APP  she.  The  Surface  Waves  Processes  Procram  (SWAP?)  experiment  took  place 
approximately  500  miles  off  the  coast  of  California  southwest  of  San  Francisco.  R/P  FLIP  was  moored  at 
35'  8.2'  N,  1 26c  59.0  W  from  February  22  to  March  18,  1990,  Depth  contours  arc  in  hundreds  of  meters. 


55  60  65  70  75 


lime  (yeardays  1990,  GMT) 

Figure  A.2a-e  SWAPP  Meteorological  Data.  During  the  SWAPP  experiment,  the  meteorological  data 
was  collected  by  the  W.  H.  0.  1.  croup  led  by  R.  Weller.  Wind  speed  (a)  and  wind  direction  (b)  were 
collected  with  the  Vector  Averaging  Wind  Recorder  (VAWR).  Barometric  pressure  is  shown  in  (c).  Sea 
Surface  temperature  (d;  dashed  line)  is  from  the  VMCM  at  2.25m  depth  and  is  ploued  along  with  air 
temperature  (d;  solid  line).  The  total  heat  flux  (e)  is  the  sum  of  the  estimated  short-wave  radiation,  long¬ 
wave  radiation,  sensible  heat  flux  and  latent  heat  flux. 
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Figure  A.3a-d  Depth  Averaged  Temperature  Contours.  From  February  20  lo  February  25  1990, 
Weller’s  Group  from  W.  H.  0. 1.  performed  an  X3T  survey  from  the  Navajo.  From  57  temperature 
profiles  taken,  30  art  selected  here  which  appear  reliable  down  to  -00  meters.  For  the  contour  plots, 
temperature  profiles  from  the  XBT  survey  are  averaged  in  four  depth  ranges:  (a)  surface  to  20  meters,  (b) 
20  to  200  meters,  (c)  100  to  200  meters  and  (d)  200  to  -00  meters.  These  depth  averaged  temperatures 
are  contoured  with  labels  in  degrees  Celsius.  The  survey  reveals  a  strong  front  Near  surface  water 
temperature  increases  by  2®  C  from  southeast  to  northwest,  while  the  water  temperature  below  200  meters 
decreases  by  1®  C. 
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0 


Potential  Density  Buoyancy  Frequency  (C.P.H.) 


Figure  A. 4  Historic  Hydrographic  data  and  Cruise  Averaged  Potential  Density.  The  historic 
hydrographic  data  is  from  J.  Ried  s  hydrosearch  on  the  computer  nemo’  which  contains  a  database  of 
physical-chemical  oceanographic  stations  from  around  the  world.  25  stations  were  found  that  occurred 
during  February  and  March  between  the  latitudes  34®  N  and  56°  N  and  the  longitudes  126®  W  and  128® 
W.  The  data  spans  34  years  from  1950  to  1984,  Most  of  the  stations  extend  to  600-1 500m  depth  while  4 
stations  contain  data  to  full  ocean  depths.  Plotted  here  are  the  potential  densities  and  buoyancy 
frequencies  from  the  hydrographic  data  (dots)  and  the  SWAPP  cruise  averaged  profiles  (solid  line)’.  The 
SWAPP  data  falls  within  the  scatter  of  the  historic  data. 
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Figure  A.5  Cruise  Averaged  Poieniial  Dcr.siry.  Th;  MILDEX  data  (dashed)  was  taken  for  14  days 
during  October/November  1983  while  FLIP  drifted  near  34'N,  126CW.  The  PATCHEX  data  (dotted) 
was  taken  for  21  days  during  October  1986  while  moored  at  34 °N,  127*“W.  The  SWAPP  data  (solid)  was 
taken  for  20  days  during  February/Marth  1990  while  moored  at  35°N,  127CW.  Cruise  averaged  potential 
density  is  plotted  versus  depth.  Note  that  MILDEX  and  PATCHEX  both  took  place  during  the  fall  and 
have  a  much  shallower  mixed  layer  depth  than  SWAPP  which  took  place  in  late  winter. 
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Appendix  B 


B.l.  Equations  of  Motion. 

In  this  section,  a  review  of  linear  internal  gravity  wave  theory  is  presented.  This 
will  give  a  general  framework  for  linear  internal  waves  used  in  the  discussion  of  the 
measured  internal  wave  field. 


It  will  be  assumed  the  flow  field  is  incompressible,  inviscid,  and  continuously 
stratified.  The  continuity  equation  in  this  case  is 


du  dv  dw 
dx  dy  dz 


(B.l) 


where  u,  v  and  w  are  the  components  for  the  fluid  velocity  in  the  x,  y  and  z  directions 
respectively.  In  a  rotating  fluid,  the  linearized  momentum  equations  for  small 
perturbations  in  pressure,  p\  and  density,  p',  are 


where  p(z)  is  the  unperturbed  density,  g  is  the  acceleration  of  gravity  and /is  the  Coriolis 
parameter.  The  linearized  pressure  equation  is 

~r~  +  w~r~  ~  0  (B.3) 

dt  dz 

By  making  the  Boussinesq  approximation  (w  varies  much  more  rapidly  in  the  vertical 
than  pj,  linear  traveling  plane  wave  solutions  to  the  above  equations  of  motion  can  be 
found.  The  solutions  are  of  the  form 
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l&L  N2W  gi(-k**>y+’u~aM'1 
geo 


(B.4) 


with  the  linear  dispersion  relation 

(/w +  !¥»(*»+;»)) 

(k’+l1  +m2) 


for  f  <  CD  <  N. 


(B.5) 


This  highest  frequency  of  oscillation  is  the  Brunt-Vaisala  frequency  or  buoyancy 
frequency,  N,  and  is  defined  as 


N(z)2  ~  -gp~] 


dp(2) 

dz 


(B.6) 


B.2.  Counter  Clockwise  to  Clockwise  Rotating  Shear  Ratio. 

Consider  an  internal  wave  traveling  in  the  positive  x  direction  (/=0),  the  phase  is 
propagating  in  the  positive  z  direction  and  energy  is  propagating  in  the  negative  z 
direction.  If  z=0  at  the  ocean  surface  and  z  increases  with  increasing  depth,  then  energy 
is  traveling  down  from  the  surface.  Since  the  interest  here  is  to  study  the  velocity  field  in 
time  and  depth,  the  horizontal  coordinate,  x ,  is  set  equal  to  zero.  The  Eulerian  velocity 
field  for  this  internal  wave  is  defined  by 


w  =  Wc  cos{mz  -  cot) 


u  -  Wocos{mz-  ox) 

k 

v  =  sin(mz-  cot) 

k  co 

p'  -  --P—N2WBs\n{mz- cot) 
geo 


(B.7) 


where  the  internal  wave  dispersion  relation  is 


m 

k 


' co2-N 2^ 
f2~co2 


for/  <  co  <  N, 


(B.8) 
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Solving  for  the  vertical  shear  field  for  this  internal  wave  yields. 


u.  = 


dii  +  ,dv 
dz  dz 


= 


r  co2  -  N2  ^ 


co- 


n- 


sin(mz  -  cot )  -  i  —  cos(mz  -  cot) 
co 


(B.9) 


Solving  for  the  vertical  shear  spectrum  obtains 
<t>u  ( m,cb)  = 

co2 -N 


W0  m 


f2-co2 


1  +  — )  S(m  +  m)5(cb-  co)  +  (l~  — 1  6(m-  m)8(cb  +  co) 
co)  V  co) 


(B.10) 


The  counter  clockwise  (CCW)  to  clockwise  (CW)  shear  variance  for  this  case  is 


ccw  _  ( co- 
cw  \co  +  f  j 


(B.ll) 


For  the  high  frequency  waves,  co»  f ,  spectral  energy  is  found  equally  at  both  (+m,-co)  and 
(-m,+co).  As  co  approaches  the  inertial  frequency,  however,  energy  is  only  found  only  at 
(- m,+co ).  Extending  this  to  a  full  spectrum  of  upward  and  downward  propagating  internal 
gravity  waves,  near  inertial  energy  is  found  only  in  the  (+m,+co)  (counterclockwise, 
upward  propagating  phase,  downward  energy)  quadrant  and  the  ( +m,-co )  (clockwise, 
downward  propagating  phase,  upward  energy)  quadrant.  High  frequency  waves  fill 
energy  throughout  all  the  quadrants  in  ( m,co )  space. 


B.3.  Shear  Variance  to  Strain  Variance  Ratio 


Strain  is  defined  here  as  the  vertical  derivative  of  vertical  displacement.  The 
vertical  displacement,  T),  can  be  found  by  integrating  the  vertical  velocity  field  to  yield 


n 


(B.12) 


From  the  linearized  pressure  equation,  rj  is  found  in  terms  of  the  perturbation  density 


r)(2,t) 


P'(z.t) 

d'Po 


(B.  13) 


The  Eulerian  strain  is  the  vertical  derivative  of  the  displacement  field  with  respect  to  the 
vertical  coordinate,  z  ,  and  is  defined  as 


(B.  14) 


where  the  density  field  is  determined  by  a  mean  state  and  a  fluctuating  component  due  to 
internal  waves,  p(z,t)  =  p'  (z,t)  +  p(z).  The  Eulerian  strain  is  calculated  by  first 
differencing  density  in  depth,  normalizing  each  profile  by  the  first  differenced  mean 
density  profile  and  subtracting  1.0.  The  Semi-Lagrangian  strain  is  the  derivative  of  the 
displacement  field  with  respect  to  the  Semi-Lagrangian  coordinate 

C  =  z~Tl  (B.15) 


and  is  calculated  by  taking  the  instantaneous  difference  between  isopycnal  displacement 
and  dividing  by  the  mean  depth 


t,  -iQ- 
V<~d£ 


ac 


AC 


AC 


(B.16) 


The  Semi-Lagrangian  strain  is  related  to  the  Eulerian  strain  by 

„  Jg.V2i.JL- 

k  UzJ  dz  1-77,' 


(B.17) 


Linear  internal  wave  theory  (using  equation  B.4b-d)  predicts  a  buoyancy  frequency 
normalized  strain  variance  to  shear  variance  ratio  of 


fl/2  j^zl  _ 

<!>„  (to) 


N 


mJ2  N2  -  to2  {  to2  +  / 


to2-f2 


(B.18) 


for  a  spectrum  of  linearly  superimposed  waves. 


105 


Appendix  C 


CTP  processing 

Fine  structure  density  estimates  are  made  using  conductivity,  temperature  and 
pressure  data  collected  with  two  repeat  profiling  CTDs.  Sea-Bird  Electronics  (SBE) 
instruments  were  deployed  during  SWAPP  experiment  in  a  similar  arrangement  as  has 
been  used  in  past  upper  ocean  experiments  on  FLIP  (Williams,  1983;  Sherman,  1989). 
The  conductivity  sensors  are  model  SED  4-01  (5X10‘5  S/m  resolution  at  12Hz),  and  the 
temperature  sensors  are  model  SBE  3-0 1/F  (5X1  O'4  °C  resolution  at  12Hz).  The  CTD 
sensors  and  a  pressure  case  containing  sampling  electronics  are  mounted  in  a  streamlined 
cage  with  a  55  kg  lead  weight  in  front  to  allow  for  high  fall  rates.  During  the  SWAPP 
deployment,  a  SBE  dual  needle  conductivity  probe  was  placed  in  front  of  each 
instrument.  The  shape  of  the  cage  allows  profiling  only  in  one  direction  since  all  the 
sensors  are  in  the  instrument  wake  while  being  pulled  to  the  surface.  The  cycle  period 
was  130  seconds  from  the  start  of  one  profile  to  the  beginning  of  the  next  and  fall  rates 
were  3.65  m/s.  The  upper  CTD  profiled  from  5-220  meters  and  the  lower  CTD  profiled 
rrom  200-420  meters.  The  temperature,  pressure  and  conductivity  sensors  are  sampled  at 
24  Hz  and  the  dual  needle  probe  is  sample  at  96  Hz.  The  upper  CTD  was  operated  from 
February  26  to  March  18,  1990  collecting  -12000  profiles  and  the  lower  CTD  was 
operated  from  March  1  to  March  16,  1990  collecting  -9200  profiles. 

A  pitot  tube  equipped  pressure  sensor,  used  to  measure  the  static  pressure  during 
profiling  and  is  used  to  estimate  depths,  z.  A  typical  profile  of  first  differenced  pressure 
shows  a  high  frequency  fluctuation  associated  with  water  turbulence  (Figure  C.l).  Using 
a  simple  transformation  from  frequency  to  wavenumber  assuming  constant  fall  rate,  a 
spectrum  for  the  first  differenced  pressure  reveals  a  +1  power  law  above  0.1  cpm.  This  is 
likely  due  to  water  turbulence  and  a  does  not  reflect  real  fluctuations  in  fall  rate.  The 
pressure  signal  is  low  pass  filtered  below  0.14  cpm  to  remove  the  turbulent  noise.  Note 
the  spectral  peak  at  0.022  cpm  corresponding  to  a  pressure  signal  from  a  surface  wave 
swell  period  of  12.5  seconds.  The  pressure  data  was  then  resampled  at  half  the  data  rate. 

Following  the  procedure  developed  by  Williams  (1985)  with  one  slight 
modification,  the  time  response  of  the  temperature  (T)  and  conductivity  (C)  sensors  are 
matched.  The  method  assumes  that  the  conductivity  fluctuations  at  scales  smaller  than  10 
meters  are  due  mainly  to  temperature  fluctuations.  By  computing  the  cross  spectrum 
between  temperature  and  conductivity,  a  transfer  function  phase  and  amplitude  function 
can  be  estimated.  The  Fourier  transform  of  measured  temperature,  tm,  and  conductivity, 
cm  ,  are  related  to  the  true  values  by 
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(C.l) 


tm(m)^t(rn)GT(m)e‘*'im) 
cm(m)  =  c(rn)Gc(m)e‘*c(m] 

where  m  is  the  vertical  wave  number,  Gt>  Gq  ,4>t  and  are  the  amplitude  and  phase 
responses  of  the  sensors.  The  temperature  spectrum  is  related  to  the  cross  spectrum 
(assuming  no  salinity  contribution)  by 

(O.)0/  =  a(,;c.)GrOce (C2) 

,  dl 

where  &  =  — — . 

aC 

The  transfer  function  is  found  by 

aV«  w  Gt  (C.3) 

H(m)  is  estimated  using  100  profiles  from  the  upper  and  lower  CTD’s.  The  profiles  are 
first  differenced  then  Fourier  transformed  to  compute  temperature  and  cross  spectra.  The 
transfer  function  is  modeled  by  fitting  a  fifth  order  polynomial  to  the  estimated  amplitude 
and  phase  (Figure  C.2).  All  profiles  of  T  and  C  are  first  differenced  and  Fourier 
transformed.  The  transfer  function  is  applied  to  match  the  spectral  response  of  T  to  that 
of  C.  Both  T  and  C  are  then  low  pass  filtered  with  a  cut  off  of  1.3  c.p.m.  before  inverse 
Fourier  transforming.  T  and  C  are  then  reintegrated  and  sampled  at  half  the  data  rate  to 
form  corrected  profiles.  Salinity,  potential  temperature  and  potential  density  are 
estimated  using  the  equation  of  state  (Fofonoff  and  Millard,  1983).  Examples  of  raw  and 
corrected  density  profiles  show  that  much  of  the  density  overturns  from  sensor  mismatch 
is  removed  (Figures  C.3  and  C.4). 

Upper  CTD  and  lower  CTD  density  profiles  are  combined  using  the  overlap 
depths  of  200-220  meters.  A  linear  ramp  is  using  in  the  region  such  that  at  205m  100% 
of  the  upper  CTD  data  is  used,  at  210m,  50%  for  each  profile  is  used  and  at  215  m  100% 
of  the  lower  CTD  profiles  is  used.  It  should  be  noted  that  there  is  a  60  sec  time  lag 
between  the  sampling  of  this  region  by  the  lower  CTD  and  sampling  by  the  upper  CTD. 
This  time  lag  will  bias  strain  variance.  As  an  illustration,  assume  vertical  advection 
moves  a  constant  density  gradient  with  no  straining.  The  density  field  will  be  displaced 
between  sampling  by  the  lower  and  upper  CTD.  Advection  towards  the  surface  will 
appear  as  a  positive  strain  and  advection  towards  the  bottom  will  appear  as  negative 
strain. 
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pressure  signal  power  spectra  (tlaslwal  line.) 


Figure  C.2  Temperature  and  Conductivity  response  matching.  Transfer  functions  (solid  lines)  between 
temperature  and  conductivity  are  estimated  from  500  profiles  from  both  the  upper  and  lower  CTD.  A  fifth 
order  polynomial  is  fit  to  the  amplitude  and  phase  of  the  estimated  transfer  function  and  is  used  to  match 
the  response  of  the  two  sensors  prior  to  salinity  calculations. 
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(offset  .1  kg/mA3) 


Figure  C3  Vncorrected  Der.siry  Profiles.  Five  Density  profiles,  calculated  from  raw  conductivity  and 
temperture  data,  from  the  upper  CTD  are  plotted.  A  mismatch  between  the  spatial  responses  of  the 
conductivity  cell  and  temperature  probe  causes  'salinitv  spiking'  leading  to  poor  estimates  of  potential 
d. 


Appendix  D 


Velocity  Biasing  Due  to  Horizontal  Beam  Separation 

Horizontal  velocity  from  Doppler  sonar  beams  are  obtained  by  comparing  slant 
velocites  from  back  to  back  beams.  This  velocity  is  only  a  good  estimate  of  the  true 
horizontal  velocity  if  the  horizontal  wave  scales  are  much  longer  than  the  horizontal 
beam  separation  (Sherman  and  Pinkel,  1991).  The  section  estimates  the  expected  spectral 
error  as  a  function  of  vertical  wavenumber  and  frequency.  From  linear  internal  wave 
theory,  the  vertical  and  east  velocity  components  for  a  single  internal  wave  are 

w  = 

u  =  -  -(cos  y  +  i  —  sin  y)W0e**1*t**"a-m)  (°' 1 } 

K  co 


where  the  horizontal  wave  numbers  k=Kc os(y)  and  /=^sin(y).  Assume  that  the  sonar  is 
located  at  the  coordinates  (x,y,z)=(0,0,0).  This  wave  is  observed  by  the  east  and  west 
beams  of  the  sonar  in  the  x-z  plane  with  a  beam  down  look  angle  of  8 .  The  slam  velocity 
is  found  by  taking  the  projection  of  the  velocity  field  along  each  beam.  The  horizontal 
locations  for  each  beam  as  a  function  of  depth  are  xt  =  z  cot  0  and  xw  -  -2  cox  8.  The 
measured  slant  velocities  are  then 


vw  =  W0e,(mI-“) 


cos0  +— (cosy  +  /— sin  y)sin0 
K  co 

Tn  f 

cosd - (cosy  +  t  —  sin  y)sin0 

K  co 


iKx  cot  yco \  8 


—  iKt  cos  y  cot  9 


(D.2) 


The  estimated  horizontal  velocity  is  obtained  by 


u  —  —  •  =  -W0e‘(mJ~cu,) 

2  sin  6 


^^cos  y  +  i  '"sin  y  jcos0  +  i  cot  6  sin  <p 


(D.3) 


where  <P  =  Kzcosycoxd 

For  a  wave  propagating  in  the  north-  south  plane  (y=±tc/2),  the  estimated  east  velocity 
equals  the  true  east  velocity.  For  a  wave  moving  in  the  east-west  plane,  there  is  a  phase 
and  amplitude  shift  in  the  estimated  velocity  field.  This  error  is  small  for  waves  with 
horizontal  wavelengths  that  are  much  larger  than  the  beam  spacing.  The  biasing  of 
spectral  estimates  of  the  velocity  field  is  assessed  by  comparing  the  measured  variance  to 
the  expected  variance  The  bias  at  depth,  z,  for  a  single  internal  wave  is  given  by 
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The  total  spectral  bias  for  a  horizontally  isotropic  wave  field  is  determined  by  averaging  b 
over  the  vertical  length  of  the  sonar  beams  and  over  all  wave  angles,  y.  Averaging 
Equation  D.4  over  z  =  1  to  256  meters  and  y=0  to  27t,  and  using  N  =  6  c.p.h.  and /=.0478 
c.p.h.  yields  Figure  (D.l),  the  total  spectral  velocity  bias  as  a  function  of  vertical 
wavenumber,  m,  and  frequency,  0).  The  total  spectral  bias  has  a  minimum  of  0.5.  There 
is  very  little  bias  found  at  low  frequencies  (a><2f)  for  all  measurable  wave  numbers.  As 
wave  frequencies  approach  N,  the  waves  have  strong  vertical  velocities  and  short 
horizontal  wavelengths  leading  to  an  overestimation  of  the  horizontal  velocity  variance. 
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Figure  D.l  Estimated  Spectral  variance  Bias  Due  to  Beam  Pair  Separation.  Horizontal  velocities 
determined  by  comparing  back  to  back  beams  art  biased  estimates  of  the  true  velocity  field  due  to 
horizontal  beam  separation.  Assuming  a  horizontally  isotropic  internal  wave  field,  the  bias  is  averaeed 
over  all  directions  and  over  the  vertical  extent  of  the  sonar  measurements  (z=0  to  256  meters;  >=.0478 
c.p.h.;  N=6  c.p.h.).  At  low  frequencies  (2f<cuc2/),  waves  have  long  horizontal  wavelengths  and  biasing 
is  negb'gible.  As  9  approaches  N,  waves  have  strong  vertical  velocities  and  short  horizontal  wavelengths 
and  variances  is  overestimated.  Variance  is  underestimated  by  as  much  as  50%  at  hich  vertical 
wavenumbers. 
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Appendix  E 


SflpaLgaia^Eflcfissiae 


E.l.  Doppler  sonar  background 


To  estimate  water  velocities,  Doppler  sonar  systems  transmit  a  pulse  of  sound  and 
measure  the  Doppler  frequency  shift  of  the  reverberation  produced  by  acoustic  scatterers 
in  a  volume  of  water.  The  acoustic  scatterers  in  the  upper  ocean  are  biological  and 
comprised  mostly  of  zooplankton.  To  the  extent  that  the  scatterers  are  advected  passively 
by  the  surrounding  water,  the  frequency  shift  reflects  the  mean  velocity  of  the 
interrogated  volume  of  water.  For  a  monostatic  Doppler  sonar,  the  sonar  transducers  are 
both  the  source  and  receiver.  Assuming  that  the  relative  velocities  between  the  water  and 
the  sonar  are  much  smaller  than  the  speed  of  sound,  c,  the  reverberation  frequency  at  the 
receiver,  /r,  is  determined  by. 


1  +  ~V 


(E.l) 


where /c  is  the  transmit  frequency  and  V  is  the  relative  velocity  of  the  interrogated  water. 
Rewriting  (E.  1)  to  find  the  relative  velocity  in  terms  of  the  Doppler  shift,  fd  -frfc  obtains 

V  =  (E.2) 


where  V  is  positive  if  the  volume  is  moving  towards  the  sonar.  For  a  simple  Doppler 
system,  a  volume  is  interrogated  with  a  pulsed  transmission  of  the  form 

X0(r)  =  Acos(27T/cr)G(t)  (E.3) 


where 


if  0<t<Tp 
otherwise 


and  Tp  is  the  pulse  duration.  The  received  signal  at  time,  /,  is  defined  as 

X(r)  =  Re{Z(t)elJ"/‘1} 


(E.4) 
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where  Z(t)  is  the  complex  envelope.  (X(r)  is  base  banded  with  a  reference  signal  and 
with  the  reference  signal  shifted  in  phase  by  90°  to  remove  the  carrier  frequency.)  Then 
the  complex  envelope  is  sampled  with  two  A/D  conveners  with  a  sample  interval  of  x. 
This  results  in  the  sampled  complex  envelope,  Z(rn).  where  trn  =/  is  the  sample  number 
and  tn  is  zero  at  the  start  of  the  acoustic  transmission  sequence.  To  estimate  a  mean 
Doppler  shift  from  the  spectrum  of  shifts  returned  by  various  scatterers  in  the  interrogated 
volume,  a  power  weighted  mean  frequency  of  the  complex  envelope  is  found.  This  is 
done  by  estimating  autocovariances  using  a  technique  developed  by  Rummler  (1968). 
The  autocovariance,  C,  is  the  product  of  the  signal  at  sample  rn  with  its  self  at  sample 
tn+L,  where  L  is  the  sample  lag  number.  The  autocovariance  is  then  averaged  in  time  to 
increase  the  stability  of  the  autocovariance  estimate.  The  averaging  over  consecutive 
autocovariances  in  time  is  called  range  averaging  and  the  averaging  intervals  are  called 
range  bins.  The  average  covariance  over  a  range  bin  is  found  by 

,  sib,* o-i 

cuo=-j-  I  zd.)z’(i.+i)  (E.5) 


where  N  is  the  averaging  interval  and  bn  is  the  range  bin  number.  This  step  is  commonly 
referred  to  as  range  averaging.  To  increase  the  stability  of  the  autocovariance  estimates, 
several  realizations  from  successive  transmissions  are  averaged  together  which  results  in 
time  averaging.  Rummler  showed  that  the  estimated  covariance  is  related  to  the  power 
weighted  mean  frequency  of  the  sampled  complex  envelope  by 


f  (*„)« 


arg{C(fr,)) 

2nLx 


(E.6) 


In  a  noise  free  sonar  with  a  symmetric  returned  Doppler  frequency  spectrum,  this  is  an 
unbiased  estimate  of  the  mean  Doppler  shift.  In  practice,  for  time  lags  much  less  than  the 
inverse  bandwidth,  f  is  a  good  estimate  of  the  mean  Doppler  shift.  Slant  velocities  are 
then  found  from  f  using  (E.2). 

The  optimum  performance  of  a  simple  Doppler  sonar  in  the  limit  of  an  infinite 
signal  to  noise  ratio  can  be  estimated  using  the  Cramer-Rao  bound  (Theriault,  1986) 

c2 

AVAR  =  — ±-=  (E.7) 

8  nfc4P 


where  AV  is  the  rms  velocity  imprecision,  AR  is  the  effective  range  resolution,  and  P  is 
the  number  of  independent  transmit  sequences  averaged  together.  The  effective  range 
resolution  is  half  the  pulse  length:  AR  =  %cTp . 

Recently,  repeat- sequence  transmit  coding  has  been  used  to  increase  the  precision 
of  Doppler  sonars  (Pinkel  and  Smith,  1992).  Repeat-sequence  codes  are  comprised  of 
several  repeats  of  a  broad  band  subcode.  The  subcode  is  selected  to  have  the  smallest 
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possible  autocorrelation  except  at  zero  lags  and  reverses  the  sign  of  the  carrier  signal  at 
regulated  time  increments.  In  this  case,  G(r),  the  transmit  envelope,  is  no  longer  constant 
but  takes  on  a  value  of  +1  or  -1.  The  subcode  is  made  up  of  L  bits  each  lasting  a  period 
length  of  x.  The  subcode  is  then  repeated  M  times.  An  example  of  four  repeats  of  a  7-bit 
Barker  code  is  shown  in  Figure  E.l.  An  ideal  code  is  one  that  results  in  an 
autocovariance  of  nearly  zero  everywhere  except  at  sample  lags  L.  By  transmitting  codes 
instead  of  discrete  pulses,  the  average  transmit  power  is  increased.  Corresponding 
autocovariances  are  calculated  at  lags  equal  to  the  code  length.  Pinkel  and  Smith  (1992) 
found  performance  is  limited  by 


A  VAR  = 


JA 


AnP 


Vi 


m. 


T 

1  mix 

V  row  A 


7  1 

1  0*1  +  _ 

T  2 


VS 


(E.8) 


where 

c 

-  ~~~  »  maximum  unambiguous  velocity  estimate, 

M  s  number  of  repeats  of  the  subcode, 

L  s  number  of  bits  in  the  subcode, 

Tmt  mNv, 

T0vi  a  (M  -  \)Lx ,  overlapped  common  range  in  the  covariance  estimate, 

/V  s—  ,  maximum  bandwidth. 

T 

f  =  the  transmission  frequency. 

Tm.*  =  the  max  of  (Tavg,Talv). 

Thus  the  rms  precision  of  the  sonar  is  enhanced  roughly  in  proportion  to  the  product  of 
the  bandwidth  and  the  duration  of  the  subcode.  These  coded  pulses  have  now  been  used 
on  several  Doppler  sonars  at  the  Marine  Physical  Laboratory.  Precision  tests  agree 
roughly  with  the  predictions  of  (E.8). 

E.2,  The  MPL  161  kHz  Doppler  sonar  system. 

The  MPL  161  kHz  Doppler  sonar  system,  designed  and  built  originally  for  the 
CEAREX  (1989)  experiment,  was  mounted  on  the  thruster  of  FLIP  at  a  depth  of 
approximately  15  meters  for  the  SWAPP  experiment.  The  system  configuration 
parameters  used  for  the  FLIP  deployment  are  listed  in  Table  E.l.  The  system  used  a 
traditional  four  beam  arrangement  where  the  beams  were  spaced  90°  apart  horizontally 
(Figure  E.2).  Each  beam  pointed  down  with  a  downlook  angle,  6q  ,  of  60°  from  the 
surface.  The  slant  velocity  for  each  beam  is  the  projection  of  the  true  velocity  field  onto 
that  beam.  The  estimated  horizontal  and  vertical  velocities  are  found  by  comparing  back 
to  back  beams  such  that 
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U02  “ 


2cos0n 


U.3  = 


(^~n) 

2  cos  dn 
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02 


■  i^and  w„  =  <^) 


2  sin  8n 


2  sin  8n 


(E.9) 


where  V0,  V;,  V2,  and  V?  are  the  slant  velocity  estimates  from  the  four  corresponding 
beams;  U02  and  u  13  are  the  orthogonal  horizontal  velocities  and  wq2  and  W13  are  the 
vertical  velocities.  Implicit  in  (E.9)  is  the  assumption  that  the  horizontal  scales  of  the 
flow  field  are  larger  than  the  horizontal  spacing  between  the  sonar  beams  (Plueddemann, 
1987).  This  beam  configuration  is  typically  called  the  Janus  configuration  and  the 
estimated  velocities  from  (E.9)  are  referred  to  as  the  Janus  velocities. 

E.3.  Variable  Gain  correction 

To  extend  the  dynamic  range  of  the  analog  to  digital  converters,  a  variable  gain 
circuit  is  used.  The  analog  signal  is  sampled  every  250  fisec  by  a  12  bit  A/D  converter 
that  has  a  dynamic  range  of  72  dB.  Without  any  additional  dynamic  range,  the  sonar's 
potential  profiling  range  would  be  limited  by  the  decay  in  time  of  the  acoustic  return. 
The  A/D  converter  samples  the  analog  signal  after  the  variable  gain  circuit  actively 
amplifies  the  true  receive  signal  This  extends  the  effective  dynamic  range  to  102  dB. 
The  profiling  range  is  thus  limited  by  system  electronic  noise.  Here,  a  model  of  the 
variable  gain  circuit  is  used  to  correct  the  sampled  covariance  and  intensity  amplitudes. 

The  sampled  acoustic  return  data,  Zs(rn),  is  related  to  the  input  signal,  Zo(rn)>  by 

ZA)  =  G%.)Z„('.)  (E.10) 

where  G(/n)  is  the  variable  gain  function  and  fn  is  the  sample  number  from  the  start  of  the 
transmit  sequence.  The  timing  of  the  variable  gain  circuit  is  locked  to  the  transmit  pulse 
so  that  G(tn)  is  applied  identically  to  every  received  sequence.  For  the  first  60  samples 
(15000  ^sec),  G(/n)  is  equal  to  one.  For  later  samples,  the  gain  is  ramped  up  to  30  dB  to  a 
final  value  of  1000,  with  60  steps  of  0.5  dB  (Figure  E.3).  Each  step  lasts  20  samples 
(5000  fisec).  The  recorded  intensity  data,  Ins,  which  is  range  averaged  over  9  samples 
and  is  related  to  the  input  signal  by 

IgV,>z0<oz;(i.)  (e.u) 

w  1. = Kb, 

where  bn  is  the  range  bin  number  and  N  is  the  range  averaging.  The  true  intensity  data, 
Ino(tn),  is 
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(E.12) 


i  Nib.*  !)-l 

Ino(^)  =  T7  E  ZoOMlO 


If  G(tn)  is  a  constant  over  each  range  averaging  interval  then  the  true  intensity  is  simply 
related  to  the  sampled  intensity  by 


In  0(bM) 


In M 

G(6J 


(E.13) 


where 


G(bK)  = 


\_ 

N 


N(b.*  !)-l 


EG^>- 


(E.  14) 


The  ideal  gain  would  be  constant  over  a  given  range  bin.  If  there  is  a  gain  change  over  an 
averaging  interval,  there  will  be  some  biasing  of  the  intensity  that  cannot  be  removed  in 
post  processing.  (In  the  case  of  the  161  kHz  sonar  discussed  here,  the  transmit  pulse 
effectively  smoothes  the  intensity  profiles  over  a  range  bin.  It  is  reasonable  to  assume 
that  this  biasing  error  will  be  small  since  Zo(tn)  w'ill  be  nearly  constant  over  a  range 
averaging  interval.)  Here,  the  biasing  errors  are  neglected  and  the  modeled  variable  gain 
function  and  estimated  true  intensities  are  found  using  (E.13)  and  (E.14).  Figure  E.4 
shows  a  typical  one  hour  averaged  intensity  profile  with  and  without  the  variable  gain 
correction  and  the  full  dynamic  range  of  nearly  80  dB  has  been  restored. 

One  side  affect  of  the  steps  in  the  variable  gain  circuit  is  the  addition  of  an 
intensity  ripple  due  to  the  mismatch  of  the  range  averaging  and  gain  step  lengths.  Figure 
E.5b  shows  a  first  differenced,  hour  averaged  intensity  profile.  During  the  period  of 
active  gain  change  (bins  7  to  122)  a  high  frequency  ripple  occurs  even  after  averaging  for 
one  hour.  This  ripple  can  also  be  seen  in  the  model  variable  gain  function  plotted  in 
Figure  E.5a.  The  same  profile  now  corrected  for  the  variable  gain  is  plotted  in  Figure 
E.5c  and  shows  much  less  variance. 

The  amplitude  of  the  estimated  covariance  data  is  also  affected  by  the  variable 
gain  function.  The  recorded  covariance  data  is  averaged  over  9  samples  and  is  related  to 
the  input  signal  by 

i  A'O-.+l)-! 

c,(f>.)  =  -  £(G(i„)G(r,  +  i))  Z,(r.)Z;(i.U)  (E.15) 

W  ,.*Nb. 


where  L  is  the  lag  number  in  samples.  The  true  covariance  data  is  defined  as 
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(E.16) 


,  N{b.*\)-\ 

C0(bK)  =  j-  X  Z0(r.)z;(r.+L)- 

N  '.=Nk. 


If  G  constant  and  equal  over  each  range  averaging  interval  then  the  true  covariance  is 
obtained  by 


C  0(bJ 


CM 

G(fe.) 


(E.17) 


where 

,  N(b.+\)-\ 

G(fc,)  =  -i  X  (G(r.)G(r.  +  /.))»■  (E.18) 

M  t,*Nb, 

In  the  case  where  G  varies  over  a  range  averaging  interval,  there  will  be  biasing  of  the 
covariance  data  that  cannot  be  removed  with  post  processing.  For  the  161  kHz  sonar,  a  L 
is  7  samples  and  a  gain  step  length  lasts  20  samples  so  only  one  step  could  occur  over  the 
averaging  interval  of  9  samples.  Since  the  gain  steps  are  small  and  the  pulse  length  has 
smoothed  Zo(tn),  it  is  assumed  that  the  biasing  will  be  small  and  the  gain  biasing  errors 
will  be  neglected.  The  modeled  variable  gain  function  for  the  covariance  and  the 
estimated  true  covariance  are  found  from  (E.17)  and  (E.18).  Figure  E.6  shows  a  typical 
profile  of  both  sampled  and  estimated  true  covariance  amplitude. 

E.4.  Noise  correction 

The  correction  applied  here  assumes  that  band  limited  additive  noise  in  the 
receiver  causes  biasing  of  the  covariance  estimates.  The  measured  covariance  is  modeled 
as  the  sum  of  the  signal  covariance  and  the  noise  covariance, 

C„(^)  =  Cl(&.)-fCn  (E.19) 

Estimates  of  f  from  the  measured  covariance  will  be  biased  further  away  from  the  true 
mean  Doppler  shift  as  the  signal  to  noise  ratio  gets  small.  The  bias  can  be  removed  by 
estimating  the  true  signal  covariance.  This  is  done  by  determining  Cn  from  the  return  at 
far  range  where  the  signal  has  decayed  away  leaving  only  noise.  The  noise  covariance  is 
assumed  to  be  constant  over  all  ranges  and  is  then  subtracted  from  the  measured 
covariance  to  find  the  estimated  signal  covariance.  Then  an  unbiased  estimate  of  the 
mean  Doppler  shift  is  found  from  the  corrected  covariance  estimates.  In  practice, 
uncorrected  slant  velocity  profiles  will  tend  towa.’  Is  a  non  zero  velocity  at  far  range.  (To 
remove  this  bias  from  the  161  kHz  sonar  data,  the  average  covariance  from  the  last  9 
range  bins  for  every  profile  is  used  as  an  estimate  of  the  noise  covariance,  Cn,  for  each 
beam.  The  noise  covariance  is  subtracted  from  the  one  minute  averaged  signal 
covariance  prior  to  calculating  slant  velocities.) 
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shape  as  the  wire  spectrum  in  the  surface  wave  band.  However,  at  frequencies  below  the 
swell  peak,  the  energy  density  increases  towards  lower  frequency  with  a  higher 
proportional  energy  than  the  corresponding  wire  spectrum.  It  is  reasonable  to  assume  that 
below  .0167  Hz  most  of  the  accelerations  are  due  to  FLIP  tilting.  Tilt  estimates  used  in 
the  sonar  processing  are  from  1  minute  averaged  accelerometer  data. 

E.5.2.  Transformation  Matrix 


There  are  three  right-hand  coordinate  systems  to  consider.  The  first  is  the 
geovertical  coordinates  where  k  points  vertically  upward  parallel  to  gravity,  j  points 
horizontally  towards  the  bow  of  FLIP  and  i  points  horizontally  towards  the  starboard  of 
FLIP.  The  other  coordinate  system  is  the  ship  coordinates  and  is  defined  by  the  unit 
vectors  i,  j,  and  k  in  Figure  H.8.  The  third  coordinate  system  is  the  geographic 
coordinate  system  where  (i  j,k)  correspond  to  (East,  North,  Up). 


A  vector  in  the  ship  coordinates  can  be  related  to  the  geovertical  coordinates  by  a 
transformation  matrix.  Let  f  =  (x,y,z)  be  a  vector  in  the  tilted  coordinates  and 
r  =  (x,y,z)  be  the  corresponding  vector  in  the  geovertical  coordinates.  The  two  vectors 
are  related  by 


r 


cosyu  cosy21  cosy31 
cosy12  cosy22  cosy32 
cos  yt3  cos  y^  cos  y33 


f 


=  Mr 


(E.21) 


The  cosine  angles  ,  y y,  are  the  angles  that  the  tilted  unit  vectors  i,  j,  k  make  with  the 
geovertical  unit  vectors  i,  j,  k .  It  can  be  shown  (Karamcheti,  1967)  that  if  each  of  the 
coordinate  systems  is  orthogonal,  then 
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£cosyvcosy,  =  1 


i= i 


=  0 


if  i  =  k 
if  i  *  k 


(E.22a) 


and 


3 

^cosy,;cosya=  1 


=  0 


if  j  =  k 
if  j*k 


(E.22b) 


In  order  to  relate  the  two  coordinate  systems,  three  angles  must  be  given.  This  is  often 
done  using  the  Euler  angles  (Teichmann,  1969)  where  the  rotation  is  carried  out  in  three 
steps.  First  a  rotation  around  the  k-axis  through  angle  0,  then  a  rotation  about  the  i-axis 
through  angle  8  and  finally  a  rotation  about  the  new  k-axis  through  angle  The 
resulting  transformation  matrix  is  defined  by 
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E.5.  Platform  Motion  Correction 

There  are  four  main  steps  to  the  tilt  correction  for  the  slant  velocities.  1) 
mechanical  alignment  errors  are  determined  using  the  translation  correction  method.  2) 
one-minute  averaged  slant  velocity  profiles  are  corrected  for  beam  translation  by 
interpolating  velocities  to  common  depth  bins.  3)  horizontal  velocities  are  estimated 
subtracting  the  back  to  back  beam  translation  corrected  slant  velocities.  4)  geographic 
velocities  are  found  by  rotating  the  horizontal  velocities  using  one  minute  averaged 
heading  data.  Included  in  this  section  is  a  discussion  of  the  tilt  rotation  correction  which 
is  important  in  vertical  velocity  estimates  and  may  sometimes  by  used  to  determine 
mechanical  alignment  errors. 

E.5.1.  Platform  Tilt  and  Accelerometer  measurements 

A  three  axis  accelerometer  package  was  mounted  on  the  hull  of  FLIP  and  used  as 
an  inclinometer.  The  accelerometers  were  sampled  at  10  Hz  and  averaged  over  1  minute. 
Tilts  are  estimated  by  assuming  that  the  horizontal  accelerations  measured  are  from  the 
earth’s  gravitational  field.  The  package  is  mounted  such  that  the  x-axis  accelerometer, 
points  starboard  and  the  y-axis  accelerometer  points  towards  the  bow.  Note  that  the 
beams  directions  are  rotated  45°  with  respect  to  the  accelerometer  (Figure  E.2).  The  tilt 
angles  are  estimated  by 

(  \  ( n\ 

a  =  sm  —  and  p  =  sin  — 

U )  U  J  (E.20) 

where  ax  and  ay  are  the  horizontal  accelerations,  a  and  [5  are  the  x  and  y  axis  tilt  angle 
estimates  and  g  is  the  local  gravitational  vector.  This  approximation  is  only  accurate  for 
low  frequency  fluctuations  in  acceleration.  In  studies  of  a  free  drifting  FLIP,  Rudnick 
(1964)  showed  that  in  the  surface  wave  frequency  band,  0.05-0.2  Hz,  most  of  the 
horizontal  accelerations  are  due  to  the  surface  wave  forcing  while  at  low'er  frequencies 
horizontal  acceleration  variance  decreased.  Equivalent  testing  has  not  been  performed  on 
a  moored  FLIP. 

During  SWAPP,  a  two  point  mooring  held  FLIP  in  position.  The  tethering  of 
FLIP  above  the  water  line  causes  a  tilting  moment  when  there  is  a  near  surface  current. 
The  strongest  near  surface  currents  occur  at  inertial  currents  and  semi-diurnal  frequencies 
which  cause  tilts  of  up  to  2.5°.  It  is  clear  that  for  these  low  frequencies,  the  measured 
accelerations  are  due  to  tilt  and  not  horizontal  accelerations.  At  what  frequency  is  the  tilt 
signal  masked?  Capacitance  wave  wire  data  from  the  starboard  boom  (Figure  E.7a)  show 
that  at  higher  frequencies,  surface  waves  dominate  the  signal  with  a  swell  peak  at  0.08  Hz 
and  a  typical  wind  wave  spectral  slope  of  0)  ‘5  (Philips,  1957).  At  frequencies  below  the 
swell,  there  is  very  little  contribution  of  sea  surface  height  due  to  surface  waves.  Instead 
the  tilting  of  FLIP  lifts  the  boom  up  and  down  causing  apparent  shifts  in  the  surface 
height.  Concurrent  x-axis  accelerometer  spectrum  (Figure  E.7b)  shows  the  same  spectral 
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(cos  <p  cos  yz-sin  0  cos  9  sin  yz)  (sin  0  cos  yz+cos  0  cos  9  sin  yz)  sin  yz  sin  6 
-(cos0sin  yz+sin0cos$cosyz)  (-sin0sin  yz+cos 0 cos 0 cos  yz)  cosyzsin# 
sin  0  sin  9  -  cos  0  sin#  cos# 


(E.23) 


The  angles  measured  by  the  3-axis  accelerometer  package  are  not  the  Euler 
angles.  The  horizontal  tilt  angles,  a  and  /?,  that  are  estimated  from  the  accelerometers  are 
the  angles  the  tilted  unit  vectors  i,  j  make_with  the  geovertical  horizontal  plane.  The 
fact  that  there  is  no  rotation  about  k  or  k  allows  the  elimination  of  the  third  angle 
needed  for  the  transformation  of  the  tilted  coordinates.  Therefore,  the  tilt  angles  are 
related  to  the  cosine  angles  by 

sin  a  =  cos(a  +  90°)  =  cos  y31  =  -  cos  y13 

sin/?  —  cos(/?  +  90  )  —  cos  y33  —  —  cos  y^.  ^  24) 

Using  (E.22)  and  (E.23),  the  resulting  transformation  matrix  can  be  found  in  terms  of  a 
and  /3 


Ma„  = 


sin2  p  +  sin2  a(l  -  K)*  sin  a  sin  p ((1  -  ff)K  - 1) 

K  K 

sin  a  sin  /?((!  -  k)*  - 1)  sin2  a  +  sin2  /?( 1  -  k)* 


K 

-sin  a 


K 

-sin/? 


sina 

sin  P 
(1-K-)* 


(E.25) 


where 


k  =  sin2  a  +  sin2/?. 


Thus  a  vector  in  the  tilted  reference  frame,  r  is  related  to  the  corresponding  vector  in  the 
geovertical  reference  frame,  r,  by 


r  =  Ma^  . 


(E.26) 


It  is  assumed  that  there  is  some  mechanical  misalignment  between  the  accelerometer 
package  mounted  on  the  hull  of  FLIP  and  the  161  kHz  sonar  mounted  about  20  meters 
above  on  FLIP'S  thruster.  The  Coordinate  system  of  the  sonar  is  defined  by  the  measured 
tilt  angles  plus  some  constant  tilt  offset  such  that 


a'  =  a  +  ea 
F  =  P  +  e, 


(E.27) 
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The  design  tolerances  in  the  mounting  of  the  two  sensors  could  result  in  misalignment 
errors  of  a  few  degrees.  The  beam  directions,  however,  are  built  to  much  tighter 
tolerances.  The  sonar  housing  was  machined  from  one  piece  of  aluminum  and  the 
transducers  are  mounted  rigidly  in  place.  Orthogonality  of  the  beams  and  the  look  down 
angles  are  fixed  to  within  a  fraction  of  a  degree. 


E.5.3.  Rotation  Correction 


Plueddemann  (1987)  showed  that  the  rotation  correction  is  only  important  for  the 
vertical  velocity  estimates.  Since  the  goal  here  is  to  study  the  horizontal  velocities,  the 
rotation  correction  is  not  directly  used  in  the  final  velocity  calculations  but  is  presented 
here  for  completeness. 

The  rotation  correction  transforms  the  velocity  vector  into  geovertical  coordinates 
using  (E.25).  A  velocity  vector,  V  =  (u,v,w),  in  the  tilted  reference  frame  is  related  to 
the  corresponding  vector  in  the  geovertical  reference  frame,  V  =  (u,  v,w) ,  by 


V  =  MapV 


(E.28) 


This  correction  is  complicated  by  the  fact  that  the  sonar  only  estimates  the  along 
beam  projection  of  the  velocity  field  in  the  tilted  coordinate  system.  A  simple  reparation 
is  to  use  the  uncorrected  Janus  velocity  components.  As  shown  before,  the  Janus 
velocities  (for  slant  velocities  positive  towards  FLIP)  in  this  case  are 


u 

02  2cos0o 

(E.29a) 

w„  =  (V»  +  V2) 

2  sin  0O 

(E.29b) 

l  o 
<■> 

C-  r  ) 
II 

=r 

(E.29c) 

W  =(V>+V3> 

13  2sin0o 

(E.29d) 

where  60  is  the  downward  beam  angle.  The  sonar  is  rotated 
package  (Figure  E.2)  and  the  estimated  horizontal  velocities 
the  accelerometer  package  are  found  by 

45°  from  the  accelerometer 
in  the  coordinate  system  of 

-  _  V2 
u  —  «  u°2  * 

(E.30a) 

.  V2 

V  =  —  <U02  +U13) 

(E.30b) 
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Substitution  of  (E.30)  into  (E.28)  and  solving  for  the  vertical  velocity  component  in  the 
geovertical  coordinate  system  yields 

=  -[usina  + v sin /?)  +  (!  -  sin5  a- sin5 /?)Kw02  (E.3 1  a) 

wi3  =  ~(u sin  a  +  v  sin /?]  +  (1  -  sin5  a  -  sin5  /?)^  w13 .  (E.3 1  b) 

E.5.4.  Translation  Correction 

Let  B  be  the  unit  vector  along  a  beam  in  the  rotated  reference  frame.  Then,  by 
using  the  dot  product,  the  vertical  coordinate  in  the  tilted  reference  frame  at  range  r  along 
B  is  written  as 

5(0  =  r[B-kJ  (E.32) 

The  vertical  coordinate  in  the  geovertical  reference  frame  is  found  by  using  the 
transformation  matrix  and  is  written  as 

Hr)  =  r[(Ma/JB)-k].  <E33> 

Using  the  downward  beam  angle,  6q  =  60°,  the  unit  vectors  in  the  four  beam  directions 
are 


(E.34a) 

(E.34b) 

(E.34c) 

b3  =  -^^j-#k 

(E.34d) 

Using  equation  (E.26)  to  find  the  geovertical  vertical  coordinate  from  each  of  the  beam 
vectors  as  a  function  of  range  and  tilt  yields 


z0(r,a,/?)  =  r 
zx(r,a,p)  =  r 
z2(r,a,p)  =  r\ 
23(r,a,/J)=  rj 


+  ^-sina  +  ~sinp~^-{l-sin2  a -sin2/?)3 


V2_; 

4 


Vj, 

2 


4 

,  V2  /j  V3  . 


-^•sina  +  —sin/?  -  li(l  -  sin2  a  -  sin  ,9) 2 


-  ^  sin  a sin/?- ^(1-  sin2  a -sin2/?) 


V2 


V5  , 
2 


+  —  sin  a  -  ~  sin/?  ~  ^  (1  -  sin2  a  -  sin2  /?) 


(E.35) 


As  an  example,  let  r  =  200m  and  assume  the  tilt  angles  are  both  0.0°.  The 
corresponding  vertical  coordinate  for  each  beam  is  z  —  -173.2  meters.  Now  assume  the 
tilt  angles  are  both  1.0°  so  that  all  of  the  tilt  is  in  the  plane  of  beams  1  and  3.  Beams  1 
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and  3  will  both  have  the  same  vertical  coordinate  and  a  very  small  shift  to  shallower 
depths.  Thus,  solving  for  the  vertical  coordinates  yields  i\  =zt  -  -173.3  meters.  Beams  0 
and  2,  however,  have  different  vertical  coordinates  and  a  significant  shift  in  depth.  In  this 
case,  solving  for  the  vertical  coordinates  yields  zo=  -168.2  meters  and  zi  =  -178.1  meters. 

In  the  example  above,  Janus  velocities  calculated  without  tilt  translation 
correction  back  to  back  slant  velocities  that  are  separated  by  10  meters  in  depth. 
Since  this  near  the  resolution  of  the  velocity  estimates,  significant  errors  will  occur  if 
there  is  any  vertical  shear  in  the  velocity  field.  Besides  Janus  velocity  errors,  the  absolute 
depth  of  the  velocity  estimates  will  be  in  error.  This  becomes  important  when  velocity 
estimates  are  compared  with  density  data  from  the  profiling  CTD.  The  tilt  translation 
correction  is  a  necessary  step  to  the  velocity  data  processing. 

The  algorithm  for  implementing  translation  correction  proceeds  as  follows.  A 
standard  range  bin  vertical  coordinate  is  set  by  the  case  of  no  tilt.  Then  1  minute 
averaged  accelerations  are  used  to  calculate  an  average  tilt.  New  vertical  coordinates 
based  on  (E.35)  are  calculated  for  each  beam  using  the  average  tilt  plus  the  mechanical 
alignment  correction.  Finally  the  slam  velocities  from  each  beam  are  interpolated  to  yield 
values  at  the  standard  range  bins.  (An  example  of  corrected  and  uncorrected  slant 
velocity  profiles  is  shown  in  Figure  (E.  10).)  At  this  point,  Janus  velocity  estimates  can 
be  made  from  the  slant  velocities  at  the  standard  range  bin  depths. 

E.S.5.  Mechanical  Alignment  Error 

The  errors  in  mechanical  alignment  between  the  161  kHz  Doppler  sonar  and  the 
accelerometer  package  are  found  using  the  translation  correction.  The  method  used  is  a 
straightforward  search  for  the  offset  angles  that  result  in  the  highest  vertical  correlation 
between  back  to  back  slant  velocity  calculations.  The  underlying  assumption  is  that  the 
horizontal  velocity  vertical  finestructure  has  horizontal  scales  which  are  much  larger  than 
the  horizontal  separation  between  the  sonar  beams.  If  this  is  the  case,  then  shear  layers 
should  lie  on  the  same  vertical  coordinate  from  back  to  back  beams  and  there  will  be  a 
strong  velocity  correlation  between  the  beams.  If  the  beams  are  tilted,  however,  then  the 
translation  of  the  velocity  estimates  will  shift  the  shear  layers  to  different  depths.  This 
will  tend  to  decrease  the  correlation  between  the  back  to  back  beams. 

The  data  set  used  to  find  the  mechanical  alignment  errors  consisted  of  990 
consecutive  slant  velocity  profiles  and  tilt  angles  averaged  over  one  minute  from  yeardav 
72.  The  range  bins  25  to  75  were  used  which  correspond  to  depths  of  approximately  100 
to  245  meters. 

The  tilt  correction  angles,  £a  and  £p ,  are  allowed  to  vary  over  a  range  of  possible 
values.  For  each  combination  of  values,  the  translation  correction  is  applied  to  each  beam 
and  33  independent  30  minute  averaged  profiles  are  computed.  The  correlation  between 
back  to  back,  30  minute  averaged  profiles  is  found  from 
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This  correlation  is  divided  by  the  correlation  between  the  same  profiles  with  no  tilt 
correction  applied  to  create  a  correlation  ratio.  This  correlation  ratio  will  be  greater  than 
one  if  the  translation  correction  has  enhanced  the  corresponding  vertical  coordinate 
estimates  and  less  then  one  if  the  vertical  coordinate  estimates  are  corrupted. 

The  average  correlation  ratio  is  found  by  averaging  over  the  correlation  ratio 
obtained  from  all  33  profiles  and  both  pairs  of  back  to  back  beams  (Figure  E.l  1).  The 
mechanical  alignment  correction  angles,  £aand  £p,  are  varied  from  -2.25°  to  2.25°.  The 
highest  average  correlation  ratio  is  1.022  and  is  found  at  £a=-1.55°  and  £p  =-0.10°.  These 
are  the  mechanical  misalignment  correction  angles  used  in  the  tilt  correction  processing. 

An  investigation  is  done  to  see  if  the  rotation  correction  (E.28)  can  be  applied  to 
determine  the  mechanical  alignment  errors  as  well.  This  determination  is  based  on  the 
assumption  that  the  depth-time  average  of  the  true  vertical  velocity  is  small  compared  to 
the  magnitude  of  the  errors.  Minimization  of  the  errors  due  to  the  translation  of  the  sonar 
beams  is  accomplished  by  using  depth  averaged  slant  velocity  estimates.  Then  time 
averaging  removes  the  vertical  velocities  associated  with  high  frequency  internal  waves. 
The  remaining  velocities  are  assumed  to  be  due  to  horizontal  currents  measured  in  the 
tilted  coordinate  system.  Mechanical  misalignments  are  determined  by  minimizing  the 
mean  and  the  variance  of  the  estimated  depth-time  averaged  vertical  velocities. 

The  data  set  for  the  misalignment  investigation  consists  of  slant  velocity  and  tilt 
data  averaged  over  1  minute  in  time  for  14  days  from  yearday  63  to  77.  The  slant 
velocity  data  were  also  averaged  over  168  meters  from  100  to  268  meters  in  depth.  The 
uncorrected  vertical  velocities  and  the  horizontal  velocities  are  calculated  using  (E.29). 
The  rotated  horizontal  velocities  are  calculated  using  (E.30).  The  resulting  velocity 
estimates  axe  then  low  pass  filtered  using  a  5-pole  Butterworth  filter  with  a  cutoff  period 
of  20  hours  and  plotted  in  Figure  (E.9a-c).  The  horizontal  velocities  have  maximum 
amplitudes  of  5  cm/sec.  There  is  also  a  low  frequency  signal  in  the  horizontal  velocity 
data  due  to  a  from  moving  through  which  reverses  the  mean  current  almost  180°  from 
day  63  to  day  77.  The  tilt  angles  have  maximum  amplitudes  of  E.2°  The  tilting  of  FLIP 
at  low  frequencies  is  driven  by  near  surface  horizontal  currents,  near  surface  winds, 
ballasting  and  mooring  line  tension.  It  is  apparent  that  much  of  the  variance  in  both  the 
horizontal  velocities  and  the  tilt  angles  is  in  the  near  inertial  band  (-21  hour  period). 

Since  the  vertical  velocity  correction  is  a  function  of  both  horizontal  velocity  and 
tilt  angles,  one  would  expect  some  variance  in  the  uncorrected  vertical  velocities  in  the 
near  inertial  band.  Theoretically,  inertial  waves  have  no  associated  true  vertical  velocities 
so  this  frequency  band  might  be  a  suitable  frequency  to  use  in  the  determination  of 
mechanical  alignment  errors.  However,  close  examination  of  the  low  pass  filtered 
vertical  velocities  (Figure  E.9c)  shows  that  is  not  the  case.  A  diurnal  signal  (24  hour 
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period)  with  a  velocity  amplitude  of  0.4  cm/sec  dominates  the  uncorrected  vertical 
velocity  signal  about  a  mean  velocity  of  -0.3  cm/sec.  Doppler  shifts  at  this  frequency  are 
dominated  by  the  nightly  migration  of  plankton  towards  the  surface  to  feed.  This  signal 
dominates  over  any  rotation  errors.  (For  example,  a  vertical  velocity  correction  of 
magnitude  0.3  cm/sec  needs  a  horizontal  current  of  4  cm/sec  and  a  tilt  of  4.3°.)  An 
analysis  of  lower  frequencies  may  not  be  useful  since  there  may  be  strong  daily  variations 
in  the  migration  of  the  plankton  which  is  not  correlated  with  the  horizontal  velocities  and 
tilt  angles.  The  rotation  correction  will  not  be  used  in  the  determination  of  mechanical 
alignment  errors  between  the  161  kHz  sonar  and  the  accelerometer  package.  The  strength 
of  the  diurnal  plankton  migration  masks  any  estimation  of  vertical  velocity  errors. 

E.5.6.  Heading  Correction 

The  heading  of  FLIP  was  constrained  by  a  2  point  mooring  which  held  the  bow 
pointing  between  210°  and  275°  TN.  Wind,  upper  ocean  currents  and  mooring 
retensioning  all  dynamically  affected  the  heading.  The  heading  correction  consists  of 
rotating  the  tilt  corrected  horizontal  velocities  to  geographic  coordinates.  The  geographic 
coordinate  system  is  defined  such  that  k  points  vertically  upward  parallel  to  gravity,  j 
points  horizontally  towards  the  North  and  \  points  horizontally  towards  the  East.  The 
rotation  of  the  horizontal  velocities  is  obtained  by 

=  u13cos(0-315°)-uO2  sin(#-315°) 
vf  =  u,3sin(0-315°)  +  uO2Cos(0-3150)  ^  ^ 

where  9  is  the  compass  heading,  u,  is  the  east  velocity  component,  v,  is  the  north  velocity 
components  and  u„  and  u„  are  the  horizontal  Janus  velocities.  FLIP'S  compass  supplied 
the  heading  data  at  a  sample  rate  of  10  Hz.  This  heading  is  then  averaged  over  one 
minute  and  the  heading  correction  is  applied  to  one  minute  averaged,  tilt  corrected 
horizontal  velocities. 
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a)  4  repeats  of  a  7-bit  Barker  code 
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Figure  E.l  Repeal  Sequence  Coded  Transmit  Pulse,  a)  Schematic  Representation  of  a  repeat  sequence 
code.  The  transmitted  waveform  consists  of  a  seven  bit  Barker  code  repeated  four  limes,  b)  The 
corresponding  autocorrelation  function.  Values  are  small  except  at  even  multiples  of  the  subcode  length, 
c)  The  frequency  spectrum  of  the  code.  The  spectrum  of  an  uncoded  sinusoidal  pulse  of  the  same  duration 
would  consist  of  a  single  peak  of  width  equal  to  the  width  of  any  individual  peak  in  the  above  spectrum. 
(Pinkel  and  Smith,  1992) 
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Table  E.l  MPL  161  kHz  Incoherent  Doppler  Sonar  System  Configuration  Parameters .  Presented  m 
this  table  are  the  Doppler  Sonar  system  parameters  as  configured  for  the  SWAPP  deployment. 


MPL  161  kHz  Incoherent  Doppler  Sonar  System  Configuration  Parameters 


Transmit  Frequency 

fc 

161  kHz 

Sample  interval  (bit  length) 

X 

250  msec 

Bandwidth 

_  Af 

4000  Hz 

Sequence  Length 

.75  sec 

Sequence  Averaging 

P 

80  seq/estimate 

Range  bin  length, 

number  of  samples  in  covariance  estimate 

N 

9  bits 

Subcode 

G(0 

7  Bit  Barker  code 

(1,1,1, -1, -1.1, -1) 

Subcode  Length 

L 

7  bits 

Number  of  Repeats  of  Subcode 

M 

6  repeats 

Averaging  Length 

T„,  =  Nr 

.00225  sec 

Overlapped  common  range  in  the  covariance 
calculation 

T„i  -  (M  -\)Lz 

.00875  sec 

Sound  Velocity 

c 

1485  m /sec 

Maximum  Unambiguous  Velocity 

V  =  c 
m“  4/rL 

131.7  cm/sec 

Averaging  Length 

=  maxCT^,!^) 

.00875  sec 

Effective  Range  Resolution 

T 

A R=c  ro" 

2 

6.50  meters 

RMS  velocity  imprecision 

AV 

0.74  cm/sec 

Beam  Down  look  Angle 

So 

60° 
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Figure  E.2  161  kHz  Beam  Directions  and  Accelerometer  Tilt  Axis  Directions.  Looking  down  from 
above  (with  the  z-axis  coming  out  of  the  page),  this  diagram  shows  the  direction  of  the  four  161  kHz 
AD  CP  sonar  beams.  The  ADCP  was  mounted  on  FLIP'S  direction  thruster.  FLIP  was  held  in  place  by  a 
two-point  mooring  and  the  thruster  was  never  used  during  SWAPP.  Q  is  the  compass  heading  from  FLIP 
and  varied  from  215°  to  270°.  The  three-axis  accelerometer  package  was  mounted  farther  down  the  hull 
and  was  rotated  45°  with  respect  to  the  sonar  beams.  The  x-axis  pointed  towards  starboard  and  the  y-axis 
pointed  towards  the  bow  in  the  same  direction  as  the  heading  vector. 
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Figure  E.3  Variable  Gain  Correction.  To  extend  the  dynamic  range  of  the  analog  to  digital 
conveners,  a  variable  gain  circuit  is  used.  The  gain  was  increased  30  dB  with  60  steps  of  0.5  dB  each. 
The  first  step  occured  at  sample  number  60  (15000  msec)  and  lasted  20  samples  (5000  msec). 
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Figure  E.4  1  Hour  Averaged  Intensity  Profile.  1  Hour  of  raw  intensity  data  from  beam  0  is  averaged 
and  plotted.  Also  plotted  is  the  variable  gain  corrected  intensity  data.  The  gain  corrected  profile  shows 
nearly  80  dB  of  dynamic  range, 
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a)  G(t)  Variable  Gain 


Figure  E.5  Hour  averaged  Intensity  and  Variable  Gain  First  Differenced  in  Range.  The  variable  gain 
ramp  consists  of  60  gain  steps  which  become  more  apparent  in  when  first  differencing  is  applied.  The 
range  bin  averaged  gain  is  plotted,  (a),  and  shows  a  ripple  during  the  gain  ramp  (bins  7  to  122).  This 
ripple  can  also  be  seen  in  the  1  hour  averaged  intensity  data,  (b),  during  the  gain  ramp.  When  the 
variable  gain  correction  is  used,  most  of  the  ripple  due  to  the  gain  steps  is  removed,  (c). 
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Figure  E.6  1  Hour  Averaged  Covariance  Magnitude  Squared.  One  Hour  of  1  minute  averaged 
covariance  data  from  beam  0  is  averaged  and  plotted.  Also  plowed  is  the  variable  gain  corrected  covariance 
data.  The  gain  corrected  profile  shows  nearly  80  dB  of  dynamic  range. 
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Figure  E.8  Geophysical  and  titled  coordinate  systems.  The  geovertical  coordinate  system,  1.  j.  k,  is 
defined  with  the  vertical  coordinate  parallel  to  gravity.  The  tilted  coodinate  system,  '•  }•  k ,  is  one  that 
rotates  with  the  movements  of  FLIP  and  is  defined  by  the  tilt  angles  a  and  b.  The  tilt  angles  are  infered 
from  horizontal  acceleration  measurements. 
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Figure  E.9  Low  Pass  Filtered  Tilt  Angles  and  Depth  Averaged  Velocities.  Slant  velocities  are  depth 
averaged  from  100  to  275  meters  then  Janus  velocities  are  calculated.  The  data  has  then  been  low  pass 
filter  with  a  cutoff  period  of  20  hours.  The  estimated  horizontal  velocities  in  the  coordinate  system  of  the 
accelerometer  package  ( u  solid;  v  dashed),  are  shown  in  (a).  The  tills  (a  solid;  b  dashed),  are  low  pass 
filtered  and  shown  in  (b).  Variance  in  the  horizontal  velocities  and  tilts  is  mostly  in  the  inertial  (21 
hours)  band.  The  estimated  depth  averaged  vertical  velocities  ( wn  solid;  w„  dashed)  are  shown  in  (c). 
The  dominated  signal  in  the  vertical  velocities  is  in  the  diurnal  (24  hour)  band  due  to  the  vertical 
migration  of  the  acoustic  scatterers  (plankton). 


a)  Unconcctcd  slant  velocities  b)  Translation  corrected  slant  velocities 
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Figure  E.10  Translation  Corrected  Slant  Velocities.  Slam  velocity  profiles  and  till  angles  arc  averaged  over  1  minute.  60  profiles  (1 
hour)  arc  averaged  together  before  plotting  versus  the  vertical  coordinate.  The  average  tilts  plus  the  ncchantcal  alignment  adjustments  arc 
a=-1.124°  and  (1=0.189°.  No  tilt  correction  has  been  applied  to  the  profiles  in  (a).  The  tilt  correction  is  applied  to  every  profile  before  the 
one  hour  average  is  calculated.  The  lilt  corrected  profiles  arc  shown  in  (b).  Note  that  tiic  vertical  structure  in  die  tilt  corrected  back  to  back 
beams  is  more  correlated  in  depth  than  those  without  tilt  correction. 
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Figure  E.ll  Contour  of  Correlation  Ratio  vs.  Mechanical  Misalignment  Angles.  The  correlation 
ratio  is  tilt  corrected,  back  to  back  beam,  slam  velocity  correlation  divided  by  the  non  corrected  slant 
velocity  correlation.  This  correlation  ratio  is  calculated  over  a  range  of  possible  mechanical  misalignment 
angles.  Contour  of  correlation  ratio  is  plotted  with  contour  values  starting  at  1.000  and  contour  intervals 
of  0.002.  The  maximum  correlation  ratio  is  1.022  and  is  located  at  Eq  =  -1.55°  and  ep  =  -0.10°.  These 
are  the  alignment  correction  angles  used  in  the  translation  correction. 
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